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FRgFACS 

The  Department  of  Aeronautics,  Naval  Postgraduate  School 
offers  several  courses  at  different  levels  on  the  dynamics  of 
flight  vehicles.  These  courses  include  studies  of  fixed  wing 
stability,  rotary  wing  stability,  and  the  stability  of 
missiles.  In  addition,  there  are  general  courses  that  cover 
the  dynamics  of  bodies  in  six  degrees  of  freedom,  including 
the  development  of  the  equations  of  motion. 

There  is  currently  no  signle  resource  available  that 
covers  the  material  for  the  airplane  and  missile  courses,  as 
well  as  the  introductory  dynamics  classes,  and,  as  a  result, 
the  instructors  must  make  use  of  large  numbers  of  locally 
prepared  handouts  and/or  assign  several  textbooks,  only 
portions  of  which  will  be  used. 

In  order  to  provide  a  solution  to  this  problem,  the 
following  material  has  been  developed  for  use  in  the  dynamics 
classes  in  the  Department  of  Aeronautics.  Although  this 
material  is,  in  places,  of  greater  depth  than  what  would  be 
covered  in  the  courses,  it  is  felt  that  the  value  of  this 
material  is  enhanced  by  its  degree  of  thoroughness. 
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CHAPTER  1 
DYNAMICS  OF  MOTION 

1.1  INTRODUCTION 

This  chapter  is  intended  to  be  a  formal  review  of 
mathematical  and  physical  concepts  that  will  be  used  to 
develop   the   theory  of   dynamic  stability  of -flight  vehicles. 

The  starting  point  is  the  basic  laws  of  Newtonian 
Mechanics,  with  a  useful  application  to  spring-mass-damper 
systems.  We  shall  find  later  how  the  motion  of  any 
conventional   flight   vehicle   is   similar   to  those    systems. 

We  then  proceed  to  rigid  body  equations  of  motion, 
described   in  non-inertial    reference   frames. 

1.2  NEWTON'S    LAWS    OF   MOTION 

In  any  system  of  units,  basic  quantities  and  units  are 
related  through  Newton's   laws   of   motion. 

NEWTON'S    1ST    LAW   OR   THE    LAW   OF   COORDINATE    SYSTEMS: 

If  there  is  no  force  acting  on  a  particle,  its  velocity 
remains  constant. 

This  law  is  often  called  the  law  of  inertia  and  applies 
to   inertial    frames   of    reference. 

NEWTON'S    2ND   LAW   OR   THE    LAW   OF   DYNAMICS: 

The  time  rate  of  change  of  the  linear  momentum  of  a 
particle,    is   proportional   to  the   net  force   acting   on   it. 
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The  vector  quantity 

^(4)  =  ^W  ufi)  (1-1) 

is  defined  as  the  linear  momentum  of  the  particle,  where  m(t) 
is  the  inertial  mass  of  the  particle  and  v(t)  is  its  velocity 
relative  to  the  reference  frame  in  question. 

The  2nd  law  is  mathematically  expressed  as 


(1-2) 


or 

If  the  mass  m  remains  constant,  we  can  write: 


(1-3) 


F.el    =1   ^j^vLi)  ^^_^^ 


or 


""^^       3c  (1-5) 

where   a(t)     is    the    acceleration   caused    by    the    net   applied 
force     f^'wet. 

The  proportionality  factor  is  expressed  as  V^c  and  its 
value   depends   on  the   system   of    units   being   used. 

In  the  English  Engineering  system,  one  pound  force 
(1-lbf)  will  cause  to  one  pound  mass  (1-lbm) ,  an  acceleration 
of   32.174    ft/sec^      Thus   substituting   into  Equation  1-5 

i-ih^    -  ~  Clbrn)    32.174    C-rVsec^J  ^^_^^ 
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Hence 

9,  =  32.174        (brv.    ^t/lbf  sec^  (l_7j 

A  few  other  units  used  for  velocity,  are  given  below  with 
their  corresponding  conversion  factors: 

1   mile  per   hour    (Imph)    =  1,4666   ft/sec 
1   knot  =   1   nautical   mile  per   hour   =  1.6878   ft/sec 

NEWTON'S    3RD   LAW   OR  THE    LAW   OF   MUTUAL    FORCES: 

If  one  body  exerts  a  force  on  a  second  body,  the  latter 
exerts   an   equal   but   opposite   force   on  the   first. 

This  law  allows  us  to  define  the  gravity  force,  or  the 
weight  of  a  body,  that  has  a  given  mass  (m) .  It  is  the  force 
that   results   due   to  the   earth's   attraction. 

Besides  the  three  fundamental  laws  of  mechanics,  the 
Newton's  Universal  Gravitational  Law  states  that  two  inertial 
masses  ml  and  m2 ,  separated  by  a  distance  R  are  subject  to  an 
attraction  force   according   to  the   equation: 

*■   ^  9=^         H"-  (1-8) 

where    G,     is    the    gravitational    constant,    which    in    the   EE 
system   is       G  =  S,52x-io""     'bf --ftVlbw^ 


25 

If    now   m i  =    Mg    =    1.32    X    10         Ibm    (mass   of    the    earth), 
then    near    earth's    surface    where   R   =    Re    =    2.09    X    10      ft 


(radius  of   the   earth),   we   can  write: 


W--3c   -W-   "^  (1-9) 
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or 

W-'YIO^  (1-10) 

where   g    =   g^GMe/Rg    =   32.174    ft/sec^,    called   acceleration  of 
gravity. 

By  Equation  1-10  we  can  see  that  1-lbra,  subject  to  the 
earth's  gravitational  field,  experiences  a  force  1-lbf  and 
hence  it  accelerates  with  32.174  ft/sec^,  near  earth's 
surface. 

We  can  also  define  another  unit  of  mass,  the  slug,  as  the 
mass  that  experiences  an  acceleration  1-ft/sec^,  if  subject 
to  a  force   of   1-lbf. 

By  this  definition 

1    slug   =   32.174    Ibm  (1-11) 

Consequently,  we  can  derive  the  density  (p) ,  as  the  mass 
per  unit   volume   usually  measured   in   slugs/ft. 

ddSiry   Cp)  = 


voluwe  (.v)  (1-12) 

People  in  the  field  use  the  (g) ,  as  a  short  hand  notation 
for  force  or  acceleration.  A  5-g  force  simply  means  a  force 
5    times   the  body's  weight. 

1.3      ANGULAR   MOMENTUM 

A  special  significance  to  the  2nd  law  is  given  by  the 
angular  momentum  approach  derived  simply  by  crossing  the 
radius  vector    r   ,    on   to  the   2nd   law. 
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r  X  — -    =  r  X    r,r)ei 

cH^  (1-13) 

We  define 

G-wei    =    r  X  F^e^  (1-14) 

the    net    torque    or    moment    of    the    force   about   the   origin   and 

h  =  rxp  (1-15) 

the  angular  momentum  about  the  origin. 

Differentiating  Equation  1-15  with  respect  to  time,  we 

get 


dh 


dt    di  ^ 


The  first  term  is  zero,  because 


(1-16) 


dr      -»        -•,      -.^  r  •*    ■*  \     ^ 

—  xp    =^)xC1n^»}    =  ^yri  iv  yv  }  ^=0 

at 


(1-17) 


oriQin 


Figure  1-1 
Moment  of  Force  About  the  Origin 


So  we  finally  get 


dk 
dt 


-  G  ^et 
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(1-18) 


which  tells  us  that  the  time  rate  of  change  of  the  angular 
momentum  of  a  particle,  equals  to  the  net  torque  or  moment 
acting  on  it. 

If  no  torque  is  acting  on  a  particle,  its  angular 
momentum  is  conserved,  likewise  its  linear  momentum  is 
conserved,  if  no  force  is  acting  on  it. 

A  simple  example 

For  a  pendulum  that  swings  in  a  vertical  plane,  having  a 
mass  (m)  ,  the  acting  forces  are,  its  weight  W  =  mg  and  the 
tension  of  the  rod  T.  The  rod  is  assumed  to  have  negligible 
mass. 


N  T 


Figure  1-2 
A  Simple  Pendulum 

If  9  is  the  angle  between  the  rod  and  the  vertical  axis, 
s,  is  the  arc  length  corresponding  to  the  angle  0  and  1,  is 
the  length  of  the  rod,  a  direct  application  to  the  2nd  law 
(Equation  1-5),  gives: 
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--mgsin©    =  yf[s  (1-19) 

and   equilibrium  condition  in   the  direction  of    the   rod,    gives: 

-'YTjacosS  =  T  (1-20) 

Equation  1-19,    can   be  written 

S  +gsin0   =0  (1-21) 

If    9   is   small,   we   can   assume 

0  ^  f  ^  5  =  29 

^  (1-22) 

and 

0  =    sine  (1-23) 

With   this   assumption.    Equation  1-21,    can   be  written 

^  (1-24) 

Equation   1-24,    is    a    linear    differential   equation,   with 

constant    coefficients,    which   describes    the   motion   of    the 

pendulum. 

The  same  result  can  be  derived,  if  we  use  the  angular 

momentum  approach. 

The  angular  momentum  of  the  pendulum,  about  the  pivot,  is 

simply  (?),  cross  the  linear  momentum  (tds)  ,   The  time  rate  of 

change  of  the  angular  momentum,  has  to  be  equal  to  the  net 

external  torque, 

di  ^      ^       5  (1-25) 

Since  i    is  constant  (df/dt  =  0) ,  we  can  write 

£ws  =  -  ^'Knjsm©  (1-26) 
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or 


or 


s  +  iY)a  sin9  =  O 


e.fe-o 


(1-27) 


(1-28) 

The  minus  sign  in  Equations  1-19  and  1-25,  indicates  that 
if  the  angle  9  decreases,  the  acceleration  i  in  Equation  1-19 
and  the  angular  momentum  in  Equation  1-25,  increase. 

1.4   UNDAMPED  HARMONIC  OSCILLATOR 

In  this  section  and  in  the  following  section,  we  will 
derive  the  equations  and  the  solutions  to  the  spring  mass 
damper  systems. 

The  undamped  oscillator  is  simply  a  spring  mass  system 
without  friction,  shown  in  Figure  1-3. 

Let  k,  denote  the  spring  stiffness,  m,  the  mass  of  the 
body,  x(t),  the  displacement  of  the  body  from  the  equilibrium 
position.  Notice  that  the  weight  of  the  body  stretches  the 
spring  permanently,  so  that  Xe  is  the  equilibrium  position, 
after  the  body  has  hanged  to  the  spring. 


Xe 


^) 


LUitmi 


! 


/VM 


P, 


Figure  1-3 
Spring  Mass  System 
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The   restoring   force   according   to  Hook's   law,    is 

Fs  =  -kx  (1-29) 

Then  according  to  the   2nd   law 

'mx   =  -  kx  (1-30) 

X   +  —  X   =0 

^  (1-31) 

As  an   example  we  will  take 

m  =   18    Ibm  (1-32) 

k   =  40.5    Ibf/ft  (1-33) 

If  we   let 

^  k 

^'  /w  (1-34) 

where   (a^tj    is    the   so   called  natural    frequency  of   oscillations. 
Equation  1-31   becomes 

X  +  ULi«  X    =0  (1-35) 

This    is    a    linear    differential    equation  with   constant 
coefficients.      Its   characteristic     equation   is   clearly 

^^cj;    --0  (1-36) 

which  has   characteristic   roots 

^  =   ijOJ^  ■  (1-37) 

The    most    general    solution    of    Equation    1-35    can    be 
expressed  as 

.tt)  =  (Me^"'"^  WMe"^""'  ,,.33, 
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i 


where    A^.    and    A_     are    constants    depending    on    the    initial 
conditions. 

Expanding    the    exponential    terras    of    Equation   1-38,    we 

obtain 

xC-tJ  =  ex  cosoUnt  *-i3  s\nuJni  (1-39) 

where 

ex  =  A+  +  A_  (1-40) 

and 

^  =  S^^^'^-)  (1-41) 

If  we  now  express 

smuJ„i    =  cos  (lUnl  +  gO")  (1-42) 

we  can  give   the   following  useful   form  to  our    solution 

><(i)    =^   X-max    COS   CuJn-t^©-)  (1-43) 

where      x^gj<     is    the   maximum   deviation   from    the    equilibrium 
position  and   Oo  ,    the   initial   phase. 

We    can    relate    the    quantities     x^ox     and    ©»  r    with    the 
initial   conditions,    as   follows: 

At    time    t=0,    the    initial    displacement   x,  ,    is    given  by 
Equation  1-43 • 

Xo  =  X^ox   cos  Qo  (1-44) 

and    the    initial    velocity    x^  ,    is    given   by   differentiating 
Equation  1-43    with   respect   to  time   and   evaluating   at   t=0. 

^o  =  -x^ax  "^r,  sme,  (1-45) 

As  an   example  we  will   take 

Xo   =   4    ft  (1-46) 

Xo  =   32    ft/sec  (1-47) 

23 


From  Equations  1-44  and  1-45,  the  initial  phase  and  the 
maximum  displacement  can  be  found  to  be 


e.  =  w  ( 


-uJk)Xo   /  (1-48) 


Note:  In  Equation  1-48,  keep  the  sign,  so  that  the 
quantity  in  brackets  will  give  the  right  quadrant  for  the 
angle    ^o  . 

For   our   example 

2  K  405 

,024    ^^'2 


/  1024    \'^ 


(1-52) 


so   that  the   response    is 

xCl)  -    -5.4*3   C05(8.5U  ^2.59  rad)  ^1-53) 

1.5      DAMPED   HARMONIC   OSCILLATOR 

Consider  the  spring  mass  damper  system  shown  in 
Figure  1-4.  Assume  that  the  friction  force  exerted  by  the 
damper,    is   a   linear   function  of    the  mass  velocity. 
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Figure  1-4 
Spring  Mass  Damper   System 


Let  k  the  spring  stiffness,  m  the  mass  of  the  body,  x(t) 
the  displacement  of  the  body  from  the  equilibrium  position 
x^,  as  it  was  defined  in  section  (1.4),  and  c,  the 
coefficient   of    friction. 

The   restoring   force    is   again 

F3  =-kx  (1-54) 

and  the   friction  force    is 


or 


Then   according   to  the   2nd   law 

omx    =-cx-i<x 


x+—  x+— X   =0 


(1-55) 


(1-56) 


(1-57) 
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As 

an 

example  we  will 

take 

m  = 

18  Ibm 

k  = 

40.5  Ibf/ft 

f 

nc 

we 
J 

let 

CO.?  =gc 

c 

(1-58) 
(1-59) 

(1-60) 


(1-61) 

where    5    is    the    so    called    damping    coefficient,    then   the 
Equation  1-57    is  written 

W  +  2Sa;„  X +- oUti  X   =0  (1-62) 

This    is    a    linear    differential    equation  with    constant 
coefficients.      Its   characteristic   equation  is   clearly 

51'  +  25u;„:^  4-  a;/  =0  (1-63) 

Three    cases    arise     in    obtaining    solutions    to    the 
differential   Equation  1-62, 

As  an   example,   we  will  take 

^o    =   4    ft  (1-64) 

^o    =   32    ft/sec  (1-65) 

Case   1:      Critically  Damped   Motion    (^=i) 
As   an   example  we  will   take 

c   =  9.52     lbf/^|.scc-'  (1-66) 

so   that 

In    this    case,    the    characteristic    equation  has   a   double 
root 

r,=-a;^  (1-68) 
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and   the  raost  general   solution  is 

xCi)  ^  (A^^AA)  e'^"^  (1-69) 

where    A+    and    A_    are   constants    determined    by   the    initial 
conditions  as   follows 

Af   =  ^o  (1-70) 

and 

A-    =    ^o  -^  XoUJn  (1-71) 

so    that    the    solution    in    terms  of   the   initial   conditions  can 
be  written 

x(0  =[xo  (Ua}„^}  +  )^lj  e'*^"  (1-72) 

For   our   example 


xU)    r  [4  (  If  8.5i-t)  +  32i]  G 


-  8.5i  i 


or 


-8.su 


(1-73) 


(1-74) 


Casg  2; Overdamped  ^^Qtion   (^>1) 

As  an   example  we  will   take 

c   =   14.28   Ibf/ft-sec'^  (1-75) 

so   that 

^  "-^^     2maJr,     ^   ^^  (1-76) 

In  this   case   the    roots   of   the   characteristic   equation  are 
real   and   negative 

^    X   -^UD„   iCOr,  (5'-i)''^  (1-77) 

and   the  most  general    solution   is 


U)  =   G  ^4-  <2  +  A.  G 


(1-78) 
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where  A^  and  A_  are  determined  by  the  initial  conditions  as 
follows 

x„  =  A+  4- A.  '  (1-79) 


x»  =  -Sujr,  (A^.  +  A.)  4-(ju„y5^-l  fA4.-A.)  (1-80) 

which  yields 


(/5M  -BJy.-^^ 


2/5M^  (1-81) 

2/5^  (1-82) 
For  our  example 

A4.  =  6.36  ft  (1-83) 

A_  =  -2.36  ft  (1-84) 

and 

-12.77-t       951^.       -9.54  4:  ^ 
x(-l)  =  e     (636  (2    -  2.56  €     i  (1-85) 

or 

f,^  -3.2fii       -22.28 -t 

yC-t)  =  6.S6  (?     -  2.36  e  (1-86) 

Case  3; Underdamped  Motion  fg<l) 

As  an  example  we  will  take 

c  =  4.76  Ibf/ft-sec'^  (1-87) 

so  that 

"^"   2rr)U)o  "  ^'^  (1-88) 

In  this  case  the  roots  of  the  characteristic  equation  are 
imaginary  with  negative  real  parts. 
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If  we  define 

00^  ^COr,  /<-S'    ({or   5<1)  (1-90) 

the  so  called  damped  frequency  of  oscillations, 

^  =  -  ^(jOr)  ±  juJi  (1-91) 

and  the  most  general  solution  is 

xlU  =  e    [^^Q        +  A.  G    j  (1-92) 

Expanding  the  exponential  terms  we  obtain 

xlij  =  G    ((^cosGL)^t  f^  sinctyi)  (1-93) 

where 

ot  ^  A+  +  A-  (1-94) 

and 

^=j(A^.-A-j  (1_95) 

The  solution  can  be  formulated  as 


x{i)    =  X    G  CDsCcOjifGo)  (1-96) 

where   X    and     Go    are    constants    determined    by    the    initial 
conditions  as   follows:      From   Equation  1-69   we   obtain 

e^""""^  x(i)   =  X  cos  Ca;d-t+  Go)  (1-97) 

Differentiating   Equation  1-97    with   respect   to  time  we   obtain 

€  "^"     [5aj„  xa)f  xdj]   =  -XcUd  sin  ((!!,{  f0o)  (1-98) 

For    t=0.    Equations  1-97    and   1-98   become    respectively 

x^  =X   cos  0o  (1-99) 
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and 

5(ju„  x„  +  Xo  =  - X  GUj  sin  0e 
So  the   constants   are  determined  as 


(1-100) 


-tt)d  Xo 


(1-101) 


and 


X   =i[ 


Xo         + 


(^(jUnXo*-Xo)' 


W<^ 


i/2 


(1-102) 


(1-103) 
(1-104) 
(1-105) 


,,.  -4.26-t 

xU)=-f.?Ge  cos  C?.3?t  +  2.11  rad) 


Note:  In  Equation  1-101,  keep  the  sign,  so  that  the 
quantity  in  brackets  will  give  the  right  quadrant  for  the 
angle   9o  . 

For   our   example 

ujj    =  8.5 ^  (\'0.'b^f^    =  7.5?    sec'^ 
0o   -    121  02°    -  2.11    rad 
yi   =    ±  7.76  ft 
so   that  the   response    is 

26-fc 

(1-106) 

The  roots  of  the  characteristic  equation  can  be  plotted 
in  the  complex  plane  as  shown  in  Figure  1-5. 

In  cases  1  and  2,  the  damping  is  so  large  that  no 
oscillations  take  place  and  the  mass  simply  returns  gradually 
to  the  equilibrium  position  x=0.  In  those  cases  the  roots 
are  real  and  negative. 

In  case  3,  the  damping  has  been  reduced,  so  that 
oscillations  about  the  equilibrium  position,  do  take  place, 
the  amplitude  of  which  dies  out  as  t  ^oo  . 
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Hz 


•Figure  1-5 
Plot   of   the  Characteristic  Roots    in   the  Complex  Plane 


Figure  1-6  shows  the  responses  in  the  above  three 
described  cases  where  it  is  assumed  that  the  initial 
displacement   is  x^  ,    and   the   initial   velocity   is    ><o=0. 


critically   damped  moiion   S*i 

oycrckamped  mo^ton     5>1 

underdo mpeei  mohon    5<1 


envelope     C  e 


5uH,t 


Figure  1-6 
Response  of  Displacement  (x)  of  a  Spring  Mass  Damper  System 
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Notice  that  in  the  case  where  the  damping  5  is  zero,  we 
have  a  spring  mass  system,  without  friction,  described  in 
section  1.4.  In  this  case,  the  roots  are  imaginary  and  u)^  is 
the   same   as    w^,    (Figure  1-5)  • 

It  is  clear  that  the  system's  physical  constants  (in  our 
case  B,m,c,k),  uniquely  define  the  roots  of  the 
characteristic  equation.  Those  roots  are  given  the  name: 
normal  modes  because  one  can  determine  from  those  roots,  the 
corresponding  characteristic  motion  (in  our  case,  undamped 
critically  damped,    overdamped,    underdamped) . 

We  have  seen  so  far  how  an  harmonic  damped  oscillator 
responds,  if  initially  disturbed  from  its  equillibrium 
position.  Then  knowing  the  initial  conditions  we  can  always 
find  the  response  of  the  displacement  or  velocity  versus 
time. 

Now  we  are  interested  in  what  happens  if  the  oscillator 
is   driven  by  an   external    input    force   F(t). 

In   this   case    the   solution  will  consist   of   two  parts: 

M   =   X^,o^  (i)  ^  KpaK  (.^)  (1-107) 

where  ^homC^)  is  the  homogenous  part  of  the  solution 
discussed  so  far,  i.e.  considering  the  input  force  zero  and 
^pur-l  C-t)  is  the  particular  solution  due  to  this  particular 
applied   force,    i.e.    considering   the   initial   conditions  zero. 

To  obtain  particular  solutions  one  can  use  the  method  of 
undetermined  coefficients,  the  method  of  variation  of 
parameters   or    the   Laplace   Transform  method. 
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In    illustrating   the   above  methods  we   are  considering   the 
following   cases: 
(a)    Sinusoidal    Inputs 

The   spring  mass  damper    is   driven  bythe   sinusoidal    force 
F(t)    =  Fo    cos  CO   t  (1-108) 

As  an   example r   we  will  take 

F(t)    =   24    cos   3.5    t  (1-109) 

The    method    of    undetermined    coefficients    assumes    a 
superposition  of    sinusoidal   solutions  of   the   form 

^p(xr\^^)    =    A^.  coscdW  A_  sm  ujI  (1-110) 

or   equivalently 

XpQK^^)   -  X  cos(ujU9o)  (1-111) 

where  X  and    ©o    are   the   coef fieients   to  be   evaluated. 

Substituting   into  the  differential   equation 

X  •»■  25(jOy,x  +  (jj^  X    -    Fo  cos  cot  (1-112) 

we  obtain 

Cu;Tj'-aj')cosCaui+0o)-25cor)a)sin  (ujU0»)  =  ^  c05oot  (1-113) 

Expanding   cosine   and   sine   and  collecting   terms,   we  get 
[  (oUm -cu*]  cos  0o  -  25a)na)sin0o]  cosu}^ 

-[(ujTf-u;")  sm9o  +  25uJna)cos6o]  sincot  =  ~  cosiul:  (1-114) 

Equating   coefficients 

CuUt? -0)')  COS0O  -  25uJnUJsm  Go  =  r^  (1-115) 

and 

(W»,^-U)')   SinOo   +  2S(JUr,WC05Go    =0  (1-116) 
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From  Equation  1-116  we  obtain 


ian  Go   = 


-2SaJnUJ 


(1-117) 

Multiplying    Equation   1-115    by   cos6o.    Equation  1-116   by 
sin  ^ ,    and  by  adding,    we  obtain 


Ca)^2-u)^)    .  ^  COS0O 
The   cos  9o   is   expressed   in  terms  of    tan  ©o  as 

V-tan  6o 


(1-118) 


COS  Oo  = 


uJ»  -  0)^ 


(1-119) 


so   that 


and 


y^  =  - 


-2§0Or7CU 


(1-120) 


0o  -  -Itin" 


OJrf-Ui^ 


(1-121) 

For    our    example    with    the    three    different  values    of  5 
considered,    the    responses   are: 

-for    5^1          XpciKC-l)  =  0.28  cos  (3.5 -t  +5.5  rod)  (1-122) 

-{or    ^^1.5        XparlCO  =0.22   cos  (S.5l  +  5.gl   rod)  (1-123) 

:(or   5  =  0.5       Vpari(0  =0,56  cos  C3.51+5.82  rad)  (1-124) 
The   total    responses   correspondingly  are: 


,      -8.51  t 

x(t]  ^    (4+(5G.04-t)  e 

+   0.28    cos   Cs.5i  +  5.5    fad) 

.  .  -3.261  -22  28  -fc 

xd)  =    6.56  e  -2.36  <2 

+    0.22    COS  C  3.5  i  ^5.31  rod) 


(1-125) 


(1-126) 
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xil)    =  -7.76  e'  ■       cos  (757^  f  2.11  rad) 

4  0.56  cos    C5,5t  +5.32    rad)  (1-127) 

(b)    Singularity   Inputs 

The  Laplace  Transform  method  can  be  very  easily  applied, 
since  we  are  dealing  with  linear  differential  equations 
having   constant  coefficients. 

Appendix  A  gives  a  short  Laplace  Transform  Table  and  some 
Laplace  Transform  Theorems  that  help  to  transform 
differential   equations  to  algebraic. 

In   illustrating   the   above  method,   we   are  considering   unit 
impulse   and  unit   step   input    forces. 
(1)    Unit    Impulse   Response 

The    spring   mass   damper    is   driven  by  the   singularity 
input 

FCl)   =  8(0  (1-128) 

The   differential   equation  becomes 

X  f  25a;„  X  t  oUnN    =   S(U  (1-129) 

Neglecting     initial     conditions    and    Laplace 
transforming,    we   obtain: 

Xp^^^Cs)     =    Cs^.250DnS    -.CUnM    =1  (1-13  0) 

or 

XpaK  (5)   =    1/sH  25aJnS4-a;^^  (1-131) 

For   our   example  with  5=1 

Xparl^O    =    -/Cs+85i)2  (1-132) 

^par^  a)  =    t   e  (1-133) 
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and  the   total    response    is 

yc[i)  =  (A^GG.OAi)  e    '       ■^iG'^'^'  (1-134) 

with  5=1.5 

Xparl  (^)  -  Vs^.  25.55s  -72.42    =  VCs^3.26Xsv22.28)  .-   0P525_  ^  0.0525 
^  s+5.26        sv 22.28  ,.     .._, 

y^ri  Ct)  =  0.0525  e"^-^^"^  '0.0525  e""^^^  (1-136) 


and  the  total  response  is 
With  5=0.5 


Ai)  =  6.S6  .-^^''-  2.SG  e^^-'^\a0525.-^-^^^a0525  i'''^'        (1-137) 


XpuKCs)   =  Vs^^  8.51s. ^2. 42     3  yCsU.26+j7.5T)(:  5.^.2G  -j  7.57} 


ao7  e  0.0?  Q 


5+ 4.26 +j  7.57  54- 4.26 -J  7.57 


(1-138) 


M)  -  0.14  g'"'-^^^  cos  (757t  +  f57  rcidj  (1-139) 


and  the   total    response    is 

xC-l)    =-776  e'  cos  C7.S7t  4-5.HI  rad]  + 

-4  261 
+  0.14  e  cos  (7. S7+  f 57  rad}  (1-140) 

(2)    Unit    Step   Response 

The    spring    mass    damper    is   driven  by  the   singularity 

input 

F(U  -  u(l)  (1-141) 

The   differential   equation  becomes: 

y  +2S(jOr)X  4-uj/x    =   u(l}  (1-142) 

Neglecting  initial  conditions  and  Laplace  transforming,  we 
obtain: 

XpaK(s)  -  C6^-2?uJnS^a3n^)  =s  (1-143) 
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or 


For  our  example  with  5=1 

5  i  8.51 


(1-144) 


Xparl  Cs)    = 


Cs+8.5i)^  S4-8.5H        Cs+8.50^  (1-145) 

-8.su  -8.51  t 


(1-146) 


and  the   total    response    is 

x[-l)=  U  +  G6.04-<.)  e         +e  -  i  (2  (1-147) 

With  5=1.5 

<  0.17  -^.17 


Xpcirl  (0    = 


(s+5.2&)Cs+22.28)  5+5-26  s+22.28  (1-148) 


„s  -526-t      -22.28t      -3.26-t      -22.28-t         /,  ^c^^ 

x(0  =  6.56  e    -2.56  e     +  0.17  e    -  \ AT  c  (1-150) 


and  the  total  response  is 
=  6,56  e  ""-  5 

With  5  =0.5 

-30.11-  50,r 

5,  0.58  G  o.58e 

/pcrUO    =   ^^^^26tj^.57)(S4-4,26-j?.57)    "    S^4.26tjT.37  '  5  +  4,2^ -j  T^?    (1-151) 

>^paH^^)  =    ri&e'''^^^  COS  (7.57  i  ^  0.55  rod)  (1-15  2) 

and  the   total    response    is 


y(U  =  -T.Y6s'  cos  CT,57-t  ^^M  rcid)  + 

-  4  26  t 
+  1.16  Q.    '        cos  C7.5Ti  +0.55  tad^  (1-153) 
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1.6      AXES    SYSTEMS 

In  describing  the  motion  of  our  flight  vehicle,  we  will 
make   the   following  two  assumptions: 

First,  we  will  establish  an  inertial  frame  of  reference 
on  the  earth  considering  the  earth  to  be  fixed  in  space.  In 
this   inertial    frame,    Newton's   laws  are  valid. 

Second,  the  vehicle  will  be  assumed  to  behave  like  a 
rigid  body  so  that  translational  motion  can  be  described  by 
considering  it  as  a  particle  with  all  its  mass  located  at  the 
center  of  mass,  and  rotational  motion  can  be  described  by 
considering  moments  about  that  center  of  mass  and  constant 
moments   of    inertia. 

To  describe  the  translational  and  the  rotational  motion, 
one  has  to  establish  a  second  non-inertial  axes  system,  fixed 
on  the  vehicle  with  origin  at  its  center  of  mass.  In  this 
reference  frame  Newton's  laws  do  not  hold,  in  general,  and  we 
have  to  add  additional  terms  to  our  basic  fundamental 
equations. 

It  is  convenient  to  use  cartesian  coordinate  systems  for 
both    reference   frames   used. 

Let's  consider  the  translational  motion  first.  We  will 
assume  that  Oxyz  is  the  inertial  reference  frame  and  0  y'lj'i' 
the  non-inertial.  As  in  Figure  1-7  shown,  r  and  7  are 
position   vectors    pointing   to  an   object   and   originating   from 
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the    origins    of    the    coordinate    axes    correspondingly.    7^ 
denotes    the    position    vector    of    the    non-inertial    origin 
measured   from   the   inertial   origin. 


/object  .2' 


Figure  1-7 

Inertial  and  Non-Inertial  Reference  Frames 
Translational  Motion 


The  relation  between  those  three  vectors  is  clearly 

-»   -•/  -» 
r  -  r   +  To 

By  differentiating  twice,  with  respect  to  time,  we  obtain 

7   =  r '  4  r;  (1-154) 

Now  since  we  have  assumed  that  Newton's  laws  hold  on  the 

frame  Oxtjz  ,  we  may  write 

^7   ^f  (1-155) 

or   substituting   Equation  1-154 

/yn? '   =    F  -  TnK;  (1-156) 

where  yr\   is  the  mass  of  the  object  and  F  the  force  applied 

onto  it. 
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Equation  1-156  merely  says  that  we  need  the  extra  term 
onr,  to  describe  the  translational  motion  of  the  object  in  the 
non-inertial   frame. 

If  To  =0  then  the  frame  O/y'z'  is  inertial  and  Newton's  laws 
hold.      In  this  case.    Equation  1-156    takes   our   familiar   form 

nm'r'  =  F  (1-157) 

which   is   similar   to  Equation  1-155. 

Let's  now  consider  the  rotational  case.  We  will  assume 
that  ''o  =0  so  that  0  is  the  same  as  0,  as  shown  in 
Figure  1-8. 


t       J 


I' 


Figure   1-8 

Inertial   and   Non-Inertial   Reference   Frames 
Rotational   Motion 


The    position   vector    r    can    be    expressed   in   the   inertial 
frame   as 


r  =  xt  +  {j  J  4-  2  k  (1-158) 


or   in   the   non-inertial    frame 


r=:XL+-l|J-l-2< 


(1-159) 
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and 


where    (x^yrZ)     and    (  x' ,  \^\  z' )    are   the    corresponding    axes 
components  of  r . 

Similarly  the   vector    ^   can  be   expressed  as 

A  =  Ay  t  +  Ay  :j  +  Az  k  (1-160) 

or 

A    ^    ^^  1'  ^  Ay  j'  +  Ar  k'  (1-161) 

Differentiating    Equation   1-160,    1-161    with    respect   to 
time,   we  obtain 

dA        dAx  ^       dAy   A       dAz  J 

di     "    di     •■   ""  "dT  ^   ""  "dT  (1-162) 

dAj         dAx   A/       dAi)  //        dAz    A/ 

dW     ^    di      ""    "■  oK    ^     "■    di  (1-163) 

where    ( dA/^  )^    denotes   that  we  are  taking   the   time   derivative 

of   A  on  the  primed  non-inertial   system  where  the   unit  vectors 

L   ,  J   ,  K     are  constants. 

We    are    interested    now    in    the    time    derivative    of    the 

vector    A    relative    to   the    inertial   system,    but  we  would  like 

to    differentiate    Equation   1-161    instead  of   Equation  1-160, 

obtaining  thus: 

df   dA'.  A/  dAi  ./   dAz  f   .  .  c\V     A '  ^y     X  '  di' 

dT"-  lir  '  '  "dT  ^  '  IT  ^  '■^^Ti'^^^—K'^^-^  (1-164) 

or 

-^  '\^)   -nx  ^-^iZ-^T  "  "'^  ~^  (1-165) 

It  can  be  shown  that  a  more  convenient  expression  for  the 
three  last  terms  is  given  by 
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I  f  ell'     i'  dj'     ./  die'      -»    r 

where  u)  is  the  angular  velocity  vector  of  the  primed  system 
which  is  assumed  to  rotate  about  an  axis  having  the  direction 
of  CD  ,  through  the  origin,  and  with  angular  velocity  oo, 
relative   to  the   unprimed   inertial    frame. 

For    proof,    assume    for    a   moment    that   A    is    a   vector    "at 
rest"   on   the  primed   system   so  that 

(<^VdO'-0  (1-167) 

We  shall  show  that 

^^Ak    =  00  X  A  (1-168) 

Using   the   time   derivative   definition,   we  can  write: 

^iTo ^^ "■  ""  (1-169) 

Figure  1-9   shows  two   successive   positions  of    the  vector  A 

which    rotates   with    angular   velocity  ou   with    respect   to  the 

inertial   frame,    so   that  the   tip  of   the  vector    traces   a  circle 

perpendicular   to  the   direction  of  w , 

To    show    Equation  1-169,    we  must   show   that   the   directions 

and  magnitudes   of   the   vectors   on  the   left   and   right-hand   side 

are   the   same. 

Notice    that    the-   direction    of    the    vector     olA/d^    ,    by 

Equation   1-169,    is    the    same    as    the    direction  of    the  vector 

A(-i+At)-  A(-t)    ,    which    is    the    same    as    the    direction   of    the 

vector    uJxA      by   looking  at  Figure  1-9. 
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ALi*A{)  -  AW 


m 


Figure  1-9 
Time   Derivative   of   a  Rotating  Vector 


The   magnitude    of    the       fe£t-hand    side   of   Equation  1-169 


is: 


Al-o 


Ai 


or 


il-o 


=  A  oj  sincj? 


a 


(1-170) 


(1-171) 

which    is    the    same    as    the   magnitude    of    the    right-hand   side 
vector,  lujy  a1  . 

Hence,  both  directions  and  magnitudes  are  the  same,  so 
Equation  1-168  holds  and  if  A  is  not  stationary  with  respect 
to  the  non-inertial  system,  we  have  according  to 
Equation  1-165 
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—  =  —      +  ojx  A 


(1-172) 

which  says  that  the  time  derivative  of  a  vector  A  relative  to 
an  inertial  frame  equals  the  time  derivative  of  the  same 
vector  measured  at  the  non-inertial  frame  plus  the  term  u)^^ 
where  c5  is  the  angular  velocity  of  the  non-inertial  as  seen 
by  the    inertial    frame. 

Combining    translational    and   rotational   motions,    one   can 
easily  arrive   at  the   equation: 

where  Z  is  the  linear  velocity  of  the  non-inertial  with 
respect   to  the    inertial   frame. 

One     step     further,     differentiating    once    again 
Equation  1-172,   with   respect   to  time,   we   obtain 

d'A        d  (SlV      -     dA      dw      7 

=  TT    +CUX— -■(-— 7-xA       = 


/d'A\'    -    (c^iy    -.    f/dAi^    ^    rl     doJ 


|U  /  J       dt  (1-174) 


or 


'''^(#)'^^^^(#)-^(=^^^^f^^' 


dV    '[~^l  ^^    idU   --^^--'^       di   ^^  (1-175) 

which    says    that    the    second    time    derivative    of    a   vector   A 
relative    to    an    inertial    frame    equals    the    second    time 
derivative    of    the    same    vector    measured   at   the   non-inertial 
frame,    plus    some    extra   terms    involving  uJ ,  A  ,    ^^/di   ,    and 


44 


( dA /di  )       where   a>    is    again    the    angular    velocity   of    the 
non-inertial   as   seen  by  the   inertial    frame. 

Combining  translational  and  rotational  motions  one  can 
easily  arrive  at  the   equation 

where   fo   is    the    linear    acceleration  of    the  non-inertial  with 
respect   to  the   inertial   frame. 

A  special  significance  of  Equations  1-173  and  1-176  is 
given  by  substituting  ^  with  the  position  vector  r  of  an 
object.      They  are  taking  now  the   form 

dr        dr        ^    -      ^ 

or 

V    =   dvY  +  a3yr  ^t  (1-178) 

where   ^  is   the   velocity  of   the   object,    and 

dV      (d'rV      ^-.    [drV     -    /-,   ^^     du)     -     - 

or 

d    ^  Ca)'  +  2  a)  )f(^/)  +•  u3  x((JL)xf )  4- uJxf  +  a^  (1-180) 

where  a  is   the   acceleration  of    the   object. 
One   can    recognize   that 

2oi3x(V')     =   corriolis   acceleration  (1-181) 

and    exists    if   the   object   has   a  velocity    (V)^    with   respect   to 
the   non-inertial    frame.      Also, 

ujyCcoxr)    =   centripedal   acceleration  (1-182) 
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which    is    always    directed   toward  the   axis   of    rotation  and   is 
perpendicular   to  that  axis,    as   shown  in  Figure  1-10. 


axis  o( 


Figure  1-10 
Centripetal   Acceleration 

Since  we  have  assumed  that  the  unprimed  system  is 
inertial,   we  may  write 

nma    =  F  (1-183) 

and   substituting   Equation  1-180 

omCd)'  =  F  -2m[aixivy]  -'yn[cli^(ti)vr)]-'m(a)xr)  -  wOo  (1-184) 

where   <vyi  is    the  mass  of   the   object   and    F   the   force   applied   to 
it. 

Equation  1-184  merely  says  that  we  need  four  extra  terms 
to  describe  the  translational  and  the  rotational  motion  in  a 
non-inertial    frame. 

It  is  clear  that  the  forces: 
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2'm[ujx(V)'J  =  corriolis  force  (1-185) 

om[d)x((Bxr)J  =  centripedal  force  (1-186) 

om((i)yf )  =  Euler  force  (1-187) 

Tna,    =  D'    Alembert  force  (1-188) 

are   not    true    forces.      They    only    exist   due   to  the   fact   that 

our   object   is  moving  on  a  non-inertial    frame. 

According  to  our  first  assumption  that  the  earth  is  the 
inertial  frame  of  reference,  we  consider  that  it  is  rotating 
with  fixed  angular  velocity  a)  so  that  we  can  write 
Equation  1-184    as: 

^CaY  =  F  -  Z-m  ai  x(V)  -  on  wvfooxr  )  -  ^oi  (1-189) 

We  have  seen,  up  to  now,  the  effects  of  considering  a 
non-inertial  frame  of  reference  and  our  next  step  will  be  to 
orient  our  flight  vehicle  in  that  frame.  The  reason  for 
selecting  a  non-inertial  frame  will  be  greatly  appreciated 
when  one  will  see  how  easy  it  is  to  formulate  equations  of 
motion  of   a  body  attached   in   this  frame   and  moving  with   it. 

For  flight  control  applications,  it  is  appropriate  to 
select  the  origin  of  the  axis  system  on  the  center  of  mass  of 
the   vehicle. 

From  then  on  we  will  orient  the  nose  of  the  vehicle 
towards  the  positive  x-axis,  which  will  also  be  a  symmetry 
axis.  The  positive  y-axis  will  point  to  the  right  and  the 
positive  z-axis  will  point  downwards  having  the  direction  of 
gravity,     if     one    considers     straight    level    flight    as 
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Figure  1-11  shows,  so  that  the  xz  plane  will  be  a  plane  of 
symmetry. 

Vector  quantities  associated  with  this  axis  system  are 
shown  in  Figure  1-11  and  are  described  as  follows: 


VERTICAL 
AXIS 


H    posi-fiie  pi4dii»| 
X     aerodynamic  ^oicw 


L  positive  reiiivj 
p    roUin^  raieci4y 

o     vcioeriy 


M  potiiive  yawing  tuoMMiif 

r  yawing  vebcMtf 

Z  aerodynamic  ftrtc 

Mf  velocity 


Figure  1-11 
Axes   System 


Velocity    of    the    center    of    mass    with    respect    to    the 
inertial    frame  measured  on  this   axis   system 

V    :rut^Yj^+w/k  (1-190) 

Angular  velocity  about    the   center   of   mass 

_^  A. 

a?=pt+qj4.rt  (1-191) 
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Aerodynamic  force  acting  on  the  center  of  mass 

F  -  X-t  ^-rj  f  Zi  (1-192) 

Moment  about  the  center  of  mass 

G  =  Lt  +  M^  ^-Nk  (1-193) 

Directions  as  shown  in  Figure  1-11 
Also,  we  can  define:   (see  Figure  1-12) 
Angle  of  attack: 

o<  =  -tan' (^/u)  (1-194) 


projection  o(  7  -^ 
0n  xz  pJant 


Figure  1-12 
Angle  of  Attack  and  Angle  of  Sideslip 


Angle  of  sideslip 


^  =  s,n-'(V\^) 


(1-195) 

The  reference  axis  system  has  been  adopted  by  B.  Etkin, 
D.  McRuer  and  others,  and  is  usually  referred  to  as  body  axis 
system  (see  references).  The  notation  used  is  quite  common 
on  most   references. 
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Another  axes  system  adopted  by  C.  Perkins  is  shown  in 
Figure  1-13.  It  is  usually  referred  to  as  wind  or  stability 
axes  system. 


Z  Cwind) 


Figure  1-13 
Wind   Axes   System 

In  this  system  the  x-axis  points  to  the  opposite 
direction  to  that  of  the  relative  wind.  Notice  that  in  the 
undisturbed  flight  the  body  axes  coincided  with  the  wind 
axes,    but   in   the   disturbed   flight   that  may  or   may  not  happen. 

1.7      RIGID    BODY   DYNAMICS 

Rigid  body  is  a  system  of  particles  in  which  the  distance 
between  any  two  of  them  is  not  changing  regardless  of  the 
acting    forces.       For    our    flight    vehicle    the    rigid    body 
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assumption  implies  that  as  far  as  its  motion  is  concerned  in 
three  dimentions,    only  negligible  errors   are   introduced. 

We  will  consider  first  the  planar  motion  of  rigid  body 
where  we  will  find  the  kinetic  energy  and  the  angular 
momentum  expressions. 

Assume  that  the  rigid  body  in  Figure  1-14  can  move  in  the 
xy-plane  and  can  rotate  about  an  axis  perpendicular  to  that 
plane  passing  through  the   center   of   mass  c. 


Figure  1-14 
Rigid  Body  -   Planar  Motion 

We  denote  7^  ,  the  position  vector  of  the  center  of 
mass  c,  T  the  position  vector  of  any  particle  on  the  rigid 
body  and  r'  the  vector  from  the  center  of  mass  to  that 
particle. 

Usually,  the  mass  of  the  rigid  body  is  not  uniformally 
distributed    over    the    entire   volume   V.      It   is   convenient   to 
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think  for  a  moment  that  the  entire  volume  of  the  rigid  body 
containing  a  mass  on,  is  projected  on  the  xy-plane  on  an 
area  A.  Let's  denote  f^i  the  mass  per  unit  projected  area  of 
the  body.  It  is  clearly  seen  that  the  summation  over  the 
entire  area  of  all  those  masses  tjcIA  gives  the  total  mass  of 
the  body.      This  mathematically   is   expressed  as: 

From  Figure  1-14  we  have: 

r   =  fo  ^7'  (1-197) 

For  a  differential  mass  -viclA  at  which  the  f  vector  points, 
the  kinetic  energy  is  expressed  as  -  nf[dA  r  .  Integration  over 
the  entire  area  A  gives  the  total  kinetic  energy  of  the  rigid 
body. 

'^    -if^'^^'^  (1-198) 

The   quantity    f'  according   to  Equation  1-197    is   expressed 


as 


f    =  r-r  =  Cr;4.p)-(r;+r')  =  r/  +  r'  f  2  r;  r'  (1-199) 

Then   Equation  1-19  8  becomes 

Ja  -'a  -'a 

The   last   term   is   zero   since 

If'^d^    -   /r'^cIA  -r;  T-tjdA     =Hooi-r;nia  =o  (1-201) 
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We    now    express    the    velocity     r'    as    a    sum    of    two 
velocities: 

^     =   dT  "^   '  ^^^  (1-202) 

where  dr'/di  is  the  radial  velocity,  9  is  the  angular 
velocity  of  the  body,  r'B  the  tangential  velocity  of  the 
particle  ,  t  and  -J  denote  the  corresponding  unit  vectors  in 
polar  coordinates. 

Noticing   that    c^r'/dl=0.    Equation  1-200    is  written 


K=  Lyt^rJ    .IqV  r'NdA 


2  ^        '      1^'"  (1-203) 


We  define 


(1-204) 
'A 

the  moment  of  inertia  of  the  body  about  an  axis  perpendicular 
to  the  plane   of   motion  through  the   center  of   mass.      Then 

K  4-1^.'  4  19^  (1-205) 

This    equation   merely   says   that,    in   addition  to  the   term 
-  mf,      that    describes    the    translational   motion,    we   need   the 

2 

extra  term  il9^  which  describes  the  rotational  motion  of 
the  body  around  an  axis  perpendicular  to  the  plane  of  motion 
through   the   center   of   mass. 

The    angular   momentum    of    the    same    differential   mass    ^[(^1^ 
about   the   origin   is   expressed  as 

X 


}xy,6(^r  (1-206) 
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Integration   over    the    entire   A   gives    the    total    angular 
momentum  of   the   rigid  body. 


k    =   /  Crxrj'Ti.dA 
The  cross  product   term  can  be  written 


(1-207) 


r)(r  =  (f^  +  r')x(r;  +  r')  =  r;  x  f;  +r'xr'  f  r\rl   ff;  xr'        (1-208) 
Then  Equation  1-207  becomes 


(1-209) 


i[  =  roxf„  Udk  +  M  rVrjTjdA  +r„x  f  r'i|dA  ^  ?„  x  /r'r|dA 

In  this  expression  the  integral  of  the  last  term  is  zero 
as  Equation  1-201  shows.  By  a  similar  argument  the  third 
term  is  also  zero.  Noticing  that  r^=f0a)  ,  Equation  1-209 
becomes 

k   Ar^t7,]w    ^  e     r'[r'xl^)^dA  (1-210) 

If  we  use  the  center  of  mass  at  the  origin  (^  =0)  we  are 
left  with: 

K    =  e   frir^hn^^  (l_2ll) 

•A 

The  cross  product  in  the  integral  indicates  that  the 
angular  momentum  vector  is  perpendicular  to  the  plane  of 
motion.  Since  Irxill  =r,  the  magnitude  of  the  angular 
momentum   is   given  by: 


(1-212) 
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So  the  magnitude  of  the  angular  momentum  of  a  rigid  body 
rotating  in  a  plane,  about  an  axis,  through  the  center  of 
mass,  is  given  by  Equation  1-212.  The  angular  momentum 
vector  points  in  the  direction  of  the  axis  of  rotation  that 
can  be   easily   found   by  the   right   hand   rule. 

In  summary,  we  have  seen  that  the  motion  of  a  rigid  body 
in  planar  motion  can  be  described  by  a  translational  motion 
of  the  center  of  mass,  plus  a  rotational  motion  about  an  axis 
through  the   center   of   mass. 

We  will  now  look  at  the  motion  of  a  rigid  body  in  a  three 
dimensional  space.  It  generally  has  six  degrees  of  freedom 
if  not  constrained.  Hence,  we  need  six  coordinates  to 
specify  the  motion.  The  three  of  them  can  describe  the 
translational  motion  of  the  center  of  mass  of  the  body  and 
the   remaining  three,    the   rotational   motion. 

We  consider  two  coordinate  systems.  An  unprimed  system, 
as  shown  in  Figure  1-15,  whose  origin  is  0,  for  simplicity  is 
fixed  on  the  rigid  body  but  its  axes  do  not  rotate  with 
respect  to  an  inertial  frame,  and  a  primed  system  whose 
origin   and  axes   are  fixed   on  the    rigid  body. 

We  also  assume  that  the  rigid  body  has  a  mass 
distribution  per  unit  volume  (density)  which,  in  general,  is 
a  function  of  the  position  vector  f  to  a  particle  on  the 
rigid  body  having  a  differential  volume  dv'  .  The  total  mass 
is   given  by: 
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Figure  1-15 
Rigid  Body  -   Space   Rotational   Motion 


om 


=  [    pCr}clv 


(1-213) 

where   p(r)     is    the    density  as   a   function  of    r   and   <^v    a   small 
differential   volume  where  r    points   to. 

For  the  differential  volume  d/  the  angular  momentum  for 
rotation  about    the   origin   of    the   unprimed   system   is   given  by 

rx  [  pcf)dv  -f]  (1-214) 

Integration  over  the  entire  volume  gives  the  total  angular 
momentum:         ^ 

where  f,  ''  are  measured  on  the  unprimed  system. 

Let  o)  the  angular  velocity  of  the  rigid  body  about  the 
origin  of  the  unprimed  system  and  f  the  position  vector  from 
the  primed  system. 
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It  is  more  convenient  to  express  7  ,  which  can  be  seen  as 
the  velocity  measured   on  an   inertial   frame,    as 

F  =  7'  +coxr'  (1-216) 

Notice   that 

'     =    X   ^'"'''^^  (1-217) 

where    t'    and    ^     are  unit   vectors   in  polar   coordinates   of    the 
primed   system. 

It  is  clear  that  d^' /cli  =0  according  to  rigid  body 
assumption,  and  ^=0  because  there  is  no  relative  rotation  of 
the  differential  volume  dv  that  can  be  seen  by  the  primed 
system.  In  other  words,  ^'  is  a  vector  at  rest  in  the  primed 
system  so  that  we   can  write 

r   =  wxr  (1-218) 

Substituting  in  Equation  1-215) ,  we  obtain 

k  =r  r   f  X  (ujxr)  p(f)  dv  (1-219) 

Jv 

Using  now  the   tripple   product    identity 

Ax(Bxc)   =  CAC)i  -  CA-B)c'  (1-220) 

we  can  write 

fx  C^xf)    =    r'(?  -  (^.f)?  (1-221) 

so   that  Equation  1-215   becomes 

k   =j    [  r^^  -  foSrjrJ  pCrjd/  (1-222) 

In  component   form  we   can  write 

h^=|    [r^w^  -C5-f)  rj  pcOdv  (1-223) 

where    1=1,     2,     or     3    denotes    the    x,    y,    or    z    component 
respectively.      The   dot  product    is   expressed  as: 
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cJ3.r     =  u;,f,  i.fj},r,  *(x)^r^   =£^  rjoOj 


(1-224) 


That  permits  us  to  write 

3 


J'i 


where  mj  =1  if  i=j  and  oij  =0  if  i^j  because 
3 


(1-225) 


We  can   now  define 


(1-226) 


^4  =l[^'^'j-nr,]pcr;d/ 


'y  -  -  •—  (1-227) 

as    the     Lj  th    component    of    the    inertia   tensor    I    about    the 
origin  of   the   unprimed   system. 

So    finally,    we    express    the    c  th  component  of   the   angular 
momentum  as : 


3 

r — 

In  the  expanded  matrix  form  this  gives: 


(1-228) 


k 

■*xx 

ixy    Ixr 

"Jx 

s 

= 

Ix, 

lyy     ^^ 

CUy 

k-z 

Ixz 

Jyz      Iz? 

C02 

(1-229) 


This  result  is  compared  with  the  result  of  the  planar 
motion  case  that  is  expressed  by  Equation  1-212,  where  the 
inertia   scalar   term   is   expressed  now  as   a  matrix. 

The  components  of  the  inertia  tensor  are  given  by 
Equation  1-227    in   cartesian   coordinates   as: 
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7//    Cy'+2')  p^r)  dxcJydz 


Similarly, 


(1-23  0) 


lyy  =/  /  (x^z^)  p(r)  dKcJydz 


(1-231) 


Also, 


Izz  ^ffi   xVl/]   p(f)  dx^y^z 


(1-232) 


Ixij  =  ^(j> 


xy    p(r)  cJxdijciz 


(1-233) 


Txz  -  Izx  "-///^^    pC*^)  dxciydz  (1-234) 

Example:  Calculate  the  inertia  tensor  of  a  cube  of 
side  a,  mass  m  and  uniform  density  p  about  one  of  its  corners 
(Figure  1-16) . 


Figure  1-16 
Inertia  Tensor-Calculation  of   a  Cube 
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U^'   ?l[[ii'*^')'>''ii'^'    ^^j  j  f    ff-^'Jdzdydx^ 


O    ■'o      » 


^7^f[^'i''^^T)^^'^^^.fiT-i''^'^^' 


'r<^    [  or  „       (^       \         m 
-  —  \   — a  +  —  a      r    — 
a*  I    5  5      i         a' 


-yyi     2a^ 


a'     3 


2  2 

—  oD  a 
5 


(1-235) 


Similarly, 


T  T  2  ^2 


(1-236) 


a  ^a  ,a 


Ix,^  -- -pjjfjT^y  d^^^y^^  -'-^{/j[   ^!iJ^^y^^  = 


^M.r  '^°^^^^  -'^i 


5     12 


Similarly, 


cm    a^    a}  <  2 

= —    =• on  a 

a3    2      2  4 


1x2    =  'T  '^^' 


(1-237) 


(1-238) 


so   that   the    inertia   tensor    is  written: 


T  =  wa^ 


2/5        -'/4  -^4 

-'/4         2/5        -V4 
-V4         -'/4         2/3 


(1-239) 
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The  eigenvalues  and  the  eigenvectors  of  an  inertia  tensor 
have  a  very  special  significance.  The  eigenvectors  determine 
the  directions  of  the  so-called  principal  axes  of  a  rigid 
body  that  pass  through  the  point  where  the  inertia  tensor  is 
calculated  (usually  the  center  of  mass).  The  corresponding 
eigenvalues   are  called   the  principal   moments  of    inertia. 

It  can  be  shown  that  if  a  rigid  body  rotates  about  one  of 
its  principal  axes,  its  angular  momentum  about  the  point 
which  the  inertia  tensor  is  calculated,  is  parallel  to  its 
angular  velocity. 

For  proof,  let  ^  a  principal  axis,  and  B  the 
corresponding  principal  moment  of  inertia.  The  eigenvalue 
problem   equation  which   has   to  be   satisfied   is: 

I^    =  B^  (1-240) 

or 

Itj  /^j     -"B/S,  (1-241) 

Having  rotation  about  a  principal  axis  simply  means  that 
u>  has  the  direction  of  the  principal  axis  which  is  expressed 
as: 

^    ^  c^  (1-242) 

where   C    is    a    constant.       In    fact,    if    /S   is   unit   vector,    C   is 
the  magnitude   of   the   angular  velocity  u5. 

Now  the   ith  component  of    the   angular  momentum  becomes 

k,   --I.jojj    .  I,j  Cc|?j)  =  citjfj  =cS^,  =  BCc^J  (1-243) 

and   finally, 

k[   =  BuJ,  (1-244) 

61 


In  vector   form 

h    =  B  uJ  (1-245) 

which  proves   that    k  is  parallel   to  OJ. 
It  can  also  be   shown  that 

1.  Any  plane  of  symmetry  of  a  rigid  body  is 
perpendicular   to  a  principal   axis. 

2.  Any  axis  of  symmetry  of  a  rigid  body  is  a 
principal   axis. 

The  density  p  may  not  be  constant.  For  example,  if  the 
xy  plane    is  a  plane   of   symmetry,    that   simply  means  that 

p(-^,y,z)    =  pCx,y,z)  (1-246) 

and    if    the    z-axis    is    an    axis   of   symmetry,    that   simply  means 
that 

p(-x,-j^,2)    =    p  Cx,y,z)  (1-247) 

3.  Since  the  inertia  tensor  of  a  rigid  body  is 
symmetric,  (i.e.  Jxy  =  Jyx  etc.),  the  principal  axes  are  always 
orthogonal  to  each  other  (Important  result  from  Linear 
Algebra) . 

4.  If  two  or  three  eigenvalues  of  a  symmetric  tensor 
are  the  same,  there  exist  two  or  three  orthogonal 
eigenvectors   associated  with   that  eigenvalue. 

Example:  The  inertia  tensor  of  a  rigid  body  has  been 
calculated  with  the  reference  point  at  the  center  of  mass  and 
has   been  found   to  be 
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\5 

0 

-iO 

0 

25 

0 

-<o 

0 

30 

1  = 


Find  the  principal  axes  and  the  corresponding  principal 
moments   of    inertia. 

Solution: 

Let  J  denote  the  eigenvector  (principal  axis)  and  A  the 
corresponding  eigenvalue  (principal  moment  of  inertia) . 

The  following  eigenvalue  equation  has  to  be  satisfied 


from  which 


(A-I)  J  =  0  ^ 


-^J 

.AJ 

A-15 

0 

W 

0 

A-25 

0 

\0 

0 

A- 50 

=  0 


or 


(A'\5)(A-2^)(A-lO)  +  lO[.{0[A-25)]   =  0 

^  (A-25)[C/\-iS)CA-iO)  -m]    =0 
=>  (A-2S)CA-lS)(A-iO)    =0 
Hence,  the  principal  moments  of  inertia  are: 
A,  ^  \0  Aj  -  25 

To  find  the  eigenvectors: 
For  A^  =10  Equation  1-249  becomes 


J.  -0 


^5=  55 


-5 

0 

W 

0 

-6 

0 

>iO 

0 

-20 

(1-248) 


(1-249) 


(1-250) 
(1-251) 
(1-252) 

(1-253) 


(1-254) 
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which  yields 

The  eigenvector  that  satisfied  Equation  1-25  4  is 


(1-255) 


.1.  = 


2 
0 


(1-256) 


It  is  more  convenient  to  express  the  eigenvector  as  a 
unit  vector 


.% 


Similarly,   we   find   that 
0 

i 


0 


(1-257) 


i 
0 


assosiaied     wi-lh        A^  =  25 


(1-258) 


and 


Vis 
0 


assosialeol     vui4(t 


A, -55 


(1-259) 


It  is  important  to  note  that  the  above  eigenvectors  are 
orthogonal   since   they  are  coming   from   a   symmetric  tensor. 

Since  now  the  principal  axes  of  a  rigid  body  are 
orthogonal   to  one  another,    or   can   be   chosen  orthogonal,    it    is 
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very  convenient  to  have  them  coincide  with  the  primed  system 
axes,  the  origin  of  which  can  be  moved  to  the  origin  0  of  the 
the  unprimed  system.  The  importance  of  this  arrangement  will 
be  shown  a  little  later. 

The  kinetic  energy  of  the  differential  volume  dv  is 
2  (r  f)p(r)dv.   Integration  over  the  entire  volume  gives  the 
total  kinetic  energy: 

^       2  }y^  ^  (1-260) 

Substituting 

f  =   (]3x?  (1-261) 

we  obtain 

K  r  if  (d)xr)(ujxr)pCfjdv  (1-262) 

Using   Equation  1-224,   1-226    and   the   vector    identity 

(AxB)(c>cP)  =  CXcJCS-D)  -  CaJjCBC)  (1-263) 

we  can  write 


=   lb    [^'^M-^^^jl^^Jf'^^  (1-264) 


so  that 

3        f      3 


K-i  1.^1.  \Mk-'^'^^^^r-^ 


(1-265) 


and   finally. 
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K  =  i  1(  Xhjoajuj,    -^}_K^v    4(hii;)  (1-266) 

u-l       J--1  t-'J 

Notice  that  if  the  body  is  rotating  about  a  principal 
axis  (Kllu))  its  kinetic  energy  is  maximum  K=  j  hu)  .  That 
happens  when  the   rigid  body   rotates   about   a  principal   axis. 

This  compares  with  the  result  of  the  planar  motion  case 
that  is  expressed  by  Equation  1-205,  if  we  neglect  the 
translational   motion. 

We  are  now  in  a  position  to  formulate  the  rigid  body 
equations  of   motion. 

The  translational  motion  equation  is  given  by  the 
equation  of   motion  of    the   center   of   mass  of    the   rigid  body 

oyir;     -   F  (1-267) 

or 

y  c^i  I  (1-268) 

The  rotational  motion  about  the  point  (0)  which  is 

conveniently  chosen  as  the  center  of  mass  (c) ,  as  Figure  1-17 
shows,  is  given  by 

dT  "  ^  (1-269) 

where  G   is  the  net  external  moment  about  the  center  of  mass 
(c)  and  ^   is  given  by  Equation  1-228. 
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Figure  1-17 
Rigid  Body  -   Space   Motion 


The    vector    oK/di     is  measured   in   the   unprimed   system,   but 
it  can  be  expressed  as: 

^   idTJ    *"''^  (1-270) 

where  ( dit/d-t  )  is  the  same  quantity  measured  in  the  primed 
system,  and  ou  is  the  angular  velocity  about  the  center  of 
mass. 

So   Equation   1-269    is  written   in   component   form   as: 


(^)    W-'^X    -G. 


where    t     refers    to    a    primed    axis    now. 
Equation  1-228 


(1-271) 
Or    better,    using 
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_3_ 

j7f^  ^'  (1-272) 

where    I  ij     is    a    component    of    the    inertia   tensor    calculated 
along  the  primed  axes. 

Notice  that  the  primed  axes  unit  vectors  are  vectors  at 
rest,  in  the  primed  system,  and  obviously  I  ij  is  not  changing 
with   time.      Hence, 

I    h    -Z     '^'^'''^^^     -<^i  (1-273) 

This  result  is  interpreted  as  follows:  Suppose  the  rigid 
body  rotates  about  a  principal  axis  with  constant  angular 
velocity  u5 ,  Rotation  about  a  principal  axis  implies  ai 
parallel  to  ^  ,  so  that  i5^l^=0.  This  simply  means  that  the 
moment   required  for   this   rotation   is   zero. 

This  compares  with  the  planar  motion  case  result  which 
comes   by  taking  the   time   derivative   of   Equation  1-212. 

k   r  lob    =G  (1-274) 

Equation  1-27  3  is  the  equation  of  rotational  motion  in 
the  general  case  for  rotation  about  a  random  axis,  where 
external  moment  is  required  even  for  constant  angular 
velocity. 

It  is  now  clear  why  we  formulate  the  equations  of  motion 
in  the  primed,  non-inertial  coordinate  system.  We  have 
avoided  the  time  differentiation  of  the  inertia  tensor  in 
Equation  1-272   by  paying  a   little   for   the   extra  term   a)xk. 
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In  summary,  we  found  the  equations  of  motion  of  the  rigid 
body  described  in  the  convenient  form  which  will  be  used  to 
derive   the   flight  vehicle's   equations  of   motion  in   Chapter  2. 

They  are: 


anl  —  t  uJx  V      =r  F 


and 


where 


di  /  (1-275) 


J--1 


3 


1^1     =   ^    lij  u^j  (1-277) 

J--1 
and    Itj      is    calculated    about   the   center   of   mass  with   respect 

to  the  primed  axes. 

Finally,    if    the   primed   axes    are  principal    axes    of   the 

rigid  body   it   is  not  difficult   to  see  that  the   inertia   tensor 

has   only  diagonal   terms   representing   the  principal   moments  of 

inertia     (eigenvalues).       This     result    comes    from    the 

diagonalization  process  of   a   square  matrix.      So,    in   this   case 

Equation  1-276    is   simplified  as: 


where 


h     ^'    ^    ('^xkj,     =G.  (1.278) 


k    -III  "^i  (1-279) 
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CHAPTER  2 
VEHICLE'S  EQUATIONS  OF  MOTION 

2.1      INTRODUCTION 

With  the  methods  of  obtaining  rigid  body  equations  of 
motion,  established  in  the  preceding  chapter,  we  now  proceed 
in  this  chapter  to  formulate  the  equations  of  motion  in  a 
most  generally  useful   form. 

The  equations  are  basically  simple  in  the  sense  that  they 
are  derived  by  making  use  of  the  basic  Newtonian  laws  for 
forces  and  moments.  Great  advantages  are  obtained,  as  we 
have  seen,  by  using  a  non-inertial  frame  of  reference, 
because  with  an  axis  system  fixed  in  the  body,  the  inertial 
properties  are  constant  and  that  simplifies  the  moment 
equations. 

Some  penalties  are  obtained  by  introducing  non-inertial 
frames,  the  most  important  of  which  is  the  addition  of  extra 
terms  in  the  basic  equations  of  forces  and  moments,  in  trying 
to  make   the  basic   laws  valid   in   this   environment. 

In  any  case,  although  the  equations  are  simple,  they  are 
highly  non-linear  and  not  generally  solvable  mathematically. 
The  formulation  is  greatly  simplified  by  the  use  of 
linearization  techniques,    and   small  pertubation  theory. 
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This  method  assumes  an  equilibrium  position  which  is,  in 
fact,  a  solution  to  these  equations.  That  allows  for  the 
changes  in  forces  and  moments  to  be  expressed  as  linear 
functions,  in  terms  of  small  pertubations  in  linear,  angular 
velocities,  and  other  parameters  involved. 

The  equations  of  motion  will  be  derived  for  the  body  axis 
system  established  in  Figure  1-11  and  will  be  left  in  their 
dimentional  form.  Non-dimensional  form  will  be  mentioned  as 
it  applies  to  body  axis  system  (Etkin's)  and  wind  axis  system 
(Perkin's) . 

Finally,  the  equations  of  motion  will  be  split  into  two 
separate  groups,  the  longitudinal  and  the  lateral  groups  of 
equations. 

2.2   FORCE  EQUATIONS 

Starting  from  the  vector  force  equation 

di  '   '  '  (2-1) 

where 

ZF   =  IX';  +  ZY^^  +IZic  (2-2) 

y  ^    ul   fyj4-wk  (2-3) 

UJ    -  -pt  +  q  J  +  r  k  (2-4) 

we  can   derive   the   force   equations   in   component   form. 
The   cross  product   term   is   expanded  as: 
(3x  V   =  (qw-yr)  i  +(-pw+ur)j+  (^v-qu  j  k  (2-5) 
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IF   =^(  J>(i3xi7) 


and  the  force  equations  of  motion  are  simply 

f  X"  =  w  (  u  f  q  w-vr  ) 
2Y  =  w  (  V  -pw  i-vr) 


(2-6) 
(2-7) 
(2-8) 


These  equations  describe  the  translational  motion  of  the 
center  of  mass. 

Consider  the  motion  of  the  center  of  mass  in  the  xz  plane 
as  Figure  2-1  shows,  in  two  successive  different  locations  on 
the  flight  path  during  time  hi. 


^■x 


Figure  2-1 
Motion  of   the  Center  of   Mass   in   the  xz-Plane 


The   acceleration  along  the  x-axis   is   expressed  as: 


a,  = 


51  (2-9) 

For  small  pitch  angle  0  ,   cosSS  =1  and  sin69  =56  ,  so  that 


a. 


hi 


(2-10) 
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letting    now     hi-^O   and    neglecting   the   higher   order   term  b\uod  , 
we  obtain 

a^  :r    U-h  wcj^  (2-11) 

and   for   the   corresponding   force   equation 

Zy   =    w  (  U^■W£j^)  (2-12) 

which   is   the   same   as   Equation  2-6  by  neglecting   the   term  vr . 

Similarly,    the   acceleration  along   the   z-axis   is   expressed 
as: 

(w+8w)  cos ^6  -  iu+hv)  smhB  -  \aj 

U.  7        =         


which  is  simplified  as: 

a,   =      -  ,2-14) 

or 

a^  -  w -au  (2-15) 

and   the   corresponding   force   equation  becomes: 

11    -   on  C  W  -qv)  (2-16) 

which   is   the   same   as   Equation  2-8  by  neglecting   the   term  pv . 

If  we  now  consider  the  motion  of  the  center  of  mass  in 
the  xy  plane  as  Figure  2-2  shows,  in  two  successive  different 
locations  on  the  flight  path  during  time  hi,  we  can  derive 
the   acceleration  along  the  y-axis   as: 

(v+§v)cos3ip  ^  (v^hv)  sm^w  -V 
^il  --    -g (2-17) 
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Figure  2-2 
Motion  of   the  Center   of   Mass   in   the   xy-Plane 


For   small  yaw  angle  ]^ ,   cosSi(7=l   and   sin^J<>  =  5TVf    so   that 


a. 


hi 


(2-18) 


letting    now    ^-♦O     and    neglecting   the   higher   order   term  huhil}  r 

we  obtain 

ay  r  Y  f ur  (2-19) 

and   for    the   corresponding   force   equation 

lY    =  ^  (v^ur)  (2-20) 

which  is  the  same  as  Equation  2-7  by  neglecting  the  term  pw. 

In  summary,  the  force  equations  are  simplified  neglecting 

higher  order  terms  as: 

IX  =  w  Cu  +  qw)  (2-21) 

2Y  =  ^  (Wur)  (2-22) 

11    =    ^   (w-qu)  (2-23) 
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2.3      VEHICLE    SYMMETRIES 

Up  to  this  point  we  have  made  two  assumptions.  The  first 
assumption  was  that  earth  is  considered  to  be  fixed  in  space 
so  that  we  can  establish  an  inertial  reference  frame  on  it. 
The  second  assumption  was  that  the  vehicle  is  considered  to 
be   a   rigid  body. 

Facts  are  simplified  further  by  considering  the 
symmetries  of  the  vehicle.  For  an  airplane  or  a  missiler  the 
xz-plane  is  a  plane  of  symmetry,  if  one  assumes  that  the 
following  holds   for    the   density  distribution. 

p(x,y,2)    =r  p  (x,-y,£)  (2-24) 

so  that 


p  (^>y»2j  yy  dxcjyc^z  ^- ///pf^-y,2)xy  -^x^y^^ 


(2-25) 

Notice   that  Equation  2-24    says  that    p   is   an   even   function 
with    respect    to   y   but    the   integrant   function  is   odd  because 
it    is    a   product    of    an    even    and    odd    function.       Since    the  y 
integration  has   symmetric   limits,    the   integral    is   zero. 
For    the   same   reason 


Jijz  -  Izu  - - 


pCXjij^z)  LjZ    clxdijc^z     ^  0 


Jyz^iZLj       j^i  I      .J,    •  J  ^  (2-26) 

The    inertia   tensor,    therefore,    for   an   airplane  becomes: 

(2-27) 


Jxy 

0 

-Ixz 

0 

Iljij 

0 

-1x2 

0 

Jzz 

75 


Ixx 

0 

-Ivz 

0 

0 

0 

lyy 

0 

J 

=   Jijy 

1 

-u 

0 

Izz 

0 

0 

Notice  also  that  the  y  axis  is  an  eigenvector  of  the 
inertia  tensor  (Equation  2-27)  because  it  satisfies  the 
eigenvalue  matrix  equation: 


(2-28) 


Hence,  the  y-axis  is  a  principal  axis,  with  Ii/y  being  the 
corresponding  eigenvalue  or  principal  moment  of  inertia. 
This  was  expected  since  y-axis  is  perpendicular  to  the  plane 
of   symmetry  xz. 

Unfortunately,  this  is  the  only  symmetry  that  exists 
generally   in   any  conventional   aircraft. 

In  addition  to  this  symmetry  another  symmetry  can  be 
found  in  gun  projectiles  and  in  most  missiles  and 
space-crafts.  They  have  a  second  plane  of  symmetry,  the  xy 
plane,    so   that  one  can   assume 

p(K,y,z)    =    p  (x,y,-2)  (2-29) 

and  hence, 

Uz  =  hy-0  (2-30) 

Then,    the    inertia    tensor    for    a    two-symmetry    vehicle 

becomes 

0  0 

lyy         O  (2-31) 


Ixv 

0 
0 


iZ2 
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Notice  that  in  addition  to  the  y-axis  being  a  principal 
axiSr  the  z-axis,  as  well  as  the  x-axis,  are  all  principal 
axes  with  corresponding  principal  moments  of  inertia  Izz,  and 
Ixx ,    respectively. 

2.4      MOMENT    EQUATIONS 

Starting   from   the   vector   moment  equation: 

or   in  component   form 


where 


dt  ^----  (2-33) 


K    -  I.J  UJ,  (2-34) 

k     ^    hx  t   +hy  5  ^  hz^t  (2-36) 


d5  -A       .    . 


A. 


,       ,       +     f"     ^ 

di  '  ^'  (2-37) 


Equation  2-34    is   expanded  as: 

ky    -  pixx   -  r  Ixz  (2-38) 

Ky    =  Cj^Iyy  (2-39) 

hz    =  --pixz  +rlz2  (2-40) 

also,    the   cross  product    terms   is   Equation  2-33  is   expanded   as 

(ujxM  -(^H^-rhy)!  ^(rh,-^hj^  +(p^-(j^k)ii  (2-41) 
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or 


"•■  ( a  p  lyy  -  c?  p  ^xx    f  rq  Ixz  )  k 


or 


(ujxk  )     =    [-^plxz     t  ir  (Tzz-Iyy)]  t 
■^  I  P^  (Jyy-Ixx)    4  ra    Ixz]  k 


(2-42) 


(2-43) 


Thus,    the  moment  equations  are: 

L   =   Ixx  p   '  J^z  f   -^pl^2    ^  (ir  (I^2-IyL|)  (2-44) 

M  =  lyy  ^    +  rp(  Ixx-Izz)    +IxzCp2-r2)  (2-45) 

•^   =    Izz  r  -  1x2  p    +  p^    (lyi/  -Ixx  )  4-  rq  Ixz  (2-46) 

2.5       EXPANSION   OF   GRAVITY   FORCES,    THE    LINEAR   AND   ANGULAR 
VELOCITY 
The    total    force    acting    on  a   flight   vehicle   can   be   split 
into    a    summation    of    aerodynamic    forces    (lift,    drag    and 
thrust)    and   a  body  force  which   is   the  gravity   force. 
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Figure  2-3  shows  an  accidental  orientation  of  the  vehicle 
with  respect  to  the  vertical.  The  gravity  force  \V  can  be 
expanded   in   components   along  the  body  axes   as   follows: 

Wx    =  -rfn^  sme  (2-47) 

Wy    =     fw\q  cos  Q  s\n  (p  (2-48) 

Wz    -      'Wg   cos9  cos 4)  (2-49) 


vtrHcot   pkme  C0tt4ainirig  x-a«(s 
(MmI  wigM  vccior 


Figure  2-3 
Expansion  of  the  Gravity  Force 

For  the  expansion  of  the  linear  velocity  vector  along  the 
body  axes  consider  in  Figure  2-4  the  inertial  reference  frame 
being  ^>^4yi2i  .  Transformation  to  the  body  axes  coordinate 
system  is  obtrained  by  rotating  the  reference  frame  in  an 
order-dependent  manner   as   follows: 

(a)  Rotate  around  reference  axis  Oz^  by  the  yaw  angle  i^ , 
which   defines    the    intermediate    reference    system     OvjyjZ^  = 
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(b)  Rotate    around    intermediate    axis     ^y^     by   the   pitch 
angle  0,   which  defines   a   second   intermediate   system      C'xjyjZi    = 

(c)  Rotate    around    final    axis    Ox    by  the   roll  angle  4> ,    to 
define   the   desired  body  axes   system    O^y^    =  Oxjyz    . 


x«x 


^'*. 


^i 


'---  /^ 


•H2  =  M, 


2^.2, 


0    -*   -pi+ch 
(t>    -»    roll 
1^     -»      \^QSU 


Figure   2-4 
Expansion  of    the   Linear  Velocity 

The  components  of  the  velocity  vector  in  the  body  axes 
system  (U,V,W)  are  defined  in  terms  of  the  vector  components 
in  the  fixed  inertial  frame  (U^ ,V,  ,W, )  by  means  of  the  angles 
®,<t?,V   and    the    intermediate  vector   components    (UifVzyWj)    and 

(U3,Vj,W3). 
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The  following  relations  hold: 

for  the  first  rotation 

V^  =  Ui  cosTp  +  Vi  sm-ip  (2-50) 

V2  =  -V^  sin-i};  +  V^  cos^  (2-51) 

W^  =   \Mi  (2-52) 

for   the   second   rotation 

Uj-  1)2  cos  e  -  W^  sir7€)  (2-53) 

^3=  ^2  (2-54) 

W3=  i;^    sme  +  W2  COS0  (2-55) 

and   for    the   third  rotation 

U  ^  U3  (2-56) 

V   -  V3    C05  4>  +■  ^5    sm43  (2-57) 

W  =  -V3    sm4;   +  Wj   cos 4?  ^2-58) 

so    that    the    velocity   components    along    the   body  axes    are 
expressed  as 

0  ^  Vi  cosipcos©  +  Vj  smii^cosQ  -  IWj  s\n6  (2-59) 

V  ^  Vj  (cosTpsmO  sin(i?  -  sir7ipGos4?)  +V^  fsmip si n0 si r? <i)  ^- cos ^p cos 4^) 

+  W,  cos9  simp  (2-60) 

W  ^Uj  (cos  IV  sm0  cos  4?  -t-smzp  siy)4))  +7^  fsm^  sm€)  cosct)  -  cos  ip  s)r)  4> ") 

-^  W,  cos0cos4>  (2-61) 

Similar  approach  can  lead  us  from  the  body  axes  system  to 
the  inertial  for  which 

^  '  I'cosQcos^  +  \/Csin(>siyi9cc7Sip  -  cc>s4>  smip ) 

+  v\/(cc7s4)  sin©  cosip  +  smcp  sin  ^  )  (2-62) 
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+  ^  (cc?S(l>  sm9  sinq;  -  5^4?  coSTp)        (2-63) 
Wj   =  -UsinG  +-  V  sm4>Gos0   +  W  cos 4) cos 0  (2-64) 

Integration  of  Equation  2-62,  2-63  and  2-64  give  the 
expressions  for  the  translational  motion  but  in  other  than 
special    simple   cases   they  are  not   integrable. 

Following  a  similar  procedure  we  can  expand  the  angular 
velocities   Prqrr   into  components   along  the  body  axes. 

With  the  vehicle  in  horizontal  flight,  one  can  see  the 
yawing-velocity  vector  having  the  same  direction  as  gravity, 
the  rolling  velocity  vector  having  the  direction  of  the 
longitudinal  axis  and  the  pitching  velocity  vector, 
perpendicular  to  both  as  shown  in  Figure  2-5.  Those 
orientations  are  obtained  according  to   right-hand   rule. 


'% 


4lort20(tWI 


vtr^i    pUwe  covrlainina 
and    uieioM  vcdor 


%-axis 


vwiu  an    yj   ptoine 


Figure  2-5 
Expansion  of  the  Angular  Velocity 
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Then   the   angular   velocity  vector  UJ  can  be   expanded   in   the 
body  axes   as   follows: 

'p  =  (p  -  li^siKlG  (2-65) 

q  -    9  cos(t) +- li^cosG  sm45  (2-66) 

^  =  -0  sin43  +  V  cosO  cosct?  (2-67) 

One    can    now    solve    the    Equations  2-65,   2-66,    2-67    for   0, 

Multiplying    Equation   2-66    by   cos  (j>  ,  Equation    2-67    by 
sin  <t)  ,    and   subtracting,   we   obtain: 

0  =  q  cos4)  -r£iY](^  (2-68) 

Multiplying  Equation  2-66  by  sin  4?,  Equation  2-67  by 
cos  (p  ,    and  adding,  we  obtain 


^  =  1 


sin  4?  cos4j 


f  r 


cos©  cose  (2-69) 

Then  substituting  Equation  2-69  into  Equation  2-65,  we 
obtain 

4)    =p  +  q   siTiij? -tan©  ¥  rcos(tj4ciY79  (2-70) 

Integration  of  Equations  2-68,  2-69  and  2-70  give  the 
expressions  for  the  rotational  motion  but  in  other  than 
special   simple   cases   they   are   not   integrable. 

We    can    now   modify    the    vehicle's    equations  of   motion  by 
considering   the  gravity  as   follows: 
Force   equations: 

X  =  /vn  r  {)  -hcnto  -vr  +  q  sm9  )  (2-71) 

Y-    'm(v+ur-pw-acosesmct?)  (2-72) 

Z  =    /w  (  \AJ  +'pv  -qu   -a  cos£)  cos<p)  (2-73) 
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where    now    the  X,Y,Z   forces    represent  only  aerodynamic   forces 
(i.e.,    contributions  of    lift,    drag  and   thrust). 
Moment  equations: 

L  =  Tyy  p   -  Ixz  r    -qplxz    ^- qi^  (izz  "lyy )  (2-74) 

M  -lyy  cj^  +rp  (Ixx  -Izz)    +  1x2  Cp'-^^)  (2-75) 

N  ^   Izz  r  -  Ivz  p  +  fc^  Clyy-lxx)  f  rdj^Ixz  (2-76) 
with   auxiliary   relations: 

'p  =^  4?  -  ij^sinG  (2-77) 

<3  -  ©  cos 4?  -f-  i];  COS0  sin(t)  (2-78) 

r  -  -  9  sir)4>  +  7J^  cos9  COS43  (2-79) 

Notice    that    the    gravity    force    contributes    zero  moments 
because    it   acts   at  the   center   of   mass. 

2.6      LINEARIZATION   OF   FLIGHT  VEHICLE'S    EQUATIONS 

The  three  force  and  the  threee  moment  equations  derived 
in  Section  2.5  must  be  linearized  for  use  in  the  stability 
and  control  analysis.  It  is  assumed  that  the  motion  of  the 
vehicle  consists  of  small  deviations  about  a  reference  steady 
flight   path. 

A  steady  flight  path  is  described  by  the  general 
equations  of  motion,  if  we  consider  all  time  derivatives 
zero.  We  denote  by  capital  letters  and  the  subscript  (o)  , 
the   steady   flight   conditions,    obtaining  thus: 

Xo    =    'YnCQoW,  -  Vo^o    4^  sin  60)  (2-80) 

Yo    -    wCUX   -  PoWo    -^cosQo  sm4>  )  (2-81) 

Zo  =    w  C?oVo   'QoVo    -^cosGo  cos<^^)  (2-82) 
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U--QoPo  I>2    +  <3o^o   Clz2-Iyy)  (2-83) 

Mo=     ^Z    (Ixx-Izz)    +Ixz    (Po'-^/j  (2-84) 

No  =    ToQc  (lyy-Ixy)    +  RoQo  Ixz  (2-85) 

By  subtracting  from  the  original  equations  the  above 
steady  flight  equations,  we  obtain  the  perturbed  equations. 
For  example,  subtracting  Equation  2-80  from  Equation  2-71  we 
get: 

X-Xo   ^'w  [  OtcjW-QoWo -vr  fVoRo  +  g  (  smG  -  smGo)]  (2-86) 

Notice   that 

Similarly 


(2-87) 


Vr  =  VoRo  +VoSr  4-^oSv  (2-88) 

Also 

sin0^  sinC0o+60)  =  smQoCoshQ  -(-cosOo  swdd  (2-89) 

For   small  perturbation  quantities.    Equation  2-89   becomes: 

sir)6  =   s\Tn9o  +(cos9o)(Se  (2-90) 

Combining  all   these.    Equation  2-86   becomes 

6y  =  ^[  l)4-a5cj^fQo5w;  -VoSr  -^oSv  fC^cosa)Se]  (2-91) 

Although   we   are  deriving  perturbed   equations,    it   is  more 
convenient   to  drop   the   differential   symbol  a  . 

Similar    work    leads    us    then    to    the    following    set    of 
linearized  equations: 
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X  =  w[0  f^ocj^  4.(3o*/   -Vor  -Rov  f  (^coseo)O]  (2-92) 

V  =  w  [i/  +  Uor  +  ^oU  -Wop  -1^  w   -C^cos6?ocoi4)o)4?f(35ir?asinJ^o)o](2-93) 
Z-  w  [w  + V„p4.t^u  -OoQ-QoU  -(-(^cosQocos  4)4>  +(qsinQ, cos A»)&J^ 2-94) 

L^Ixxp  -Ixzr  -(Qop +?o^)Ixt  f(Qor  +  Ro^)(lK-Iy!/)  (2-95) 

M=Iyy(^  +  (R.p  +  ?or)(Tyx-l22)  +I>tz  C  2Pop  -  2l?o  r  )  (2-96) 

t^  =  Jzzr  -  Ixzf)  +(Po<}^+0.p)(Iyi|-Ixx)4C^o<^fQor)Jy2  (2-97) 

In   addition,    for    the    auxiliary    relations   described   by 

Equations    2-11,    2-78,    2-79,    the    corresponding    perturbed 

equations  are: 

^^  ^  -Ti;s)Y)Oo  -C^o  cos0o)0  (2-98) 

q  =  0  cos  6o  -  C^o  s)T) Go  sm  4)9  +  (4  cos  0o  cos  4  -  0»  sin  4 )  4?  -«- 

+  Tj)  cw a  sin  Cbo  (2-99) 

r  =  i^zos  Qo  cos  4>o  -  C^o  cos  ©o  sin  4  +  0o  cos  4)  4^  -  0  sir)  4 

-Cll)o  sm0o  c:os4?o)e  (2-100) 

Although   we    have    obtained    linear    equations,    they    are 
extremely   complex  and  not   suitable   for   control   analysis. 

A   further    simplification  is   obtained   if  we   consider   the 
reference   flight   path   to  be   straight,    i.e. 

IVo   -  ©o    =0  (2-101) 

with   no  angular  velocities 

X    ^  Qo    -  ^o    ^  0  (2-102) 

wings   level   flight,    i.e. 

(t?«=0  (2-103) 
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symmetric,    i.e. 

%  =0  (2-104) 

Vo   =  O  (2-105) 
With    this    assumption,    the    reference    flight    path    is 
described  by  the   equations: 

X,    =   ong  sin  6o  (2-106) 

Yo    =   0  (2-107) 

2o    =  -onQcos€)o  (2-108) 

lo    =    Mo    ^   No     -  0  (2-109) 
while   the   perturbed  equations  become 

X  =  TM  [  i)  +  WoCj^  +    ^cosOc  9j  (2-110) 

Y  =  /m[v+l)or  +  V\;op-g  cosa  4>]  (2-111) 

Z  =   tf)  [ij  -VoO^f  as\r)QoG]  (2-112) 

L  =  Jvx  p  -  Ixz  r  (2-113) 

M  ^    Jyy  <i^  (2-114) 

M    =    Iz2  f  -  1x2  p  (2-115) 
with   auxilliary   relations: 

p=  4)  --zpsinOo  (2-116) 

^=  ®  (2-117) 

r.7i)cos0.  (2-118) 

We  find  that  the  equations  are  separable  into  two  groups. 
The  first  group  associates  the  equations  that  describe  the 
motion  in  the  xz  body  plane,  called  longitudinal  equations 
and  the  second  group  associates  the  equations  that  describe 
the  motion  in  the  xy  body  plane,  called  lateral  directional 
equations.  87 


We  find  also,  the  force  equations  divided  by  mass  and  the 
moment  equations  divided  by  the  corresponding  moment  of 
inertia. 

In  the  longitudinal  equations  the  independent  variables 
are: 

V  =  x~axis  linear  velocity 
w  =   z-axis   linear  velocity 

q  =  pitching  angular  velocity 
In    the    lateral-directional    equations    the    independent 
variables   are: 

V  =   y-axis   linear  velocity 
p  =   rolling  angular  velocity 

r  =   yawing  angular  velocity 
with    all    these    quantities    being    perturbations    from    the 
initial   flight   conditions. 

We  arrive   then  to  the   following   form  of   equations: 

Longitudinal   equations 


—  X    =    iJ  +  \A/o  q   i-  g  COS  9o  0 


(2-119) 


(2-120) 


—^    -    %  (2-121) 


Lateral -directional  equations 


^Y  =  Y.^r  -Wop  -^co50o4>  ^2-122) 
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I      ,  .         Ixz    . 

ixx  ^         Ix-x 


(2-123) 


1     M        •        lyz    . 
Izz  Iz2   '^  (2-124) 

For   further  simplification,    it   is  convenient  to   introduce 

here    the   wind   axes    system,   which  as  we  have   seen   is   produced 

by   a    rotation   about    the   y-body  axis   so   that   the  x-wind   axis 

points    in    the  opposite  direction  of   the   relative  wind.      This 

axis    is    always   tangent   to  the   flight   path   and  makes   an   angle 

of     flight    path     ^„     with    the    horizontal,    as    shown    in 

Figure  2-6. 


"ffe^fie 


boriz«i4al 


Figure   2-6 

Direction  of    the   Body  Axes  with   Respect    to  Wind  Axes    in 
Straight   Level   Perturbed  Condition 


Since 


9  -  V  +  ot 


(2-125) 
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in  the  steady  situation 

0a  =  ^o  +  c^o  (2-126) 

where 

(Xo  =0  (2-127) 

In    fact,    in    order    to   maintain   straight   level    flight,  cko 
has   to  be   small.      The   appropriate   relation  for   expressing   the 
angle    of    attack    <*    in    terms   of   the  w-velocity   if  obtained  by 
Equation   1-194    if    one    considers    small  a  and  w ,    so   that  one 
can   express: 

oi  -  ^/do  (2-128) 

From  Equation  2-127    then  we  deduce: 

^o  =  Vu.    =0  ^      Wo=0  (2-129) 

Notice  that  by  this  alignment,  since  the  wind  axes  are 
fixed  for  fixed  relative  wind,  we  can  measure  the 
perturbation  quantities   about    j^^    in   the  wind   frame. 

It  is  further  desirable  to  express  the  Lateral 
Directional  equations  in  terms  of  the  sideship  angle  |S 
instead  of  the  v-  velocity.  The  appropriate  relation  for 
expressing  the  angle  ^  in  terms  of  the  v-  velocity  is 
obtained  by  Equation  1-195,  if  one  considers  small  ^  and  v, 
so   that  one   can   express 

^  =    VUo  (2-13  0) 

So   introducing 
from   Equations  2-126    and   2-127 

0o    =  V.  (2-131) 
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from  Equation  2-129 

Wo  =0  (2-132) 

from  Equation  2-130 

v=  f!>^o  (2-133) 

we  obtain 

Longitudinal  Equations 

^  Z  =  W  -  a^  +^siny„0  (2-134) 

Lateral  Directional  Equations 
i 


-;^T  =    v  +  Ur   -gcos^„  4> 


I 


XX 


L  =  f  -  ^    r  (2-135) 


J.Z2  J-ZZ 

with   auxilliary   relations: 

p  -   i   -  tJ;  sir?  y„ 

q    ^   0  (2-136) 

r    =    li^  cos  ^„ 

The  only  thing  that  remains  is  to  express  the  left-hand 
side  forces  and  moments  in  series  of  terms  involving 
aerodynamic  derivatives  and  small  perturbation  quantities, 
that   in   fact   generate   the   forces   and  moments. 
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2.7      EXPANSION   OF    FORCES    AND   MOMENTS 

In  Section  2.6,  it  was  assumed  that  the  motion  of  the 
vehicle  consists  of  small  deviations  around  a  reference 
steady  flight   path. 

Small  disturbance  theory  applied  on  the  perturbed 
equations  of  motion  assumes  that  the  aerodynamic  forces  and 
moments  are  functions  of  the  disturbance  velocities  (angular 
and  linear),  the  control  surface  deflections  and  their 
derivatives. 

In  general,    for   a   function 

K  C:^;  fi,  V,  ...  )  (2-137) 

the  Taylor   series  expansion  about   a   reference  value 

KoC.1c,  pio,Vo  ,....  )  (2-138) 

assuming    that    the    function   has    derivatives    of    all    orders 
defined,    is  written  as 

9^ 


<=K..aK-..(|f)„a.(Ss,.(flS. 


'^'k 


im)A^'^'^^iM'''^^m''^'] 


^y 


+ 


(2-139) 


where  the  derivatives  must  be  evaluated  at  the  reference 
value  denoted  by  the  subscript  (o) . 

Now,  making  the  assumption  that  the  deviations  S^,  <^p  ,  ^ 
ov  f    are  so  small  that  product  terms  are  essentially  zero,  we 
obtain 


^^-('^^^^(lfl^.^(a^ 


(2-140) 
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Using    this   method  and  provided   that  the  above   assumption 
is   valid,    we   succeed   to  express   the  perturbation  quantity  S^^ 
as    a    linear    combination   of   the  perturbation  parameters   from 
which   it   depends   upon. 

People  in  the  field  like  to  denote  the  partial 
derivatives  as  Kj^  ,  K^  ,  Kvr  etc.,  instead  of  (c)  1^/3:1  )„  , 
(9l^/3p  )o  ,  {^^/^v  )„  ,  etc.  In  addition,  according  to  the 
convention  met   in   Section  2.6   of   dropping  the  5 ,   we   obtain: 

K^  Kj^ji  ^  K^  \i  ^  Ky  V  +■ (2-141) 

where  K^ ,    Ky  ,   Ky ,    are  called   stability  derivatives. 

Fortunately,    not    all    the    stability   derivatives    in    an 

expansion   of    a  force   or   a  moment  have   a  value.      Some   of   them 

are  zero   or   essentially  zero. 

Expansion   then   of    forces    and   moments,    in    terras    of    the 

significant   stability  derivatives  yields: 
For    longitudinal   equations: 


4rX  =  x,u  +  y^^  +Xs5  (2-142) 

^I    =   Zo^   ^Z^^  ^I^^   ^Zc^l^h^  (2-143) 

1 

J-  M   =    Mu  y  +  Mw  w  ^-  My,  w  t^0|^%  +  M^^  (2-144) 

where  it  is  clear  that 


w     1   oX       I 

^^'-^■%-  ^•■^■^-  (2-145) 
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and         denotes    the    control    input    due    to    control    surfaces 
deflections. 

For   the   lateral -directional   equations: 

;^Y   =  YvV  tYvV    tYpP  .Yrr  .YsS  (2-146) 

^L  =  L(s(^  4-  L^MLpF  ^Lrr  ^-US  (2-147) 

j^  N  =  N(j(J  +  N^^   f  Npp  +Mrr  t  Ns5  (2-14  8) 

where   it   is   clear   that 

Y     1  ~  r  -  lH 

""  '  m   "hv  '  P  ~  rw    9p  (2-149) 

Therefore   the   equations  of   motion  become: 
Longitudinal    Equations 

i;  +  Xy  u  -  X^  w    +  I  cos  jj',,  0     ^  Xg^ 

-Zyu   4-  ^  -  2,^  ^  -  Z^  ^   ''^"^  ~^%%  "^  3  ^'^(fo®    ^  ^5^  (2-150) 

Lateral  Directional  Equations 

V- Yv  V  -Y^/  -  Ypp  ^:  U,r  -Yrr  -qcos^ocj)  --  Y§  ^ 

-N(3(i  -  N^jg  -  ^""Vj^J   -  Npp  +  r  -Nrr  =  N^J  (2-151) 

with  auxilliary  relations  of  Equation  2-136. 

It  is  more  convenient  to  Laplace  Transform  the  equations 
and  arrange  them  in  matrix  form.   Further  augmenting  the 
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equations  by  the  auxilliary  relations,  using  Equation  2-130 
for  the  lateral  directional  set  and  dividing  the  Y-equation 
by  "^o    we  obtain  the   form 

Longitudinal    Equations 


5-Xy 

-2. 

0 


s-sZ^  -Zw 
-sMw  -Mw 

0 


0 
1 


8^°^^o 

u 

X6 

a^'^^ 

w/ 

25 

0 

% 

MS 

-s 

e 

0 

(2-152) 


Lateral -Directional  Equations 


6-sYv-Yv 

-V 

0 

-sL^-L^ 

s-Lp 

3xz     1 

0 

-sN^-l^ 

-jS-^p 

S-Mr                     0 

0 

0 

1 

0                     -5 

sCsm^o) 

0 

0 

i                        0 

-sCcos^o) 

here 

Yp*  =  y^/Uo 
y;  =  Yr/a 

Yi*    =   Ys/i)o 

^ 

Y/ 

? 

15 

r 

= 

N5 

i) 

0 

^ 

0 

h 


(2-153) 


(2-154) 
(2-155) 
(2-156) 

We  will  further  simplify  the  equations  of  motion  in  the 
following  chapters  in  a  form  that  is  more  suitable  for 
obtaining  aerodynamic  transfer  functions,  and  computer 
solutions. 

Our  interest  first  will  be  focused  on  examining  the 
aerodynamic  stability  derivatives.  They  will  be  left  in 
their  dimensional    form   rather   than  non-dimensional. 
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CHAPTER   3 
NON-DIMENSIONAL    STABILITY    DERIVATIVES 

3.1      INTRODUCTION 

Stability  is  a  very  important  concept  in  successful 
vehicle  flight.  We  have  encountered  the  terms  "stability" 
and  "equilibrium"  but  we  did  not  emphasize  particularly  those 
concepts. 

In  the  last  chapter  we  mathematically  obtained  equations 
of  motion  and  linearized  them  about  a  referenced  straight 
flight  which  was   assumed   to  be   an   equilibrium  condition. 

But  what  is  equilibrium?  Equilibrium  is  a  state  of  rest 
or  uniform  (unaccelerated  motion)  in  a  straight  line  where 
all   forces   and  moments   acting  on  a  body  are  balanced   out. 

Stability  is  a  property  of  an  equilibrium  state.  We  will 
consider  two  types  of  stability  which  are  of  great  interest 
in  flight  vehicle  dynamics,  defined  as  follows:  Static 
stability  is  defined  by  the  direction  and  magnitude  of  the 
initial  tendency  to  return  to  equilibrium  condition  after  a 
small  disturbance  has  occurred.  It  is  distinguished  by 
positive  static  stability  if  the  body  tends  to  return  to  its 
equilibrium  state,  negative  static  stability  if  the  body 
diverges  from  its  equilibrium  state,  and  neutral  static 
stability   if   the  body   remains   in   the  disturbed  position. 
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As  an  example,  the  pendulum  considered  in  Section  1,3  has 
a  positive  static  stability  at  its  ©"*-'''  position  but  has  a 
negative  static  stability  at  its  6=180"  position,  while  its 
motion  remains  invariant  if  the  pivot  is  moved  in  any  other 
position  and   is   said  to  have   a  neutral   static   stability. 

Dynamic  stability  determines  the  resulting  motion  in  time 
of  the  body  if  initially  disturbed.  It  is  distinguished 
similarly  in  positive,  negative,  and  neutral  dynamic 
stability. 

Various  conditions  of  dynamic  behavior  were  considered  in 
Chapter  1,  in  analyzing  spring  mass  damper  systems. 
Figure  3-1  summarizes  all  possible  time  responses  of  any 
dynamical   system. 

Dynamic  stability  implies  the  existence  of  static 
stability,   but   the   opposite   is   not  necessarily  true. 

Now  the  geometry  of  the  vehicle  comes  to  play  a  unique 
and  important  role  in  specifying  any  particular  mode  of 
dynamical  motion.  All  surfaces  exposed  to  the  airstream 
contribute  forces  and  moments  about  the  center  of  gravity, 
which  are  expressed  as  we  have  seen  in  the  last  chapter  as 
linear  combinations  of   stability  derivatives. 

The  physical  meaning  of  each  of  the  stability  derivatives 
is  straight  forward.  Thus,  Xu  is  the  rate  of  chance  of  the 
force  generated  in  the  x-direction  due  to  a  v-velocity 
perturbation,    divided  by  the  mass  of   the  vehicle,   while    Xu  o 
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Figure   3-1 
Various  Responses   of    a  Dynamical    System   to  a 
Small  Disturbance 
a,b,c  Non-Oscillatory  Modes,    d,e,f  Oscillatory  Modes 
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is   the   force   generated   in   the  x-direction  due   to  a  v-velocity 
perturbation  divided  by  the  mass  of    the  vehicle. 

In  this  chapter  we  will  give  more  insight  in  evaluating 
the  stability  derivatives  in  terms  of  the  vehicle's  geometry 
and  airstream  properties. 

3.2      LONGITUDINAL    STABILITY    DERIVATIVES 

Longitudinal  stability  derivatives  are  encountered  in  the 
longitudinal  equations  and,  hance ,  they  uniquely  specify  the 
X-force,    the   Z-force   and   the  M-pitching  moment. 

Dimensional  analysis  approach  indicates  that  we  can 
express  the  aerodynamic  forces  and  moments  in  terms  of  their 
corresponding  non-dimensional  coefficients  which  are  usually 
deduced   from  wind   tunnel   tests   of   models. 


feMivf  vMiM^ 


Figure  3-2 
Forces   and  Moments  Considered   in  Longitudinal   Motion 
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In  particular,   we  can   express  the   longitudinal    forces   and 
moments   as 


lift 


drag 


L  =  C,  ipV^£  (3-1) 


1>-C,-,f'^'^  (3-2) 


pitching  moment      M=Ct^jp^^<^  (3-3) 

where  Ci.  ,C]>  ,Cm  are  the  corresponding  lift,  drag,  and  pitching 
moment  coefficients,  p  the  air  density,  V  the  velocity  of  the 
vehicle,  S  is  a  reference  area  usually  taken  the  wing  area, 
and  c,  a  reference  length,  usually  taken  the  length  of  the 
mean   aerodynamic  chord. 

The    Longitudinal    Stability  Derivatives   can   be   evaluated 
as   follows: 

1.      The    Xt)  Derivative 

Is    the    rate    of    change    in    the    X-force    due    to    a 
u -velocity  perturbation  divided  by  the  mass  m. 

According   to  Figure   3-2,    the  X-force    is   expressed  as: 
X  =  L  siy?«  -  i)  coso^  ^-  "^  cos  Wy  (3-5) 

where  T  denotes   the   thrust. 

For    oc  small,    oc-r     essentially  zero,    and  making   use   of 
Equation  3-2   we   obtain: 
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X  .-L^r^c,  .T 


Notice  that  because  a  small  change  only  in  u-velocity 
occurs 

V^-Cllo+i;)'f  v^+w^  (3-7) 

and   since   products   of   perturbations  are  essentially   zero 

V^  =  l^J-^ZVoV  (3-8) 

where   Uo    is    the    reference   flight   forward  velocity  which   for 
practical   purposes    is   assumed  equal    to  V, 
Then   Equation  3-6   becomes 

Neglecting  the  effect  of  the  -u-velocity  on  the  thrust 
T,  by  differentiating  Equation  3-9  with  respect  to  u,  and 
dividing  by  the  mass  m,   we   obtain 


^       ^       2        'Wi  ^u  (3-10) 

Or   since  V=Uo 


Finally,    if   we  define 


"■  "  ^  [^  ^   Z   '  "bo  J  (3-11) 


^"  "  2   3(j  (3-12) 

we  arrive  at  the  following  final  form  of  the  Xu  derivative 

Au  =  -  "7^  (^^D  ^LdJ  (3-13) 
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The  quantity  Cj,^     is   essentially  zero   for    low   subsonic 
speed,     but     sometimes    near    the    critical    mach    number 
(0.8<M<1.0)     where    a    large    increase    in    drag    occurs,    C-du 
obtains  a   large  positive   value. 
2.      The   Zy   Derivative 

It  is  the  rate  of  change  in  the  Z-force,  due  to  a 
D-velocity  perturbation  divided  by  the  mass  m. 

7         L^ 

^^    "   'rn    'Bu  (3-14) 

According   to  Figure  3-2,    the   Z-force    is   expressed   as: 

Z  -  Wcos©  -  Icosoi-  -DsiYlOt  -  Tsirxxr  (3-15) 

which   is   simplified   to: 

Z=^Wcos0-L  (3-16) 

Neither    W    or   9     depend    on   't>-veloc  ity .       Since    by 

similar   approach  as   in   the  Xo  derivative,   we   obtain: 

7     -     ^(C      C    ) 

^"   '         /m      V^L^^Lo)  (3-17) 

where 

•-^  "  T     OiT  (3-18) 

The  quantity  Cl^  is  essentially  zero  for  low  subsonic 
speed,  but  sometimes  near  the  critical  mach  number  may  reach 
a  considerable  value  and  its  sign  can  be  changed  depending  on 
the  airframe  geometry,  mach  number,  dynamic  pressure  and 
aeroelastic  properties. 

For  subsonic  straight  level  flight  (j^o=o) 
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r  w 

*■  "    pa'S/2  (3-19) 

and   Zu    is   expressed   simply  as 

^^  Vo  (3-20) 

3.      The   Mv  Derivative 

It  is  the  rate  of  change  in  the  pitching  moment  M  due 
to  a  ^-velocity  perturbation,  divided  by  the  ly^  moment  of 
inertia. 

'"  "   I,,    ^u  (3-21) 

The  Mu  derivative  is  evaluated  simply  by 
differentiating   Equation  3-3   with   respect   to  u. 


^c 
where 


'^'^    "       Ihu       (^•^^^m^  (3-22) 


"   "  7     Ov  (3-23) 

The  quantity  ^My  is  very  sensitive  upon  the  variation 
of  the  u -velocity  and  its  sign  can  change  depending  upon  the 
geometry,  aeroelastic  properties,  mach  number  and  dynamic 
pressure.  Thrust  can  also  affect  '^^v  for  propeller  driven 
aircraft. 

4.      The   Xw  Derivative 

It     is    the    rate    of    change    in    X-force,    due    to    a 
w-velocity  perturbation  divided  by  the  mass  m, 
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y       1  5i     ■ 

*»   '  w    -^w  (3-24) 

A  w-velocity  perturbation  causes  a  small  change  in 
the  angle  of  attack  «.  Thus,  in  this  case  we  cannot  neglect 
the   angle   of   attach  0/  which   is   approximately  given  by: 

Vo  ^3-25) 

Therefore,    the  X-force   is  expressed  as: 

X=  L  s\noi  -Dcoso^  +T  (3-26) 


or 


y  =   -^  p  V  5  (Q  SWOC  -  Cx)  COS  c/]  +  T 


(3-27) 

Neglecting  the  effect  of  the  w-velocity  on  the  thrust 
T,  by  differentiating  Equation  3-27  with  respect  to  w  and 
dividing  by  the  mass  m,   we   obtain 

^^^^  Q7   =  -^CCLa)So^^U^sim^Q,smcy-q)^cos<y;  (3_28) 

Now  assuming   ^   small  we   express 

pus  /-  V 

Xw    =  -^    (Cl-Cd<^J  (3_29) 

where 

^"^^  "  %r  (3-30) 

The  quantity  C-^^  represents  the  change  in  drag 
coefficient  due  to  angle  of  attack  changes.  When  the  angle 
of   attack   increases,    the   total   drag   increases   and   this 
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quantity  has  normally  a  positive  value.  It  is  affected 
primarily  by  the  wings,  and  a  little  by  the  tail  and  the 
fuselage.  In  aircraft  dynamics  it  is  not  very  important 
parameter  mainly  because  the  changes  in  angle  of  attack  are 
small  and,  hence,  its  value  is  very  small. 
5.      The   Zw  Derivative 

It  is  the  rate  of  change  in  Z-force  due  to  a 
w-velocity  perturbation  divided  by  the  mass  m. 

7         15?- 

^^  '  m)    '^w  (3-31) 

The   Z-force    simillarly   is 

Z=  Wcos9  -Icosoi  -Dswoi  (3-32) 

or 

Z  =  WCDS9  -  ^pV^^CQcosoc  -Cpsincy)  ^3^33^ 

Neither  W  or  9  depend  on  w-velocity,  since  by 
differentiating  Equation  3-33  with  respect  to  w  and  dividing 
by  the  mass  m  we  obtain: 


Now,    assuming    ot  small  we   express 


(3-34) 


Zw  =  -  -^  LCl^  -^Cp)  (3-35) 


where  _ 


Cl.  = 


'<i 


-^  (3-36) 
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The  quantity  C  ^^  represents  the  change  in  C^ 
coefficient  due  to  angle  of  attack  changes.  Graphically,  it 
is  well-known  as    "lift   curve   slope",    as   shown   in   Figure  3-3. 


Figure   3-3 
Lift  Characteristic  Curves   for   Different  AR's 


The  changes  in  Q  ,  are  linear  functions  of  oe  ,  below 
the  stall.  Aspect  ratio  (Ai^),  is  defined  simply  as  the  ratio 
(span/chord)    or   in  more  general   case   as: 


k^^ 


(span)' 


\} 


Cvoin^area)  ^  (3-37) 

This  ratio,  determines  the  slope  of  the  lift  curve.  As  shown 
in  Figure  3-3,  low  A^  •  s  produce  smaller  lift  coefficient  C 
for    the   same   angle   of   attack  ^, 

The  Ci^  derivative  is  always  positive  in  the  linear 
region.  It  is  affected  primarily  by  the  wings,  while  the 
tail  and  fuselage  contribute  a  little  in  aircrafts,  but 
largely   in  missiles.      Aeroelastic  effects    in   addition  may 
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cause    bending   or   twisting   of   the  wings   that  change   the   angle 
of   attack  and  hence   the   value   of    C^^  , 
6.      The   Mw  Derivative 

Is  the  rate  of  change  in  the  pitching  moment  M,  due 
to  a  w-velocity  perturbation,  divided  by  the  lyy  moment  of 
inertia. 


lyy     ^w  (3-38) 

A  w-velocity  perturbation  is  associated  with  an  angle 
of   attack  o^  which   as  we  have   seen   is   given  by 

<y  =   VUo  (3-39) 

When  the  angle  of  attack  changes  from  its  equilibrium 
position,  the  lift  of  the  wing  changes  in  the  same  sense,  and 
that  produces  a  positive  or  negative  pitching  moment 
depending  on  the  center  of  gravity  location.  In  addition, 
the  lift  of  the  tail  changes  and  that  contributes  to  the 
total  pitching  moment.  These  contributions  together  with  the 
fuselage  contribution  make  up  this  derivative. 

We  desire  a  change  in  angle  of  attack  to  produce  a 
pitching  moment  that  tends  to  restore  the  original  angle  of 
attack.  This  requires  a  zero  pitching  moment  in  the 
undisturbed  reference  flight,  a  condition  known  as  "trim 
condition". 

The  pitching  moment,  as  we  have  seen,  is  expressed  in 
terms  of  the  non-dimentional  pitching  moment  coefficient  as 
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^-a-P^^S'^«  (3-40, 

Because  a  small  change  in  w-velocity  doesn't  affect 
the  V  velocity,  by  differentiating  Equation  3-40  with  respect 
to  w  we   obtain: 

'^w    "    ^o<  ■ '^w    "  -Uo '^   "    2Uo    f^  2ci  (3-41) 

and   finally  by   letting    V-U^    we  get 

p  S  U)  c 
^^  "  ~2%"  ^^'^  (3-42) 

The  center  of  gravity  location  plays  a  very  important 
role  in  ^^^  derivative.  In  fact  Cm^  is  proportional  to  the 
distance  between  the  center  of  gravity  and  the  aerodynamic 
center,    the   center  where   the   lift   force   acts. 

If  the  center  of  gravity  coincides  with  the 
aerodynamic  center.  Cm*  is  zero;  if  it  is  ahead  of  the 
aerodynamic  center,  C/m^  is  negative  and  the  airframe  is  said 
to  be  statically  stable;  and  if  it  is  aft  of  the  aerodynamic 
center,  C/m^  is  positive  and  the  airframe  is  said  to  be 
statically  unstable.  The  more  ahead  the  center  of  gravity  is 
located,  the  steeper  the  ^»oi  curve  becomes.  All  the  above 
cases   are   shown   in   Figure  3-4. 

It  is  clear  that  a  stable  situation  like  the  center 
of  gravity  at  B  or  A,  an  increase  in  angle  of  attack  from  a,, 
to  o'  produces  a  negative  pitching  moment  as  Figure  3-4  shows, 
which  tends   to   reduce    the   angle   of   attack. 
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Figure  3-4 
Cm„  for   Different  Center   of   Gravity  Locations 

Figure  3-4  also  shows  that  for  stable  situations,  for 
nearly  small  angles  of  attack,  where  the  lift  is  nearly  zero 
a  mechanism  must  be  provided  to  generate  a  positive  ^Mo  . 
Here  comes  the  effect  of  the  tail  which,  being  a  negative 
angle  of  attack,  generates  a  small  lift  force  which  though 
multiplied  by  a  larger  lever  arm  produces  a  balancing 
positive   pitching  moment. 

In  other  designs  we  find  the  tail  being  ahead  of 
wings  with  positive   angle   of   attack.      This   is   a   canard 
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configuration  most   usually   found   in  missiles,    also  being   the 
mechanism  to  generate   a  positive    ^Mo  . 

In    Figure    3-4      Cw^    curves    are    shown    as    linear 
functions   of    (X.      This   is  not  true   for    large   variations   in  (x, 
but   can  be  nearly  true   for   small  variations  around  (V^  • 
7.     The   Zw  Derivative 

Is  the  rate  of  change  in  Z-force,  due  to  a  rate  of 
change   in  w-velocity  perturbation  divided  by  the  mass  m 

7        L^l. 

^*^   "   w  '  QiM  (3-43) 

The   relation    o^=  w/Uo,    leads   us   to  the   expression 

'h<fi    "   l^o     Q«  (3-44) 

According  to  Equation  3-16  the  Z-force  is  simplified 
to 

Z  =  Wcos0-L  (3-45) 

Neither  lA^  or  0  depend  on  rate  of  change  in  angle  of 
attack  (i  ,  since  differentiating  Equation  3-45  with  respect  to 
o<  we  obtain 


5L  .    2L 

^«         '  "hk  (3-46) 


where  from  Equation  3-1 


5L  -  1 D  v's  5^ 

^■^  2    P  '    "^    9^  (3-47) 


Therefore, 


^"^   ^^■'3w    '  m^o   9*       "mifo    9o^    "  /mUo    2   P  "^i'  (3-48) 
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To    form   a   non-dimensional    coefficient,   we  multiply 
and    divide    Equation   3-48   by  c/2Uo   so  that  letting  V  =   Uo  we 


or 


where 


'^Vo     2  ^  2Vo      'liioic/ZVo)  (3-49) 

pSc 
^>^    "   "  "^   ^Lor  (3-50) 


''^   "    ^C<yc/2(;«)  (3-51) 


The  quantity  ^l^  represents  a  change  in  lift 
coefficient  due  to  a  variation  in  the  rate  of  change  of  the 
angle  of  attack  « . 

If  the  angle  of  attack  changes  rapidly,  the  pressure 
distribution  on  the  wing,  tail  and  fuselage,  does  not  adjust 
itself  instantaneously  to  its  equilibrium  value.  This,  a 
change  in  lift  and  drag  occurs  to  the  sudden  changes  in  angle 
of  attack. 

Its  sign  is  positive  for  low  speeds  but  can  be 
positive  or  negative  on  high  speeds  depending  on  aeroelastic 
effects. 

Also,  a  sudden  change  in  angle  of  attack  alters  the 
downwash  field  as  seen  by  the  tail  and  a  time  delay  is 
elapsed  before  the  tail  senses  this  alteration. 

All  those  effects  are  combined  to  give  rise  in  the 
^^^   derivative. 
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8.   The  Mw  Derivative 

Is  the  rate  of  change  in  the  pitching  moment  M,  due 
to  a  rate  of  change  in  w-velocity  perturbation  divided  by  the 
^yij    moment   of    inertia, 

^^   ■  In     'Bw  (3-52) 

This   derivative    '^^/Sw    can  be   expressed   as: 

-IW   ^  y«    Qi  (3-53) 

According   to  Equation  3-3 

2    ^  (3-54) 

Since  V   is   independent  of    the   rate   of  change   in   angle 

of   attack  c< ,    by  differentiating   Equation  3-54  with   respect   to 
a   we   obtain: 


^6f  2  f^  Qc^  (3-55) 


Or.  7     1  r\-. 

Therefore, 


^yy      'SW  Ujyy      ^^  Uolj^y       2    f^  '^ck  (3-56) 

To  form  a  non-dimensional  coefficient,  we  multiply 
and  divide  Equation  3-56  by  c/2Uo  so  that  letting  V  =  Uo,  we 
get 

^^    ^    7^   ^^^  (3-58) 
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where 

^^    ~  9Co.c/2Uo)  (3-59) 

For  quantity  ^m,;,  represents  a  change  in  pitching 
moment  coefficient  due  to  a  variation  in  the  rate  of  change 
of   the  angle   of   attack  «• 

The  change  in  lift  that  is  produced  by  sudden  changes 
in  angle  of  attack  discussed  already,  multiplied  by 
corresponding  lever  arms,  generates  changes  in  pitching 
moment   and   that   is  why    ^m^   arises. 

Same  lag  effects  in  the  downwash  field  discussed, 
contribute  to  this  derivative,  as  well  as  aeroelastic 
effects,  that  also  disturb  the  air  stream  of  the  wing  and 
tail,  so  that  the  pressure  distribution  over  the  aerodynamic 
surfaces   does   not  adjust   instantaneously. 

The  sign  of  this  derivative  is  negative  and  is 
relatively  important  in  longitudinal  dynamics.  Large 
absolute  values  of  this  derivative  are  desired  because  it 
causes   damping, 

9,      The   Zq  Derivative 

Is  the  rate  of  change  in  Z-force,  due  to  a  q-pitching 
velocity  perturbation  divided  by  the  mass  m 

^%^  m    0^  (3-60) 

According  to  Equation  3-16,  the  Z-force  is  simplified 
to: 

Z  =  WcosG  -L  (3-61) 
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Neither  W  or  0  depend  on  the  pitching  velocity  (q)  , 
since  differentiating  Equation  3-61  with  respect  to  q,  we 
obtain: 

'^    ^  '  \  (3-62) 

where  from   Equation  3-1 

"^l    '  2  P  -^  (3-63) 

Therefore: 

i''>*'  "}%    "'^''3^  ^''wlP^'^  2i  (3-64) 

To  form  a  non-dimensional  coefficient  we  multiply  and 
divide  Equation  3-64  by  c/2Uo  so  that  letting  V=Uo  we  get: 

H-     ^    2   ^^'^   ^^     ^Clc/2Vo  (3-65) 


or 


where 


^i      ^^   '-'■^  (3-66) 


c,    =   ^^^ 


^       '^Cc^c/2Vo)  (3-67) 

The     quantity      G^      represents    a    change    in    lift 

coefficient  due   to  a  variation   in   the   pitching  velocity. 

If    the    airframe   pitches   about    its   center   of   gravity, 

the    angle    of    attack   changes    and    a    change    in    left    of   wing, 

tail    and    fuselage    is    developed    that   gives    rise    to   this 

derivative. 
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If  has  a  small  positive  value  in  low  speed  flight  but 
may  have  large  values  positive  or  negative  in  high  speeds  due 
to  aeroelastic  effects. 

Curved  longitudinal  flight  paths  give  rise  also  to 
centrifugal  forces  on  all  components  of  the  airframe  that  due 
to  aeroelastic  effects  cause  a  bending  of  aerodynamic 
surfaces   and,    hence,    a  change   in   lift. 

In  the  past  without  computer  solutions  available, 
this  derivative  was  considered  zero  to  simplify  the  analysis. 
10.   The   Mq  Derivative 

Is  the  rate  of  change  in  the  pitching  moment  M,  due 
to  a  pitching  velocity  perturbation,  divided  by  the  lyy 
moment  of    inertia. 

^     Jaa    2%  (3-68) 

According  to  Equation  3-3, 

M  r  ipV^ScCM  (3-69) 

Since  V  is  independent  of  the  pitching  velocity  q,  by 
differentiating  Equation  3-69  with  respect  to  q,  we  obtain: 


i!l  -l-prs 


c 


H    '  ^  H  (3-70) 

Therefore, 


(3-71) 
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To  form  a  non-dimensional  coefficient,  we  multiply 
and  divide  Equation  3-71  by  c/2Uo,  so  that  letting  V=Uo  we 
get: 


%     2luu   ^  3U 


QCm 


[jj^   ^  SUo     3C^c/2i;o)  (3-72) 


or 


^  "  lliir  ^""i  (3-73) 


where 


/-  'Xm 


%       1)C^c/2a,)  (3-74) 

The  quantity  Cm<^  represents  a  change  in  pitching 
moment  coefficient  due  to  a  variation  in  the  pitching 
velocity  q. 

The  change  in  lift  that  was  produced,  due  to  pitching 
velocity  perturbations  discussed  already,  multiplied  by 
corresponding  lever  arras,  generates  changes  in  pitching 
moment  and   that   is  why    ^Mo    rises. 

The  tail  is  the  main  contributor  to  this  derivative 
which  is  also  referred  to  as  "pitching  damping  derivative". 
Aeroelastic  effects  and  centrigugal  forces  generated  due  to 
curved  longitudinal  flight  paths  also  contribute  to  this 
derivative. 

In  longitudinal  dynamics  ^f^a  is  very  important.  It 
has   a   negative   value,    and   the   larger   it    is,    the  more   the 
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damping  is  increased;  a  condition  that  is  highly  desired.  In 
high  speed  flights  Oaq  may  have  positive  or  negative  value 
depending  on  aeroelastic  effects.  It  is  also  noted  that 
tailless  airplanes   have   a  poor   damping. 

3.3      LATERAL    -   DIRECTIONAL    STABILITY    DERIVATIVES 

Lateral -Directional   stability  derivatives   are  encountered 

in    the    lateral    directional    equations    and,    hence,    they 

uniquely    specify    the   Y-force,    the    L-rolling  moment  and   the 

N-yawing  moment. 

Dimensional    analysis    approach    indicates    that   we    can 

express    the    lateral-directional    aerodynamic    forces    and 

moments    in    terms    of    their    corresponding    non-dimensional 

coefficients   as: 

Y-force  r^C^^^p'^^S:  (3-75) 

L-rolling  moment        L  =  Cg5pV'Sb  (3-76) 

N-yawing  moment  Ni  =  Q,  ^  pV^Sb  (3-77) 

where    Cy  ,    Q    ,   Q      are    the    corresponding    Y-force,    rolling 

moment    and    yawing   moment   coefficients,    and   b   is   a   reference 

length,    usually  taken   the  wing   span. 

The    Lateral-Directional    Stability   Derivatives    can    be 

evaluated  as   follows: 

1.      The  Yv  Derivative 

Is    the    rate   of   change   in  Y-force   due   to  a  v-velocity 

perturbation  divided  by  the  mass  m: 
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Y  1    ^ 

'^    "  AW     9v  (3-78) 

A  v-velocity  perturbation  causes  a  small  change  in 
the  angle  of  sideslip  ^  which  is  expressed  according  to 
Equation  2-133,    approximately  as: 

fi  =  ^/Vo  (3-79) 

so  that 

Y  —    ^ 

■  mOo    Q^  (3-80) 

Assuming  that  the  V-velocity  is  independent  of  the 
angle  of  sideslip^,  differentiating  Equation  3-75  with 
respect   to  ^ ,   we  obtain: 

2L  .  i  p  v^s  ^ 

0^    2  ^     9^  (3-81) 

Therefore,  letting  V=Uo  we  get: 

^'  ^  ~2^   ^^/i  (3-82) 

where 

<^y^  -  "^^9^  (3-83) 

The  quantity    Cy^    represents   the   change   in   the    G^    side 

force   coefficient   due    to  a  change   in   the   angle   of   sideslip  /^ . 

The    main    contributor    of    this    derivative    is    the 

vertical    tail.      Secondarily  the   fuselage,    and  also  the  wings 

contribute  a   little. 
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Figure  3-5  shows  the  forces  that  the  airframe 
experiences  due  to  a  v-velocity  perturbation.  First  an  angle 
of  sideslip  builds  up  and  all  aerodynamic  surfaces  are 
sensing  a  change   in   the   approaching   free  stream. 


Figure  3-5 
Forces   Experienced  Due   to  a  v-Velocity  Perturbation 

As  a  consequence,  a  force  F^  due  to  the  fuselage  and 
a  force  Fvi  due  to  the  vertical  tail  appears.  The  wings  also 
experience  a  different  lift  force.  For  the  configuration 
shown  in  Figure  3-5,  the  right  wing  experiences  more  lift 
than  the  left  wing  because  it  senses  a  greater  V-normal 
component  of  the  approaching  free  stream.  Those  lift  forces 
may   contribute   to  a   small  change   in   the  Y-force. 

All  these  contributions  make  a  rise  to  Cij^  derivative 
which   is  very   important   in   lateral -directional   dynamics.      It 
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introduces  a  damping  that  tends  to  eliminate  the  angle  A, 
because  the  resultant  Y-force  appears  acting  on  a  point  aft 
of  CG  location  and  that  generates  the  appropriate  moment  to 
do   this. 

Its  value  is  negative  and  it  might  be  reasonable  to 
desire  a  large  value,  for  the  ^  angle  to  die  out  quickly. 
Although  high  maneuverability  airplanes  must  have  a  small 
negative  value  for  this  derivative,  if  they  were  to  perform 
easy  roll  maneuvers  and  banked  turns. 
2.      The   Lv   Derivative 

Is    the    rate    of    change    in    L-rolling   moment  due   to  a 
v-velocity  perturbation  divided  by  the   Ixx  moment   of    inertia. 


L     2L 

Ix/  9v  (3-84) 


1  1  dL 


It  is  also  clearly  expressed  as: 

^'"'U.Iyx   Q^  (3-85) 

Differentiating  Equation  3-76,  with  respect  to  ,^,  we 
obtain 


Q^  "  2  P      9(S  (3-86) 


Therefore,  letting  V=Uo  we  get 


^'   "  ~2ir  %  (3-87) 
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where 

^/S    "     g^  (3-88) 

The    quantity     Q^     represents    the    change    in    the   Cg 
rolling   moment    coefficient    due    to   a    change   in   the   angle   of 
sideslip  ^. 

It  is  usually  referred  to  as  the  effective  Dihedral 
Derivative. 

Dihedral  is  often  used  as  a  mechanism  to  improve 
lateral-directional  stability.  Figure  3-6(a)  shows  a  head-on 
view  of  an  airplane  whose  wings  are  tilted  up  at  some 
dihedral   angle  with   respect   to  the   horizontal. 

We  have  already  discussed  that  due  to  a  v-velocity 
perturbation,  according  to  the  configuration  shown  in 
Figure  3-5,  where  a  positive  sideslip  angle  |3  is  shown,  right 
wing  generates  more  lift  than  the  left.  This  is  true  for  the 
swept  wings.  If  the  wings  are  not  swept,  the  mechanism  to 
produce  the  same  effect  is  the  dihedral.  Usually  even  swept 
wings  are  given  a  small  dihedral  to  improve 
lateral -directional   stability. 

To  understand  the  dihedral  effect,  consider  the 
aircraft  shown  in  Figure  3-7  with  non-swept  wings  but  with 
some  dihedral  angle.  The  aircraft  is  distributed  by  an  angle 
of   sideslip  ^  . 
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diMnil  <m^ 


(a) 


^  ■«_  component  o(  BrtgM 
dc^n^  h  ocBue 
nd«<({p 


V         '   'ff  Ce) 


Figure  3-6 
Dihedral  Effect  on  Lateral-Directional  Stability 


Figure  3-7 
Dihedral  Effect  on  Non-Swept  Wings 
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The  perpendicular  components  of  the  free  stream 
velocity  are  equal  but  the  angle  of  attack  that  the  right 
wing  projects  to  this  perpendicular  component  is  greater  than 
the  angle  of  attack  that  the  left  wing  projects,  as  an  effect 
of   the  dihedral   angle. 

As  a  consequence,  the  right  wing  generates  more  lift 
than  the  left  wing  as  shown  in  Figure  3-6(c).  Those  lift 
forces  contribute  to  a  change  in  the  rolling  moment 
coefficient. 

Figure  3-6(d)  and  (c)  show  the  effect  of  the 
high-wing  placement  which  has  the  same  effect  as  the  dihedral 
angle.  It  is  clearly  shown  that  in  the  low  wing  placement 
configuration,  the  fuselage  does  not  permit  the  vertical 
component  of  the  free  stream  to  act  upon  the  right  wing  to 
produce  a  stabilizing  moment,  whereas  in  the  high  wing 
placement  configuration  the  fuselage  effect  is  nearly  zero 
and   that  makes   a   stabilizing   rolling  moment   to  show  up. 

In  addition,  the  force  acting  on  the  vertical  tail 
shown  in  Figure  3-5,  multiplied  by  the  corresponding  vertical 
arm  relative  to  the  x-longitudinal  axis,  contributes  also  to 
the   same   effect. 

Cjp  is  nearly  always  negative  in  sign  and  the 
corresponding  rolling  moment  L,  is  negative  for  a  positive 
sideslip  angle  ^  ,  It  is  a  very  important  derivative  in 
lateral   directional   dynamics   because    it   causes   a  damping  and 
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improves    the    stability.      Although   high   negative    values   of 
this    derivative  will  tend   to  decrease   the  maneuverability  of 
the    airplane    in    roll   movements.      In  missiles   a   low  negative 
value   of   this  derivative   is  desired. 
3.      The   Nv  Derivative 

Is  the  rate  of  change  in  N-yawing  moment  due  to  a 
v-velocity  perturbation  divided  by  the   Izz   moment  of    inertia. 

^  '  Jzz     9v  (3-89) 

It   is   also  clearly  expressed  as: 

N   -^.^ 

'      tt.I„      op  (3-90) 

Differentiating  Equation  3-71,  with  respect  to  ^,  we 
obtain: 

-^   ^   1  ^^  ^^   ^  (3-91) 

Therefore,    letting  V=Uo  we  get: 

^'   '      Qlzz      '-■"!»  (3-92) 

where 

C-n^-^Oi/g^  (3-93) 

The  quantity  Cna  represents  the  change  in  the  Cn 
yawing  moment  coefficient  due  to  a  change  in  the  angle  of 
sideslip  ^. 
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It  is  usually  referred  to  as  the  Weathercock 
Derivative.  The  main  contributor  of  this  derivative  is  the 
vertical  tail  and  secondarily,  the  fuselage.  If  the  airframe 
experiences  a  v-velocity  perturbation  and  an  angle  of 
sideslip  ^  appears,  the  free  relative  stream  strikes  the 
fuselage  and  the  vertical  tail  in  such  a  manner  that  the 
force  ^f  ,  shown  in  Figure  3-5,  produces  a  destabilizing 
yawing  moment  and  the  force  ^vi  a  stabilizing  moment.  Wings 
contribute  a  little  in  the  positive  (stabilizing)  sense.  The 
total  value  of  this  derivative  is  positive,  signifying  static 
directional   stability. 

The  C-n^  derivative,  is  very  important  in 
lateral-directional  dynamics.  The  vertical  tail  must  be 
properly  designed  to  provide  a  sufficient  weathercock  effect. 
At  large,  sideslip  angles  a  stall  may  occur  due  to  flow 
separation  and  a  catastrophic  sideslip  divergence  may  result. 
A  low  aspect    ratio  prevents   stalling  but   decreases   stability. 

Usually  a  large  positive  value  of  the  ^^  derivative 
is  desired,  by  extending  the  vertical  tail.  A  large  value 
helps  the  pilot  in  good  coordinated  turns  and  prevents 
sideslip. 

4.      The  Yv    Derivative 

Is  the  rate  of  change  in  Y-force,  due  to  a  rate  of 
change   in  v-velocity  perturbation,    divided  by  the  mass  m. 

Y         jL    21 

^^  "  m     %  (3-94) 
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It   is   also  clearly  expressed  as: 


Vom      9^  (3-95) 

Differentiating    Equation    3-75    with    respect   to  ^  we 
obtain: 

9,3         2^-9^  (3-96) 

To  form   a  non-dimensional   coefficient  we  multiply  and 
divide   Equation  3-96   by  b/2ao. 


0|j         2   f^         2Vo    2C(i6nVo)  (3-97) 

Therefore,    letting  V=Uo  we  get: 


pSb 

4w      ^5^4  (3-98) 


where 

^/^    '  Q(:(ib/2IA,)  (3-99) 

The    quantity    Cy-     represents    the    change    in    the 

coefficient   with  variation  in   rate  of   change   of   the   sideslip 

angle   j3. 

If   the   sideslip  angle    /3  changes    rapidly,    the  pressure 

distribution    on    the    vertical    tail    and    fuselage    does    not 

adjust   itself   instantaneously  to   its  equilibrium  value. 

This    effect,    as    well    as    aeroelastic    effects,    give 

rise    in     ^y-j      derivative    but    its    value   is   very   small   and   is 

usually  neglected. 
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5.     The  Ly  Derivative 

Is  the  rate  of  change  in  L-rolling  moment,  due  to  a 
rate  of  change  in  v-velocity  perturbation  divided  by  the  Ixx 
moment   of    inertia. 


where 


Ix*      Qv  (3-100) 


I  •    -  J_ 


obtain: 


It   is   also  clearly  expressed  as: 

^'  =  aix/  ^  (3-101) 

Differentiating   Equation  3-76,   with   respect   to  /3,   we 


^  ^  rP^  ^  i;r  (3-102) 

To   form   a  non-dimensional   coefficient 

0^     ~   2   ^  SVo      "^C^hnVo)  (3-103) 

Therefore,    letting  V=Uo 

I-  -1^  C 

^^   ~    -^Ixx     ^^  (3-104) 


'^a  ~  QC/^b/2yo)  (3-105) 

The  quantity  Q:  represents  the  change  in  the  C^ 
rolling  moment  coefficient  with  variation  in  rate  of  change 
of   the   sideslip  angle    &, 
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The  change  in  the  side  forces  on  the  vertical  tail 
and  the  fuselage  produced  due  to  sudden  changes  in  the 
sideslip  angle  /3 ,  discussed  already,  multiplied  by 
corresponding  lever  arms  provide  a  total  change  in  the  Q 
rolling  moment  coefficient. 

This    derivative    has    a   very    small   value   because   the 
lever   arms   are  very   small   and   it   is   usually  neglected, 
6,     The  Ny   Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment,  due  to 
a  rate  of  change  in  v-velocity  perturbation  divided  by  the 
Izz  moment  of    inertia. 

N^  -  i;;     9v  (3-106) 

It   is   also   clearly  expressed  as: 

^'   "  ^llz'  3(S  (3-107) 

Differentiating  Equation  2-71,  with  respect  to  ^,  we 
obtain: 

^=2P^^^^  (3-108) 

To   form  a  non-dimensional   coefficient: 

^    =     2    P  ^    ^  ''   m       W^^jm  (3-109) 

Therefore,    letting  V=Uo  we  get: 

Nv  =  Ti^r  ^'^^  (3-110) 
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where 

"^     "     9C/2b/2a)  (3-111) 

The    quantity    C^ha     represents    the    change    in    the   G» 
yawing  moment  coefficient  with  variation   in   rate  of   change   of 
the   sideslip  angle  ^. 

The  change  in  the  side  forces  on  the  vertical  tail 
and  the  fuselage  produced  due  to  sudden  changes  in  the 
sideslip  angle  /^  ,  discussed  already,  multiplied  by 
corresponding  lever  arms  provide  a  total  change  in  the  0» 
yawing  moment  coefficient. 

The    lever    arm  of    the  vertical   tail   this   time   is   long 
and    hence f    this    derivative   may   be    considerable.      Positive 
values    of    this   derivative    increase   the  damping   in   the  yawing 
motion  and   are  desirable   in  general. 
7,     The  Yp    Derivative 

Is  the  rate  of  change  in  the  Y-force,  due  to  a 
p-rolling  velocity  perturbation  divided  by  the  mass  m. 

Y    -  i-    51 

'P  "  ^     9p  (3-112) 

Differentiating  Equation  3-75,  with  respect  to  p,  we 

obtain: 


Op  ^  2  P      2p  (3-113) 
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To   form  a  non-dimensional   coefficient 


i^^P^'s 


<9p        2  ^  2Vo      2Cpbl2Uo)  (3-114) 

Therefore,    letting  V=Uo 

^    =        4^        ^yp  (3-115) 


where 


c     -        "^^^ 

yP  ~    QCpb/2yo)  (3-116) 

The    quantity    Cy^     represents    the    change    in    the    Cy 
Y-force   coefficient  with  variation   in   the   rolling  velocity  p. 

When  an  aircraft  has  a  roll  velocity  p,  a  force  is 
acting  on  the  vertical  tail  that  usually  opposes  the  roll 
movement.  This  consists  the  main  contribution  to  this 
derivative   and   is   small   enough  to  be   neglected. 

The    force    acting  on  the  vertical    fin  may  be   affected 
by  the   sidewash   from   the  wing   that  gives   a   slightly  different 
direction  of   the   airstream. 
8,      The    Lp    Derivative 

Is  the  rate  of  change  in  the  L-rolling  moment,  due  to 
a  rolling  velocity  p  perturbation  divided  by  the  Ixx  moment 
of    inertia. 
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Differentiating  Equation  3-76,  with  respect  to  p,  we 
obtain: 

'^  =  1  p  V'  S  b  ^ 

Qp   2  f^       9p  (3-118) 

To   form  a  non-dimensional   coefficientf 

r)p    ~  2   r  20o      '3Cpb/2a,)  (3-119) 

Therefore,    letting  V=Uo 

P    "       mxx        ^f  (3-120) 

where 


Q.    =         ^ 


P        QCpb/2U,)  (3-121) 

The  quantity  Qp  represents  the  change  in  Q  rolling 
moment  coefficient  with  variation  in   the   rolling  velocity  p. 

This  derivative  is  also  referred  to  as  "roll  damping 
derivative",  because  it  has  usually  a  negative  or  damping 
effect  coming  mainly  from  the  wings  and  secondarily,  from  the 
vertical   and  horizontal   tails. 

When  the  aircraft  rolls,  the  force  acting  on  the 
vertical  tail,  shown  in  Figure  3-8,  multiplied  by  the 
corresponding  lever  arm,  produces  an  opposing  rolling  moment. 
In  addition,  on  the  down-going  wing,  the  effective  angle  of 
attack  is  increased  and  on  the  up-going  wing,  decreased. 
This  makes  the  lift  on  the  down-going  wing  greater  and  thus 
an  opposing  rolling  moment  is  produced.  Similar  effect,  on  a 
smaller   scale,    is   produced  by  the   horizontal   tails. 
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Of-^o\nq  mn^ 


Figure  3-8 
Forces   Due   to  Rolling  Velocity    (p) 

This  is  true  as  long  as  the  down-going  wing  is  not 
flying  near  the  stall  due  to  a  large  angle  of  attack 
experienced  in  some  cases.  In  such  cases  a  positive  effect 
will   result   in   autorotation. 

The  derivative  ^?p  is  quite  important  in 
lateral-directional  dynamics.  The  value  of  this  derivative 
directly  affects  the  ailerons  design.  A  small  negative  value 
of   this   derivative  will   result   in   a  better   aileron  input 
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response    while   a   large  negative   value  will  make   the  airframe 
heavily  damped   in   roll  maneuvers. 
9.      The    Np   Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment,  due  to 
a  rolling  velocity  p  perturbation,  divided  by  the  Izz  moment 
of   inertia. 

M    -  i_.5!l 

P  "  IzE      9p  (3-122) 

Differentiating  Equation  3-77,  with  respect  to  p,  we 
obtain: 

<^    "    2    P  2y  (3-123) 

To   form  a  non-dimensional   coefficient, 

—   -  i  oV^Sb    —  ^^^ 

9^    "  2   r  2i;  '  'dC'phJ2U>')  (3-124) 

Therefore,    letting  V=Uo 


where 


^P   ■   ~7l^   ^^  (3-125) 


U.  =         ^^"^ 


?        0Cpb/2i;<,)  (3-126) 

The   quantity    Crip    represents   the   change   in  Cn    yawing 

moment  coefficient  with  variation  in  the  foiling  velocity  p. 
This    derivative    usually   has    a    negative   or   positive 

effect     coming     from    the    wings,     the    vertical    tail    and 

secondarily,    from   the   horizontal   tail. 
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When  the  aircraft  rolls,  the  force  acting  on  the 
vertical  tail,  shown  in  Figure  3-8,  multiplied  by  the 
corresponding  lever  arm  produces  a  yawing  moment  towards  the 
same  side  of  the  roll  motion,  i.e.,  roll  to  starboard  causes 
yaw  to  starboard.  This  contributes  a  positive  effect  and 
depands  upon  the  vertical  tail  geometry  and  the  sidewash  from 
the  wings. 

In  addition,  the  contribution  of  the  wings  can  be 
considered  to  be  made  of  two  parts,  one  opposite  the  other, 
the  total  effect  being  negative.  When  the  airframe  rolls,  on 
the  down-going  wing,  the  effective  angle  of  attack  increases 
and  on  the  up-going  wing  decreases,  as  shown  in  Figure  3-8. 
As  a  consequence,  the  lift  and  drag  are  greater  on  the 
down-going  wing  than  on  the  up-going  wing.  The  projections 
of  the  drag  forces  on  the  xoy  plane  usually  produce  a  yawing 
moment  towards  the  same  side  of  the  roll  motion.  This 
produces  a  positive  effect  while  the  projections  of  the  lift 
forces  produce  a  yawing  moment  towards  the  opposite  side  of 
the  roll  motion,  and  this  produces  a  negative  or  damping 
effect. 

The  overall  value  of  this  derivative  is  usually 
negative,  i.e.  roll  to  starboard  produces  a  yaw  to  port.  At 
higher  speeds  the  numerical  value  decreases  and  may  become 
positive. 
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This  derivative  is  fairly  important  in 
lateral-directional  dynamics,  as  far  as  roll  maneuvers  and 
their  effect   to  yaw. 

10.   The   Yr    Derivative 

Is  the  rate  of  change  in  the  Y-side  force,  due  to  a 
yawing   velocity  perturbation  r,    divided  by  the  mass  m. 

Y   -  1   ^ 

"  ^  -w     9r  (3-127) 

Differentiating  Equation  3-75,  with  respect  to  r,  we 
obtain: 

"^r        Z    ^  9r-  (3-128) 

To  form  a  non-dimensional   coefficient. 


5l  =  ioV^S  —        ^^^ 

9r        2   "  20,      '3Crb/2U.)  (3-129) 


Therefore,  letting  V=Uo 


where 


^'-^  — JT"  ^yr  (3-130) 


QCrb/2a,)  (3-131) 

The    quantity    ^y^     represents   the   change   in   ^^    Y-force 

coefficient  with   variation   in   the  yawing  velocity   r. 

When   the   aircraft   has   a  yaw  velocity   r,    as   Figure  3-9 

shows,    a   force    is   acting   on  the   fuselage   that  usually  has   a 
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negative    or    damping    effect   on  the  yaw  motion  and   a   force    is 
acting   on  the  vertical   tail   that  usually  has   the   same   effect. 


Figure  3-9 
Forces  Due   to  Yaw  Velocity    (r) 


Both   contributions  are  small  enough  to  be   neglected, 
therefore,    this    derivative    is    of    little    importance    in 
lateral-directional   dynamics. 
11.   The    Lp   Derivative 

Is  the  rate  of  change  in  the  L-rolling  moment,  due  to 
a  yawing  velocity  r  perturbation,  divided  by  the  Ixx  moment 
of    inertia. 


(3-132) 
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Differentiating   Equation  3-76,   with   respect    to   r,   we 
obtain: 


QL 


To  form  a  non-dimensional   coefficient. 


(3-133) 


Qr    "  2    P^      2Uo       'DCrb/2y.;  (3-134) 


Therefore,    letting  V=Uo 


4lxx         'r  (3-135) 


where 


'^Crb/2a)  (3-136) 

The  quantity  Qr  represents  the  change  in  Q  rolling 
moment  coefficient  with  a  variation  in  the  yawing  velocity  r. 
When  the  aircraft  is  yawing  about  its  vertical  axis, 
the  outer  wing  moves  faster  than  the  inner.  As  a 
consequence,  the  lift  on  the  outer  wing  is  greater  than  that 
of  the  inner  and,  hence,  a  rolling  moment  is  produced  towards 
the  same  side  of  the  yaw  motion.  This  contributes  a  positive 
effect,    i.e.    yaw   to   starboard  produces    roll  to  starboard. 

In  addition  to  this  major  wing  contribution,  the 
force  acting  on  the  vertical  tail  multiplied  by  the 
corresponding   lever   arm,    produces   the   same    rolling   effect. 
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The  total  value  of  this  derivative  is  usually 
positive,  but  it  is  not  so  important  in  lateral -directional 
dynamics. 

12.   The    ^r   Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment,  due  to 
a  yawing  velocity  r  perturbation,  divided  by  the  Izz  moment 
of    inertia. 

'  "  3z2       2)r  (3-137) 

Differentiating  Equation  3-77,  with  respect  to  r,  we 
obtain: 

"hr     '  2   ^  ^r  (3-138) 

To  form  a  non-dimensional   coefficient, 

^r         2    ^  2a      '^Crb/2Vo)  (3-139) 

Therefore,    letting  V=Uo 

^'  ■         4lzz  ■-  (3-140) 

where 

^  c)Cri 

^-^  ■'    ^Crb/2Uo)  (3-141) 

The    quantity    6r7^    represents   the   change   in    C^,   yawing 

moment    coefficient,     with    a    variation    in    the    yawing 

velocity   r. 
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This  derivative  is  known  as  the  "yaw  damping 
derivative",    because    it   usually  has   a   negative   effect. 

The  main  contributor  is  the  vertical  tail.  When  the 
aircraft  is  yawing,  the  force  acting  on  the  vertical  tail 
multiplied  by  the  corresponding  long  lever  arm,  produces  a 
considerable  yawing  moment  that  opposes  the  initial  yawing 
disturbance. 

In  addition,  the  drag  on  the  outer  wing  is  greater 
than  that  of  the  inner  wing.  That  produces  a  yawing  moment 
that  opposes  the  initial  yawing  disturbance.  The  force  on 
the   fuselage   also   contributes   an   opposing  effect. 

This  derivative  is  very  important  in 
lateral-directional  dynamics.  Its  total  effect  is  negative 
and   large  values   are  desired  for   more  effective   damping. 

3.4      CONTROL    SURFACE    DEFLECTIONS 

We  have  discussed  so  far  all  the  important  stability 
derivatives  that  rise  due  to  the  airframe  geometry.  They 
uniquely  specify  the  forces  and  moments  about  the  center  of 
gravity  provided  that  they  are  multiplied  by  the  correct 
factor  that  depends  upon  the  airstream  and  airframe 
characteristics. 

A  stability  analysis  is  possible  now  for  determining  the 
longitudinal  and  lateral  directional  modes  of  motion, 
corresponding  to  an  initial  disturbance,  i.e.  homogeneous 
solution  to  the  differential   equations. 
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In  addition,  some  means  of  control  in  longitudinal  and 
lateral  directional  motion  are  available  to  the  human  or 
automatic  pilot.      They  are   shown   in   Figure  3-10. 


elerciW     cot  Vol 


rnddtr  cm^ni 


Figure  3-10 
Basic  Longitudinal  and  Lateral  Control  Systems 
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Elevator,  controls  the  longitudinal  motion  of  the 
airframe  providing  nose  up  or  nose  down  moments  as  a  response 
to  corresponding   control   surface   deflections. 

Aileron  and  rudder,  control  the  lateral  directional 
motion  providing  rolling  and  yawing  moments  correspondingly 
as    a    response    to    the    corresponding    control    surfaces 

deflections. 

The  symbols  used  for  the  corresponding  angles  of  rotation 
and   control   surface   deflections  are  shown  below: 


axis 

Ox 

Oy 

Oz 

aYig?e    of  roiailon 

4) 

e 

(pikh) 

IP 

(yaw) 

COTilrol   sur-face   olefl. 

(ai  IgKon) 

(eievaior) 

5 

C  rwdder) 

In  the  longitudinal  and  lateral  directional  equations  of 
motion,  Equations  2-152  and  2-153,  respectively,  we  have 
denoted  by  the  symbol  <^  the  control  input  due  to  control 
surfaces   deflections. 

The   right-hand   side   of   the   longitudinal    Equation  2-152    is 

expanded  as  follows: 

Zs    =  Zrx 

M5  =    Mti  +  Mti 


(3-142) 
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and    the    right-hand    side    of     the    lateral    directional 
Equation  2-153    as: 

Y5  =  Ys 

U  =   Li  4.Lg   Us  (3-143) 

Ns  =  Ng  ^.  N5  +  N5 

In  addition,  to  the  right-hand  side  terms  that  take  care 
of  the  main  control  surfaces,  there  may  exist  terms  rising 
from  the  auxiliary  surfaces  such  as  flaps,  stabilizers,  dive 
brakes,    etc. 

A  general  approach  suggests  to  examine  the  subscripted 
derivatives  instead  of  the  particular  derivatives  dedicated 
for  a  specific  control  surface  deflection,  i.e.  the  total 
effect  by  all  the  possible  combinations  of  the  control 
surface  deflections.  The  approach  is  similar  to  that 
followed   in   Sections  3.2   and   3.3. 

1.      The    Xs    Derivative 

Is    the    rate    of    change    in    the    x-force    due    to   a 
control   surface   perturbation  divided  by  the  mass  m. 

^       m     'BS  (3-144) 

By   5    we    denote  a  general   control   surface   deflection. 

For    longitudinal    control    it   might    be    elevator,    stabilizer, 

flap,    slat,    dive  brakes,    etc. 

According  to  Equation  3-27,    the  X-force    is   simplified 

to: 
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We  can  assume  that  neither  V  or  T  depends  on  the 
control  deflection.   Thus  differentiating  Equation  3-145, 
with  respect  to  ,   we  obtain: 


'^a     2  t^   "  96  (3-146) 

If  we  let  V=Uo  and 

Cds   =   '3Q>/9<5  (3-147) 

we  arrive  at  the  expression: 

^  2rM        ^^5  (3-148) 

The  quantity  w?s  represents  a  change  in  the  drag 
coefficient  due  to  a  variation  in  a  longitudinal  control 
surface  deflection. 

A  control   deflection,    positive   or    negative,    increases 
the    drag    due    to  the  greater  projected   frontal   area,    but   this 
sometimes    is    desirable   because    it   provides    damping    in    an 
oscillatory  perturbed  motion. 
2.      The   2h    Derivative 

Is  the  rate  of  change  in  the  Z-force  due  to  a  § 
control   surface   perturbation  divided  by  the  mass  m. 

7  -  i_  5L 

^  '   'W     9S  (3-149) 

According  to  Equation  3-16,  the  Z-force  is  simplified 
to: 

Z  =  Wco5  9  -L  (3-150) 
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Neither  W  or  9  depends  on  the  ^  control  deflection. 
Thus  differentiating  Equation  3-151,  with  respect  to  ^  we 
obtain: 

^h    '       2  P  9S  (3-152) 

If  we   let  V=Uo  and 

Q5  =  ^^/05  (3-153) 

we  arrive  at  the  expression 

7         l^r 

^  '  Sm      ^'■<5  (3-154) 

The  quantity  Q^  represents  a  change  in  the  lift 
coefficient  due  to  a  variation  in  a  longitudinal  control 
surface  deflection. 

This  derivative  is  a  measure  of  the  effectiveness  of 
the  elevator  mainly,  in  producing  a  change  in  the  lift  force 
acting  on  the  tailplane.  A  positive  elevator  deflection 
produces   an   increase   of    lift,    i.e.    negative   Z-force. 

It    has    a    very   small  value   on  aircrafts  with   tail   but 
a   large  value   on  tailless  aircrafts. 
3.      The    Ms  Derivative 

Is  the  rate  of  change  in  the  M-pitching  moment  due  to 
a  control  surface  perturbation  divided  by  the  lyy  moment  of 
inertia 
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According  to  Equation  3-3 ^  the  pitching  moment  is 
expressed  as: 

M=  jpV-^ScC^  (3-156) 

Differentiating  Equation  3-156,  with  respect  to  S  we 
obtain: 

Q5    2  t"      2S  (3-157) 

If  we  let  V=Uo  and 

CMg   .-   "^^^/dS  (3-158) 

we  arrive   at  the   expression 

'^^  "       5Iyy       ^'"^  (3-15  9) 

The  quantity  Cw§  represents  a  change  in  the  pitching 
moment  coefficient  due  to  a  variation  in  a  longitudinal 
control   surface   deflection. 

It  is  commonly  referred  to  as  the  elevator  power 
derivative  because  it  is  a  measure  of  the  elevator  control 
effectiveness. 

The  change  in  Z-force  produced  by  the  elevator 
control  surface  deflection,  multiplied  by  a  rather  long  arm, 
provides  a  pitching  moment  in  the  negative  sense,  i.e.  a 
positive  elevator  deflection  results  in  a  nose-down 
(negative)    pitching  moment. 

That  is  the  primary  function  of  the  elevator  and  is 
considered   to  be   the  most   important  of   all  the   control 
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surface  functions.  It  adjusts  the  angle  of  attack  and  it 
trims  the  vehicle  in  its  longitudinal  motion,  by  providing 
opposite  pitching  moments  to  accommodate  disturbing  moments 
and  center   of   gravity  relocations. 

For    a   design  point  of   view,    the   larger   the  center   of 
gravity    range    is,    the    greater   the  elevator   pitching  moments 
must  be  provided,    i.e.    larger    Ci^^   values. 
4.     The  Y5    Derivative 

Is    the   change   in  the  Y-side   force   due   to  a    5  control 
surface   perturbation  divided  by  the  mass,   m. 


OY 
Differentiating    Equation  3-75,   with   respect   to  S   we 


^  "  ^   -as  (3-160) 


obtain: 


%         2  ^  7)h  (3-161) 


If  we   let  V=Do  and 

C|,5    --   ^Q/Q§  (3-162) 

we  arrive   at  the   expression: 

^    '        S-rvt       ^^5  (3-163) 

The  quantity    ^yj    represents   the   rate  of  change   in   the 

y-side    force    coefficient    due    to    a   variation  in    a    lateral 
control   surface   deflection. 
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This  derivative  rises  mainly  from  the  rudder,  A 
positive  rudder  deflection  is  defined  to  produce  a  positive 
side  force.  Therefore,  the  value  of  this  derivative  is 
positive  but   small  enough   to  be   neglected. 

Rudder    deflections   also    cause    an    increase    of    the 
total  drag  of   the  vehicle  but   that  has   a  negligible  effect   on 
the   longitudinal   motion,    and   is  neglected, 
5.     The   U    Derivative 

Is  the  rate  of  change  in  the  L-rolling  moment  due  to 
a  0  control  surface  perturbation,  divided  by  the  Ixx  moment 
of    inertia.         * 

^       Ixx     dh  (3-164) 

Differentiating  Equation  3-76,  with  respect  to  5,  we 
obtain: 

^  "  TP^    ^*^    IT  (3-165) 

If  we  let  V=Uo  and 

Q^  .   QQ/9^  (3-166) 

we  arrive   at  the   expression: 

/^    ^        2Ixx        ^^^  (3-167) 

The    quantity    Qj     represents    the    change    in    the  Q 
rolling   moment    coefficient    due    to   a   variation   in   a   lateral 
control   surface   deflection. 
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The  side  force  generated  by  a  rudder  deflection, 
produces  a  rolling  moment  in  the  positive  sense,  i.e.  a 
positive  rudder  deflection  generates  a  positive  rolling 
moment.  In  airplanes,  with  large  vertical  tails,  this  effect 
is  negligible  and   is  often  neglected. 

The  main  contribution  to  this  derivative  is  from  the 
ailerons,  and  it  is  a  measure  of  the  aileron  effectiveness. 
A  positive  aileron  deflection  is  defined  as  the  one  to 
produce  a  positive  rolling  moment,  i.e.  left  aileron  down  and 
right  aileron  up,  therefore,  the  value  of  this  derivative  is 
positive. 

In  lateral-directional  dynamics,  this  derivative  is 
probably  the  most  important  because  it  establishes  the 
maximum  available  rate  of  roll  of  an  airplane  which  is  a  very 
important   consideration  in   combat  environment. 

It  also  plays  an  important  role  in  take-offs,  and 
landings,  where  adequate  and  rapid  lateral  control  is  highly 
desirable. 

Rapid    rate    of    change    of    aileron   deflections   gives 
also    a    rise    to   this    derivative   due   to  lag  effects,    i.e.    not 
immediate  adjustment   of   the   airstream. 
6.     The    Ms   Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment  due  to  a 
d  control  surface  perturbation,  divided  by  the  Izz  moment  of 
inertia. 
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Differentiating   Equation  Z-11 ,   with   respect   to  5  ,   we 


^^   '  Izz      3&  (3-168) 


obtain: 


'35     "  J  ^  'as  (3-169) 

If  we  let  V=Uo 

Crjg  -  9C^/9S  (3-170) 

we  arrive  at  the  expression: 


^v:s 


c 


^^  "   —JuT  "-"J  (3-171) 

The    quantity    Cn^     represents    the    change    in    the  G 
yawing   moment   coefficient   due    to   a   variation   in    a  lateral 
control   surface   deflection. 

The  side  force  generated  by  a  rudder  deflection, 
produces  a  yawing  moment  in  the  negative  sense,  i.e.  a 
positive  or  to  the  left  rudder  deflection,  produces  a 
negative  yawing  moment.  This  derivative  is  a  measure  of  the 
rudder   effectiveness. 

Rapid  rate  of  change  of  rudder  deflections  gives  also 
a   rise   to  this   derivative   due    to  lag  effects. 

This  derivative  plays  an  important  role  in  crosswind 
takeoffs,  and  landings,  and  counteracts  adverse  yaw  in 
rolling  maneuvers. 
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Finally,  aileron  deflections  give  also  a  rise  to  this 
derivative  since  each  aileron  introduces  a  drag  force  of 
different  magnitude.  This  effect  may  be  positive  or  negative 
depanding  on  the  rigging  of  the  aileron  and  the  angle  of 
attack  of   the   airframe. 

3.5      DETERMINATION   OF   THE    STABILITY    DERIVATIVES   AND   OTHER 
CONSIDERATIONS 

All  the  stability  derivatives  discussed  so  far  rise  from 
the  flight  vehicle's  geometry,  being  highly  dependent  upon 
the  design  layout.  We  have  referred  mostly  in  this  chapter 
to  the  conventional  airplane  with  horizontal  tail  and  we  have 
seen  that  the  derivatives  are  made  up  of  a  number  of  separate 
contributions. 

We  had  strictly  remained  in  the  dimensional  form  of  the 
stability  derivatives  and  we  simply  expressed  them  in  terms 
of  airstream  characteristics,  airframe  geometry  and 
non-dimensional   derivatives. 

Dimensional  form  has  the  advantage  of  obtaining  real  time 
solutions  rather  than  working  with  normalized  time  mass  and 
inertia  components  that  must  also  be  non-dimensional  in 
considering   that  form. 

The  advantage  of  non-dimensional  derivatives,  i.e. 
independent  of  size  and  airstream  characteristics  is 
maintained  because   dimensional   derivatives   are  expressed  with 
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the    corresponding    non-dimensionalSf    some   of   which   though  do 
depend  upon  the   airstream  and   airframe  characteristics. 

Thus,  some  of  the  dimensional  expressions  and  some  of  the 
non-dimensionals  are  constants,  but  they  can  be  expressed  in 
polynomial   form. 

As   an   example,    for    the  METEOR  7: 

Qy    =  -0.052  -   0.009  C, 
C^r    =   0.0295   +   0.159  C, 

Crip    =   -0.042  Cl 
Cr^   =   -0.182   -    0.052    C/ 

For  others,  plots  may  exist  indicating  their  variation 
with   respect   to    Q    or   mach  number. 

During  the  development  of  an  aircraft,  it  is  important  to 
know  the  numerical  values  of  the  derivatives  under  various 
conditions  of  flight.  In  the  first  stage,  theoretical 
estimates  can  be  made  supported  by  wind  tunnel  tests  and 
model   experiments. 

By  then,  simulation  in  digital  computers  can  be  made  to 
obtain  response  characteristics  to  initial  disturbances  or 
control  deflections  and  determine  whether  the  values  of  the 
derivatives   are  satisfactory. 

When  the  aircraft  has  been  built,  flight  tests  are 
carried  out  to  confirm  the  theoretical  with  the  actual 
responses. 
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There  also  exist  flight  test  techniques  for  determining 
the  values   of    the  derivatives. 

Tables  3.1  and  3.2  summarize  the  longitudinal  and  lateral 
directional  derivatives  including  definition,  unit, 
expression  and   auxilliary   relations. 
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TABLE  3-1 
LONGITUDINAL  STABILITY  DERIVATIVES 
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TABLE  3-2 
LATERAL  DIRECTIONAL  STABILITY  DERIVATIVES 
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CHAPTER  4 
LONGITUDINAL  DYNAMICS 

4.1   INTRODUCTION 

Having  specified  all  the  vehicle's  dynamic  properties,  we 
will  now  examine  the  dynamic  motion,  i.e.  the  vehicle's 
response   to  disturbances   and   control    inputs. 

Since  any  flight  vehicle  is  a  dynamic  system,  described 
by  differential  equations  of  time,  we  must  expect  in  a 
particular  situation,  oscillatory  or  exponential  converging 
or   diverging   response. 

Any  particular  response  will  be  associated  by  its 
corresponding  characteristic  root  location  in  the  complex 
plane,  as  we  have  clearly  seen  in  the  dynamic  response  of  the 
spring  mass  damper   system. 

The  characteristic  polynomials  in  the  Laplace  variable 
(s)  involved,  are  usually  of  the  third  degree  or  higher  and 
thus  hard  to  solve  by  hand. 

In  the  past,  a  lot  of  approximations  were  introduced  and 
even  reduction  of  airframe's  degrees  of  freedom  technique  was 
required  to  analyze  coarsely  a  particular  mode  of  motion,  by 
reducing  in  effect  the  degree  of  the  characteristic 
polynomial. 

In  this  chapter,  we  will  analyze  the  longitudinal  motion 
of    a    conventional    airplane.       We   will    first   discuss    the 
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Laplace  Transform  solution  by  obtaining  the  roots  of  the 
associated  characteristic  polynomial.  Then  we  will  apply  the 
reduction  of  degrees  of  freedom  technique  to  isolate  and 
analyze  a  particular  mode  of  motion.  Finally,  the  equations 
of  motion  will  be  rearranged  in  a  form  that  is  more  suitable 
for  computer  solutions. 

As  an  example,  the  longitudinal  dynamic  motion  of  the 
aircraft  shown  in  Figure  4-1  will  be  analyzed  with  the 
parameters  indicated. 

4.2   CHARACTERISTIC  POLYNOMIAL  OF  LONGITUDINAL  MOTION 

Longitudinal  equations  of  motion  are  farther  simplified 
by  neglecting  the  effect  of  the  Zq  and  Zw  derivatives. 
Experimental  results  have  shown  that  the  effect  of  neglecting 
them  is  only  of  secondary  importance.  Further  a  straight 
level  flight  is  assumed  as  the  reference  flight  condition 
(i.e.  }^<7  =0)  ,  hence.  Equation  2-152  can  be  reduced  to: 


S-Xy 

-Xw 

0 

^ 

u 

Xl 

-Iv 

S-Zw 

-Uo 

0 

w 

21 

-My 

-sMw-Mw 

s-M^ 

0 

% 

MS 

0 

0 

i 

-5 

6 

0 

(4-1) 


The  characteristic  polynomial  is  found  by  evaluation  of 
the  above  determinant.'  The  Equation  4-1  can  be  reduced  in  3 
equations  of  the  original  y,w,q,  longitudinal  variables  by 
avoiding  to  write  down  the  definition  equation  for  q,  but 
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Figure  4-1 

Conventional  Airplane  and  Parameters 
Used  for  the  Numerical  Example 
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that  doesn't   simplify   the   hand   calculations.      The   detenninant 
comes   out   to  be : 

As^  +Bs^  +Cs^  +  Ds  +  E    =0  (4-2) 

where 

A-I  (4-.3) 

'%'-Ha-VoMii-Zw-yo  (4-4) 

C  =  ZwM^  -UoMw  -XvuZu  ^-XoCM^+y„MW  fZw)  (4-5) 

D^-Xu  CZ«>,M(|^-UMw)+2oCXwM(^f3Mw)-Mw(i;oXv«-|)  (4-6) 

E  =  I  (ZuMw  -MyZw)  (4-7) 

Equation  4-2  is  referred  to  as  the  longitudinal  stability 
quartic.  The  roots  in  this  case  for  any  conventional  vehicle 
come  out  to  be  complex  with  negative  real  partSr  which  are 
associated  with   two  oscillatory   convergent  motions. 

The  overall  longitudinal  motion  of  the  vehicle  is  a 
superposition  or  a  combination  of  these  two  oscillatory 
modes. 

One  of  the  modes  usually  is  associated  with  a  pair  of 
complex  roots  located  near  the  origin  of  the  negative  complex 
plane.  This  mode  is  called  "Phugoid  motion"  and  is  a  fairly 
lightly  damped  oscillation  with   long  periodic  time. 

The  other  mode  usually  is  associated  with  a  pair  of 
complex  roots  located  far  from  the  origin  of  the  negative 
complex  plane.  This  mode  is  called  "Short  period  motion"  and 
is  a  heavily  damped  oscillation  with  very  short  periodic 
time. 
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The  longitudinal  oscillatory  modes  are  shown  in 
Figure  4-2. 


/ 


Figure  4-2 
The   Longitudinal   Oscillatory  Modes 


For    the  given  example  we   find 
A  =   1.0 
B  =   4.2177 
C   =   18.2962 
D   =    0.1814 
E   =    0.0722 
and   the   characteristic  polynomial    is 

6^  +  4.2177  s^  ^  ^82962  s' +  0.1814  s  f  0.0722   =0 

which  has    roots 

-0.0045  ±  j    0.0627 
-2.1043   ±  j    3.7184 


(4-8) 

(4-9) 

(4-10) 

(4-11) 

(4-12) 

(4-13) 

(4-14) 
(4-15) 
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The  pair  of  roots  -0.0045  ±  j  0.0627  correspond  to  the 
phugoid  motion  and  the  pair  of  roots  -2.1043  +  j  3.7184 
correspond  to  the   short  motion. 

The  fourth  degree  polynomial  can  be  factored  in  two 
second  degree  polynomials.  Each  one  of  the  second  degree 
polynomials   is   associated  with   the   corresponding  mode. 

Further  it  is  more  convenient  to  put  each  one  of  the 
polynomials   in   the   familiar   form 

52  4.2?0Or,S  ^-0Ov,'  (4-16) 

so  that  the  natural  frequency  cOti  f  the  damping  5  are  figured 
out  by  inspection  while  the  damped  frequency  can  be 
determined  by: 

oiy  =:  oor,  (^-5^)  ^  (4-17) 

If  we  denote  with  the  subscript  (SP)  the  short  period 
quantities  and  with  the  subscript  (PH)  the  phugoid 
quantities,  the  characteristic  polynomial  may  be  written  in 
the   form: 

(s'.25^(0„„s.co4)(s^.25„u)„„st<.)  -0  (4-18) 

Those  quantities  and  their  relation  is  shown  in  the 
complex  plane   in   Figure   4-3. 

The  roots  of  the  characteristic  polynomial  are  clearly 
shown  in   Figure  4-3    to   be 


-5^,a)T,p,    ±  joJd,,  (4-19) 

-5,pa;r,,p    ±  ju;^^  (4-20) 
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or 


fU 


(4-21) 
(4-22) 


i^I- 


u»<i^-u»va%)'*'" 


uw^-uM^a-sJ,)''*^ 


Figure  4-3 

Longitudinal   Roots   and  Characteristic  Quantities 
in   the  Complex  Plane 


Other    quantities    of    importance    in    the    longitudinal 

dynamics    are:      The    phugoid  period  denoted  by  T.^^  ,    the   short 

period  period  T^p     given  by: 

Tp«   =  2n/cAJd^^         (sec)  (4-23) 

Tsp  =  2x\/uj^^^         CsQc)  (4-24) 

and  the  time  for  half  amplitude  for  phugoid  and  short  period 

given  by: 
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^'''''     ^    5wa)r,p.  (4-25) 


XV2 


0.6^  Si 
SP     '     5spU3nsp  (4-26) 


The    Equations   4-25    and    4-26    come   from   the   fact   that  the 

oscillatory    response    decays    exponentially  described  by  the 

-Ket 
envelope       e  where     Re     is     the    real    part    of     the 

characteristic    root.       Therefore,    to    obtain    half    of    the 

amplitude   the   relation 

e  =  0.5  (4-27) 

must   exist.      By  taking  the  natural    logarithm  of   both   sides  we 

obtain  the   above   relations. 

For    our    example,    the    characteristic    polynomial    of 

Equation  4-13    can  be   factored   in  the   form 

(s^   +   0.0090    s   +   0.00395) (s^  +   4.2086s   +   18.2546)  (4-28) 

From  this,   by   inspection  we  can   determine: 
Wnp^  =    (0.00395)'^^       =   0.0628   rad/sec  (4-29) 

"Jnsp  =    (18.2546)''^        =   4.2725    rad/sec  (4-30) 


'?*! 


0.009)0 
2*0.0625 

4.2086 


-     a07l7  (4_3i) 


--     -     -^7^25       -     '-'^'^  (4-32) 

and  also 

^'^2,u     =     0.0717  '  0,OG2S       '    ^^^^^     ^^  (4-33) 
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,  O.G95i 

'^sp  0.4925  *  A.ni'b  (4-34) 

(%p„    =  0. 0628(1-0. 0717^)^^^  =   0.0627    rad/sec               (4-35) 

C4<sp    =  4.2725(1-0.4925^)*^^  =  3.7184    rad/sec                (4-36) 


2n 
0.0627      ~    '""•"'     "■'"  (4-37) 


TpH     =     .,?^,     =    >I(30.2I     5ec 


(4-38) 
It   is  worthwhile   to  note   that: 

(a)  The  damped  frequency  is  indicated  by  the 
imaginary  part   of   the   corresponding   characteristic   root, 

(b)  The  real  part  of  the  characteristic  root 
indicates  the  exponentially  decaying  envelope  of  the 
corresponding   oscillatory  motion,    i.e.    the   damping. 

(c)  The   following   relations  hold 

^  =  2Sv«aJnp^  +25spCDt>sp  (4-39) 

C    =   CUripM   -t- ^  Sp„  (JUrjpH  5sp  cOt,sp    +0)^5^  (4-40) 

1)   =   25^,4,  a)Tip«  cu^sp   +  2Sspct;Ti^p  Wnp„  (4-41) 

E   =   uu^p^   uu^gp  (4-42) 

The  two  types  of  dynamical  longitudinal  oscillations  are 

shown  in  Figure  4-4. 

The  forward  velocity  \J  and  the  pitch  angle  ©  are  usually 
excited  in  the  phugoid  mode  while  the  w  velocity  and,  hence, 
the  angle  of  attack  ot  remains  nearly  unchanged. 
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Figure   4-4 
Dynamic  Longitudinal   Oscillations 

In  contrast,  the  angle  of  attack  is  heavily  excited  in 
the  short  period  mode  as  well  as  the  pitch  angle  9,  while  the 
forward  velocity   remains  nearly  unchanged   in   that  mode. 

The  pilot  can  generally  control  the  phugoid  oscillation 
although  because  of  his  slow  reaction  it  is  possible  to 
worsen  to  unstable  due  to  out  of  phase  action  of  his  control. 
Phugoid  oscillation  must  be  highly  considered  in  take-offs 
and  landings.  Short  period  is  usually  out  of  pilot  control, 
but   is   of    importance    in   automatic   control. 

The  requirements  that  an  aircraft  should  satisfy  are 
written  in  appropriate  manuals.  It  is  usually  imposed  that 
the  oscillation  should  not  take  more  than  one  cycle  to  damp 
to  half  amplitude. 
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Flight  test  can  be  carried  out  to  determine  the 
oscillations  characteristics  under  different  conditions  of 
flight.  Gyroscopic  instruments  can  be  placed  and  recordings 
of  the  longitudinal  variables  can  be  made  possible  to  obtain. 
The  aircraft  is  trimmed  to  the  selected  speed  and  height  and 
longitudinal  oscillation  is  initiated  by  an  abrupt  elevator 
deflection* 

4.3      PHUGOID   AND   SHORT    PERIOD   APPROXIMATIONS 

The  fact  that  the  angle  of  attack  and  hence,  the 
w-velocity  remains  nearly  unchanged  in  the  phugoid  mode,  gave 
the   idea  of   determining  the  phugoid  quadratic  by: 

(a)  eliminating  the  w  terms. 

(b)  eliminating   the  pitching  moment  equation  since   no 
considerable  pitching  moment  variation  occurs,    and 

(c)  letting  q=s9,    obtaining  thus   the   determinant: 

-Zi>  -Uos  (4-43) 

The   characteristic  phugoid  quadratic   comes   out   to  be: 

For   example,    given,    the   characteristic  quadratic  is 

3^0. 0097s   +   0.0047    =    0  (4-45) 
which   has    roots 

-0.0049  ±  j    0.0681  (4-46) 
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These    roots   are  approximately  equal   to  the   exact   phugoid 
roots. 

From  Equation  4-44  it  is  clear  that  the  natural  frequency 
of  the  phugoid  is  heavily  dependent  upon  the  Zo  derivative 
while  the  damping  of  the  phugoid  is  dependent  upon  the  ^u 
derivative.  It  is  also  obvious  that  high  speeds  tend  to 
reduce  the  phugoid  frequency  and  to  increase  the 
corresponding  period. 

Since     for     subsonic    straight    level    flight    the   Zu 
derivative   is   simply  expressed  as: 

Zy  =  -2^  I  Do  (4-47) 

the   following  approximate   relations  hold  for   the  phugoid 

^**  Vo  (4-48) 

f 


'^*^     "    '  a\ll  ^  (4-49) 

Aircraft  designers,  in  order  to  achieve  larger  values  of 
phugoid  damping,  must  increase  the  value  of  the  ^u  derivative 
but  this  derivative  is  proportional  to  the  drag. 
Fortunately,  in  take-offs  and  landings  where  phugoid  motion 
is  undesired,  high  drag  devices  such  as  flaps  are  used 
anyway . 

A  more  precise  approximation  for  the  phugoid  quadratic  is 
given  by  leaving  the  w  terms  and  the  Mu  term  which  in  modern 
crafts  is  never  zero.  This  is  a  three  degree  of  freedom 
phugoid  approximation. 
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-Zu  S-Zva         -Vo  (4-50) 

-  Mo  -  Mw  0 

The   characteristic  phugoid  quadratic  becomes   then 


Mu  CfoXw  -j)  1         ^ 


My 


which  gives 


(4-51) 


and 


3    z'  _        Mu    _    X  T  ^/z 


^"^«-=[^^z„-^z«)] 


5    - 

-    ^             Mu  (UoXw-a)   - 

2f  ^  fz  -  '^^   7  V 

V2 

(4-52) 


(4-53) 


Short  period  approximations  can   be  made   possible  by: 

(a)  eliminating    the    v-equations    since    very    small 
v-velocity  change   occurs. 

(b)  setting    all   X-force    derivatives    to    zero    since 
very   small   X-force   variation  occurs, 

(c)  letting   q=s9,    obtaining  thus   the   determinant: 

s-Zw  -Do 

-sMw-Mvu  s-M(^  (4-54) 

The    characteristic    short    period    quadratic  comes   out   to 


be 


s^  +  (-UoMw  -Mfl  -Zw)s  +  (-aMw  +A>loZw)  =  0  (4-55) 
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For    the   example  given,    the   characteristic  quadratic   is: 
s^   +   4.208s   +   18.2556    =   0  (4-56) 

which  has   roots 

-2.104  ±  j    3.7187  (4-57) 

These    roots    are   approximately   equal    to  the   exact   short 
period   roots. 

In    summary,    for    typical    flight    conditions,    the   phugoid 

mode    can    be    considered    to    consist    of    changes    in    the 

u-velocity  and   the  0-pitch   angle,   while   the   short  period  mode 

can    be    considered    to   consist    of    changes    in    the   w-velocity 

(angle   of   attack)    and   the   pitch   rate  q. 

For    the    short    period   mode,    the    following    approximate 
relations  hold: 

OOtIsp    =    CM(|^Zw  ~VoM^]*'^  (4-58) 

'     "      2CM<^Zw  -UoM«,y'2  (4-59) 

A  critical   condition  occurs  when  the   center  of   gravity   is 
in   that  position  for   which 

MqZw    -  UoMv»/  (4-60) 

For   this   case    the   characteristic  equation  becomes: 

s  [  s^  (-(^oMw  -M(|^-Z,,)J  =0  (4-61) 

which   has    roots 

0     ;  OoMw  +  M(^  ^-Zw  (4-62) 

The    non-zero    root    is    usually  negative  which  yields    in   a 
non-oscillatory  and   heavily  damped  exponential    response. 
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4,4       SOLUTION  TO  THE    LONGITUDINAL    EQUATIONS 

The  longitudinal  equations  form  a  set  of  linear  first 
order  differential  equations  with  constant  coefficients. 
They  are  accompanied  by  a  set  of  equal  numbered  initial 
conditionsr  i.e.  initial  deviations  from  the  equilibrium 
flight   conditions,    and/or   a  number  of   control    inputs. 

Many  ways  have  been  developed  for  solving  sets  of 
differential  equations.  Laplace  Transform  method  is  usually 
the  simplest  although  it  is  not  suitable  for  computer 
solutions. 

In  the  development  of  the  solution,  we  would  like  to 
rearrange   the   equations   in   the  matrix  form: 

Hi)    =   A   ylt)    ^.  B  ull)  (4-63) 

where  i(t)  is  a  n-dimensional  time  varying  column  vector 
named  state  vector,  ii(t)  is  a  r-dimensional  time  varying 
column  vector  named  control  input  vector,  ^  is  a  nxn 
constant  matrix  and  £  is   a  nxr  constant  matrix. 

Taking  the  Laplace  Transform  of  Equation  4-63  according 
to  theorem  3   of   Appendix  A,   we   obtain: 

sX(s)    ^  /(oj  +  AXCs)  4-B  UCs)  (4-64) 


or 


or 


(si  -A)  X(s)    =  x(o)  +  -l  U(5)  (4-65) 


X(s)  ^  (sI-Ap  x(o)  ^  (sI-A)"  ?l((s)  (4-66) 

where    X    is    the    nxn    identity    matrix    and    i(0)     is    the 
initial   condition  column  vector. 
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The  first  term  of  the  transformed  solution  corresponds  to 
the  initial  conditions  and  is  the .homogeneous  solution  and 
the  second  term  corresponds  to  the  control  input  and  is  the 
particular  solution. 

By  inverse  Laplace  Transforming  Equation  4-66  we  obtain 
the  time  solution: 

xCl)  =[i"'CsI-Af  ]  xCol  ^  /^[(sI-AT'?U(5)]     (4-67) 

To  implement  this  solution,  in  a  digital  computer,  time 
domain   formulation   is   followed. 

The  key  is  to  multiply  the  matrix  differential 
Equation  4-63   by  exp(-At)    obtaining  thus: 


-^^  -  1  All 


(4-68) 


or 


Integrating  between  0  and  t.  Equation  4-69  gives 


Q      x(y   =  xCo)  f  J    e       8  i^CT)dT 


(4-70) 


or 


x(^:)  =  e'     xCo)  f 


•^     Aa-r) 


?  uCr)di 


Jo  (4-71) 

The   exponential   term   is   expanded   in   Taylor    series   as: 

(4-72) 

In  reality,  a  series  truncated  to  ten  or  twenty  terms  is 
often  quite  sufficient  for  normal  accuracy  requirements. 


e       =1  +At  ^^A-t  ^  ■■■■  "r  -^  A  i 
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The  computer  program  determines  the  values  of  the  state 
variables   at  time ^intervals   selected  by  the   user. 

We  now  go  back  to  the  unt ransf ormed  longitudinal 
equations  and  try  to  arrange  them  in  the  form  of  the  matrix 
Equation   4-63 . 

For  convenience,  the  unt ransf ormed  equations  are 
replicated   below: 

v  -  (Xu)u  -(Xw)w   fae    =(X5)S  (4-73) 

-{Zo)v    f  W   -(Z«,)w   -CO.)^   =  (Z5)5  (4-74) 

-(Mu)u    -(Mw)w    -(Mw)w;    +i-{^%)(l    ^  (^^)h  (4-75) 

Q=    %  (4-76) 

Equation  4-73    is  written: 

u  r  (Xu)u  fCXv«)w  -(^)e4Cys)S  (4-77) 

Equation  4-74    is  written: 

«0   r   (2v)v   +(Zw)w    irCOo)c^  i.Cz^)S  (4-78) 

Equation  4-75    is  written: 

^  r  CMOu  ^-CMw)w  4(Mw)w  +(M<^)^+CM5)5  (4-79) 

or    substituting   Equation  4-78   for  i/   we   obtain: 

^  =CMui-MvOZu)L>   1-CMvtf  f  MwZw)a/  +CMwUo)<^  fCM&+MwZ5)5        (4-80) 
or    if  we   define: 

Mif  =    My  A-UijZo  (4-81) 


-* 


Mw    =    Mw  V  Mw  Zw  (4-82) 

M^    =    ^q^  M«U,  (4-83) 

Mi*  -     Ms  ^MwZ5  (4-84) 
we  will   have: 

(5    =  CM;T)u   ^C^t>)w  +CM|)<^  ^Cm1)§  (4-85) 
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It  is  further  desirable  to  include  the  height  of  the 
vehicle  as  another  variable  of  interest.  The  equation  for 
the  height  is  obtained  from  the  original  Z-equation 
Equation  2-134   which   for   ^o  =0   becomes 

i-  Z    =  v^  -l^oq  (4-86) 

The  left-hand  side  of  this  equation  is  just  the  z-axis 
acceleration,    i.e. 

-a^  ^-ii*    .  W-U,^  (4-87) 

Integrating  the  last  equation  with  respect  to  time  we 
obtain 

h   =  -  \M  +  1^9  (4-88) 

the  constant  of  integration  being  the  initial  w-velocity 
perturbation  which  is  zero.  Equation  4-88  describes  the 
height  perturbation  around   the  nominal    value. 

Collecting  Equations  4-77,  4-78,  4-85  and  4-88  in  matrix 
form  we  obtain  the  so  called  state  variable  representation  of 
the  longitudinal  dynamics  model,  the  state  variables  being, 
u,w,q,9,    and  h. 


V 

% 
9 
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O 
O 


Zvu 


0 

-1 
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<,       M'. 
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0 

0 
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0 
0 
0 
0 
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xs 

w/ 

Z5 

•  % 

4- 

Nf5 
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0 

k 

<3 

(4-89) 
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state  variable  representation  form  is  perfectly  suitable 
for  computer  solution.  Many  ordinary  differential  equation 
solver  software  routines  are  available  for  solution.  Most  of 
these  are  provided  with  plotting  capabilities  of  the  state 
variables. 

The  solution  that  follows  corresponds  to  the  given 
example  with  initial  conditions  and  control  inputs  indicated 
in  the   summary  of   specifications  tables. 

It  can  be  seen  from  the  plots  that  the  v-velocity,  the 
©-pitch  angle  and  the  h-height  are  heavily  excited  in  the 
phugoid  motion,  while  the  w-velocity  and  the  q-pitch  rate  are 
heavily  excited   in   the   short  period  motion. 

For  the  phugoid  approximation  the  following  model  may  be 
used: 
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-ZuM 

0 


0 


0 

V 

yh 

0 

Q 

+ 

Ih/Vo 

o 

k 

0 

(4-90) 


while   for    the   short   period  approximation; 


\w 


% 


Zv 


K       m{ 


vO 


% 


Z5 

Nfa 


(4-91) 


4.5       LONGITUDINAL   AERODYNAMIC    TRANSFER   FUNCTIONS 

Up  to  now  we  went  through  the  analysis  of  the  dynamical 
behavior  of  the  flight  vehicle,  by  making  use  of  the 
mathematical  model  that  adequately  characterizes  the 
longitudinal   characteristic  motions.      In  any  dynamical   system 
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two  types  of  variables  exist,  the  input  and  the  output 
variables.  The  input  variables  influence  the  output 
variables  in  a  manner  determined  by  the  dynamics  of  the 
system. 

The  objective  of  this  section  is  that  of  determining  the 
nature  of  those  input  -  output  relations  called  transfer 
functions. 

For  the  longitudinal  model,  the  input  variable  was  the 
elevator  control  deflection  §  and  perhaps  the  initial 
conditions  or  disturbances  that  can  be  treated  as  inputs, 
while   the   output   variables  were  u,w,q, 0  and  h. 

The  corresponding  transfer  functions  are  obtained  by 
solving  simultaneously  at  a  time,  the  transformed  equations 
of  motion  for  the  output  variable  of  interest,  keeping  all 
initial   conditions  zero. 

The   transfer  function  ©(s)/5(s)    is 


9(s) 


-Mw       -sMw-Mw       s.M 


% 


I 


V§ 
MS 

0 


As^4.'Ss^^Cs^+3s+6 


(4-92) 


where 


A^    =    Ms  i-U^Ls  (4-93) 

'B©    =    X5CZ^Mw+Mu)^ZaCM,.-X„Mv;,)  -MiCXy  +  Zw)  (4-94) 

Ce   ^     XS(Zi;Mv^-Zu/Mu)+Z5CMuXw-M„)(w)f  MS(Z»/Xt;-Xv«2u}         (4-95) 
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or 

S(5)     '       ;^s^4.Bs5  4.  Cs^  +l?i^■E  (4-96) 

where     ^,  ,   Z©^       ^^^®    ^^^^    roots    of    the    numerator    polynomial 
called   zeroes   of   the   transfer   function. 
Similarily  we  can   derive 

y:^)       sBCs)  Aesrs+Ze.)(s+z©,) 


^(s>  5(5)      '       A.^.^.B^-Cs^+J^  +  E  (4-97) 


5(5) 

As^^-^s^Cs^+i^  +  E 

wCs) 

Aw  S^^^w  S^  +  CwS  +Dw 

SCs)  As^^Bs'  +  Cs'+3>s+e  (4-98) 


where 


or 


Aw  :-Z§  (4-99) 

'Bw  =  )<5Zy   -Zi()(v^M^)  ^  UhVo  (4-100) 

Ca)   =  X§  CUoMu -ZuM(^)^.  ZSXu/Ufl   -OoUilXv  (4-101) 

3w   =    oCZsWy  -M5Zu)  (4-102) 

This   transfer   function  may   take    the   following   forms 

wCO  Aw  Cs+zw.)  Cs^+-25w"^»/s+ oj^) 

IC^    "  A?7i?Tc?^iITi  (4-103) 


^^;(s)  Aw  Cs-(-2w,)  C  S+Z^j)  fs+Zu;,) 


5Cs)  As^  +5s^^Cs^4.Ps  +  E  (4-104) 

We    can    therefore   have    three  zeroes,    two  complex  and  one 
real   or    three   real. 

Other   longitudinal   transfer   functions  are: 

IcsT    "    "a^^T^vTCs^T^TT?  (4-105) 
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where 

Au  =  y^  (4-106) 

^u  =  -XS  CZyyf  Ma  +i;<,Mw)  +ZsXw  (4-107) 

Cu  =   XS  (ZwM^j^-UoMw)  -  Za(yvjlM<j^f|jK^)  ,iUSfUoXv.-_gJ  (4-108) 

1>y  =  Q  CMsZw  -2s  Mw)  (4-109) 
This   transfer   function  may  also   take    the   forms 

oCs)  Ao  Cs+Zu,)Cs^^2^uaJuS+aJ„^) 


or 


SCs)    "  /s^  +  Bs^Cs^+DsfE  (4-110) 


UCs)  Au    Cs+Zi/jCs+ZyjCs+Zyj) 


5(s)  As'^^BsUCs^+Tls  +  E  (4-111) 

Similarily 

S(sJ    '     s(-4s'^+'Bs'<-Cs^^IlsvE)  (4-112) 
where 

\r-Z&  (4-113) 

^H   =  -  VaZy    -HZ6    CM(^+a,Mii;  ^Yo]  (4-114) 

Ck  =    XSZo  CM^^a,/U;,)-26[yw()Ui{^^U,Mw}-aiWw]  -  MSDo^vu            (4-115) 

3k  =  -^s^o(ZwMw  -MwZy)  +Z5[iMi/CUoXvs/-|)  -y<,MwXu] 

+  MS  [UoZwXu  -  Zu  Cu.yw-^)J  (4-116) 

and  may  take  for  forms 

h(s)     Ak(s^-Zh,)Cs'^2Sh(UnS  fcu^f  ) 


or 


K(s)   _   4  Cs^^Z^,.)(s^^ZHJ^S-^2w^) 

S(s)  "    3CAs^^Bs^Cs2^-a?s4.£)  (4-118) 
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Another    transfer    function   of    interest   is   that   relating 
the   z-axis   acceleration  to  the   control    input,    i.e. 

^^^  "      5rs)    ~  AlVB^TcvT^rri  (4-119) 

or 

azfeJ    _      ylh  s  Csf2^„)CstZh,)Cs4.^^,,■) 

(4-120) 
Those  transfer  functions  are  used  in  automatic  or  manual 
flight  control  systems,  since  they  specify  the  input  -  output 
relationships  of  the  variables.  For  example,  to  find  the 
open  loop  response  ©(t)  due  to  a  unit  step  elevator  input  we 
use   the   transfer   function  Equation  4-96 


^s'-nss'+^^s'-t  j/'sf  c  (4—121) 

Taking    the    inverse    Laplace    Transform   of    Equation   4-121   we 
obtain   the   time   response. 

4.6      LONGITUDINAL    EQUATIONS    IN    NON-DIMENSIONAL    SYSTEMS 

The    equations    of    motion   of    a    flight    vehicle,    can    be 

written    in    a    number    of    different    forms,    depending    on    the 

particular    axes    system    selected    for    the    definition   of 

variables   and   the   stability  analysis. 

We    have    just    described    the    body    axes    system    in    the 

dimensional    form,    in  which  roost   stability  work   is   carried   out 

[Ref .    1]  . 
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In  most  aerodynamic  problems  it  is  more  convenient  to  use 
non-dimensional  coefficients  to  represent  the  parameters 
involved.  By  this  non-dimensional  system,  the  equations  of 
motion  can  be  described  in  the  body  axes  or  wind  axes  system 
as  well. 

The  non-dimensional  coefficients  have  the  advantage  that 
the  effect  of  speed,  size  and  air  density  are  automatically 
accounted.  However,  it  is  required  to  express  the  mass  and 
the  time  in  non-dimensional  forms  as  well.  Therefore,  one 
cannot   obtain   real   time   solutions. 

In  the  non-dimensional  body  axes  system  the  longitudinal 
equations  of   motion  are  given  as    reference    [Ref.   2]. 


2us  -Cx 


-C; 


^ 


2Q^-Q,         2fjs.sCz- -Cz^         -3ps-sC2^ 
-Cwy  -sCm-  -Cm^  ij.  s^  -  sCm 


I 


V 
ot 
9 


0 
Czn 

s Cm-  f  Cwr, 


Yi  (4-122) 


where    the    non-dimensional    mass,     inertia    and    time    are 
expressed   as: 

j/  r  ^  /pSe  (4-123) 

2   being    the    characteristic    length    taken   as   the   half  of    the 
aerodynamic   chord 


U 


I 


ya 


(4-124) 
(4-125) 
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V 

0 

u 

■= 

Czr. 

0 

sCw  •  tCrw 

yi      (4-126) 


In    the    non-dimensional   wind   axes    system   the   longitudinal 
equations  are    [Ref.    3] 

where    the    non-dimensional    mass,     inertia    and    time    are 
expressed  as: 

p  =  w/pS£  (4-127) 


-y 


=  2 


T  =  om  /p  SDo 


(4-128) 
(4-129) 
(4-130) 
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TABLE  4-1 

SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE 
LONGITUDINAL  EQUATIONS  DEDICATED  FOR  PHUGOID 
MOTION  DUE  TO  INITIAL  CONDITIONS 


A.  Variables  and  Initial  Conditions 
u  =  5.0  ft/sec 

w  =  2.5  ft/sec 
q  =  0.05  rad/sec 
d  =   0.075  rad 
h  =  10.0  ft 
t  =  0.0 

B.  Stability  Derivatives  and  Constants 
Xu  =  -0.0097    Zo    =   -0.0955 
Xw  =  0.0016     Zw  =  -1.43 
Xa  =  0.0       Zs    =  -69.8 


C.   Special  Functions 


Mo    =    0.0 
Mw   =    -0.0235 
Mw   =    -0.0013 
M<i    =    -1.92 
Mi    =   -26.10 

g   =   32.174 
Uo    =   66  0.0 

D.      Control 
§  =   0.0 

Input 

M  *  ^ 

Mu    +  Mw  Zy 

m:  . 

Mw  ■«•  Mw  Zw 

"i  = 

A^^^  Mw(;o 

M8*  = 

Ms  *  UwZ5 

E.  Derivatives 

q-    My  u+ Ww  w  +-<M^(^  <-/^|  5 

F.  Outputs 

u,    w,    q,  9  f    h  vs   time   at   interval    0.0625 
end   calculation  when   time  > 250.0 
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Figure   4-5 
u-Velocity  vs  Time   for   Specifications  of   Table   4-1 
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q-Pitch  Rate  vs  Time  for  Specifications  of  Table  4-1 


183 


j -j 1 1 

; -1 \ j 

\- 1 1 1 

-. 

; :.._ i : 

1 

1 

in 

o 

; ; 

■ — k 







1 1 \ \ 

-                                                           : 

1 f- \ \ 

i        ;        ^       i 

\ 1 ; 1 

MM 

1   M   i 

o 
o 

)                        ■ 

;             :              :              : 

;            •■           i            : 

M-fl 

\--Ar---\ -i 

..'..'..         '              ■ 

\ 1 i- \ ■■ 

If) 

o 
o 

V  .     .     ; 

"        ^        ^ 

:Mt 

-j \ i i/-- 

^-hi ; 

•i i 1 j 

MM 

1    i    i  j/' 

"iX  1 

m\  M 

1    pT  M 

M 

!       \     1 

o 

□ 
o 

\ / ' \ : 7 

\: 

:OGi              ;           I 

DOS;              y^                  010;         ;        \         ; 

D15; 

/ 

]20i 

\  ' 

MA; 

\                      '               M       i\l 

; 

/ 

i  \^ 

MlM: 

;        /                           M       1      \ 

"S 

in 

■  r 

M    M\j 

o 
o 

1 

j 

;         :         ■         :         1         . 

1 

THETR    VS    TIME 


X-SCflLE  =  5M0E  +  01    UNITS    INCH, 
T-SCflLE  =  5M0E-02    UNITS     INCH. 

Figure   4-8 
9-Pitch  Angle   vs   Time    for    Specifications   of   Table   4-1 
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TABLE    4-2 

SPECIFICATIONS   FOR  COMPUTER   SOLUTION   OF   THE 
LONGITUDINAL    EQUATIONS   DEDICATED   FOR   SHORT    PERIOD 
MOTION   DUE    TO    INITIAL   CONDITIONS 


A.  Variables   and   Initial   Conditions 
V  =   5.0 

w  =   2.5 
q  =    0.05 
6  =    0.075 
h   =   10.0 
t   =   0.0 

B,  Stability  Derivatives   and  Constants 


Xy     =   -0.0097 
Xw   =   0.0016 
X5    =0.0 


Zu    =   -0.0955 
Zv^  =   -1.43 
Zi    =   -69.8 


C.      Special   Functions 
My    =     Mv  +  MwZw 

Ml    =    Ms  <-  Mw  Zs 
E.      Derivatives 


Mu    =    0.0 
Mw   =   -0.02355 
Mw  =   -0.0013 
n,^   =   -1.92 
Ms    =   -26.10 


g   =  32.174 
v^=   66  0.0 


D.      Control    Input 
S  =   0.0 


^1  =  -uj  +  ue 


0    =- 


i 


F.      Outputs 

V,   Wr    q, 0  ,    h  vs   time   at   interval    0.00125 
end   calculation  when   time  >5.0 
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U-Velocity   vs   Time   for   Specifications  of   Table   4-2 
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w-Velocity  vs  Time    for    Specifications  of   Table   4-2 
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q-Pitch   Rate   vs  Time    for    Specifications  of   Table   4-2 
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9 -Pitch  Angle  vs  Time  for  Specifications  of  Table  4-2 
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TABLE  4-3 

SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE 

LONGITUDINAL  EQUATIONS  DUE  TO  INITIAL 

CONDITIONS   AND   CONTROL    INPUTS 


A.   Variables  and  Initial  Conditions 
V  =   5.0 
w  =  2.5 
q  =  0.05 
9   =  0.07  5 
h  =  10.0 
t  =  0.0 


B.   Stability  Derivatives  and  Constants 

Xu  =  -0.0097    Zo    =  -0.0955    Mu  =  0.0 

Zw  =  -1.43 
Zj  =  -69.8 


X*,  =  0.0016 
Xs    =0.0 


C.   Special  Functions 

Mu  -  Wi,  +  MwZu 
l^w  =  Mw  +  Uw  Zw 

E,   Derivatives 


Mw    =   -0.0235 
M;,   =   -0.0013 


Ma    =   -1.92 
Ms    =   -26.10 


g   =   32.174 
Vo     =66  0.0 


D.      Control    Input 

S  =   0.01    for    0<t<10.0 

S-   -0.01  for  10.0<t<20.0 

5  =  0.0  for  t>20.0 
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F.   Outputs 

V,   w,  qr  9  f  h  vs  time  at  interval  0.0625 
end  calculation  when  time  > 250.0 
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U    V5    TIME 

X-SChL'E  =  5„  OGE  +  Ol  UNITS  INCH, 
T-SCHLE-2. uGE+01  UNITS  INCH. 

Figure   4-15 

v-Velocity  vs  Time   for   Specifications  of   Table   4-3 
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X-3CflLF-5. GOE+Ol  UNITS  INCH= 
T-3ChL[E  =  5=  UGE  +  00  UNITS  INCH. 

Figure   4-16 

w-Velocity  vs   Time    for    Specifications  of   Table   4-3 
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Q    VS    TIME 


X-3CRLF=5» QGE+Ol  UNITS  INCH. 
T-5CRLE-E. UGE-G2  UNITS  INCH. 

Figure  4-17 
q-Pitch  Rate  vs  Time  for  Specifications  of  Table  4-3 
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HETn    VS    TIME 


X-DLflLE=5. OOE+Ol  UNITS  INCH,  ■ 
T-5CHLF=5» GCE-02  UNITS  INCH. 

Figure   4-18 

©-Pitch  Angle  vs  Time   for    Specifications  of   Table   4-3 
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Figure   4-19 
h-Height   vs   Time   for    Specifications  of   Table   4-3 
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CHAPTER  5 
LATERAL  DIRECTIONAL  DYNAMICS 

5.1   INTRODUCTION 

In  this  chapter  we  will  analyze  the  lateral  directional 
motion  of  a  conventional  airplance,  by  folliwng  the  same 
treatment  as  in  the  preceding  chapter  on  longitudinal 
dynamics. 

In  the  lateral  directional  case  also,  any  particular 
response  will  be  associated  as  in  the  longitudinal  case,  by 
its  corresponding  characteristic  root  location  in  the  complex 
plane. 

Further  simplified  set  of  equations  will  be  developed 
that  apply  to  the  individual  modes  of  motion.  Finally, 
computer  solutions  will  be  emphasized  as  well. 

As  an  example,  the  lateral -directional  motion  of  the  same 
aircraft  as  in  the  preceding  chapter  will  be  analyzed  with 
the  lateral-directional  parameters  indicated  below: 
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dtlilude     C-fl)  20.000 

w/eight         Clb)  30.500 

omacb  01  umber  O.G58 

«it"  speed   C(+Acc)  G60 


Yv  = 

-  0.0829 

L^  - 

-4.77 

W^    =    3.55 

Y5  = 

^.656 

Lp   - 
^5    - 

'  -I.&95 
0.4776 
Q7.25 
O.G66 

Np   =    -0.0025 
Nr    =    -0.09c,7 
M5  -    -OGIS 
N^  =    --I.58S 

^«/rxy   = 

:    0.0665 

I«/l22 

=  0.0370 

5.2      CHARACTERISTIC    POLYNOMIAL   OF    LATERAL-DIRECTIONAL 
MOTION 

Lateral-directional  equations  of  motion  are  further 
simplified  by  neglecting  the  effect  of  Y^ f  Yp ,  Yr  r  ^v r  and  Ny 
derivatives  which  are  of  no  importance  in  lateral  directional 
dynamics.  Furthermorer  a  straight  level  flight  is  assumed  as 
the  reference  flight  condition  (i.e.  ^„  =0)  .  Hence, 
Equation  2-153    is   reduced   to: 


s-Yv 


0 

0 


0 

s-Lp 

1 

o 


1 

•s U 

Jxx 

S-Wr 

O 

< 


0 

o 

-s 

o 


o 

O 

o 
o 

-S 


^ 

Ys* 

p 

Ll 

r 

•= 

Ni§ 

4> 

0 

1^ 

0 

(5-1) 
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The    characteristic   polynomial,    as    known,    is    found    by 
evaluation  of   the   above  determinant  which  comes   out   to  be: 

As-^^-  Bs'  +  Cs^  -^Ds  -HE    =0  (5-2) 

where 

I'- 

A  =  ^  - 


IxxTzz  (5-3) 

B-=  -  Yv  0-  ~r^—]  -i-P  -  Nr  -  4^  Np  -  ^  U 

*■  Ixx  Izz  I  Ixx       ^        Izz     "^ 


(5-4) 


(5-5) 

-^  =  -N^Lp.Yv(VjpLr-LpWr).NpL^-^(L^-^N^)  (5_g^ 

E=   ^^    (l^Mr-N^Lr)  (5-7) 

Equation  5-2  is  referred  to  as  the  lateral-directional 
stability  quartic.  The  roots  in  this  case  for  any 
conventional  vehicle  are  two  real  negatives  or  positives 
which  are  associated  with  two  convergent  or  divergent  motions 
and  one  pair  of  complex  roots  associated  with  an  oscillatory 
motion. 

The  overall  lateral  directional  motion  of  the  vehicle  is 
a   superposition  or    a  combination  of    all  three  modes. 

The  complex  pair  of  roots  correspond  to  the  as  named 
Dutch  Roll  oscillation,  while  the  large  real  root  corresponds 
to  the  Roll  Subsidence  mode  and  the  small  root  to  the  Spiral 
mode. 
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For  the  numerical  example  given,  we  find: 

A  =   0.9975  (5-8) 

B  =   1.867  (5-9) 

C   =   3.6841  (5-10) 

D   =   6.2637  (5-11) 

E   =   -0.0086  (5-12) 

and   the   characteristic  polynomial    is: 

0.<3976  s'^    +  lSC7s^+  S.GS4I  s^  4  6.2657  s-  C3.CX3&S   =0  (5-13) 

which  has    roots: 

-0.0465   +    jl.8784  (5-14) 

+0.014  (5-15) 

-1.7801  (5-16) 

The    pair    of    roots,    -0.0465  ±  jl.8784,    correspond   to  the 
Dutch    Roll   mode,    while    the    root   +0.0014   corresponds   to  the 
Spiral   mode   and   the   root  -1.7801   to   the  Roll   Subsidence   mode. 
Further,    it    is  more  convenient   to  put   the   characteristic 
polynomial    in   the   form: 

CsS25,,aJ.„,s+aJ,'^J(s-rJ(s-reO  (5-17) 

so  that  the  natural  frequency,  the  damping  and  the  real  roots 
are  obtained  by  inspection  while  the  damped  frequency  of  the 
Dutch   Roll   can  be  determined  by: 

'^d^  =  <"-.«(  (-5^..)"^  (5-18) 

The    root   plotting   in   the   complex  plane  will  appear   as    in 
Figure   5-1. 
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Figure  5-1 

Lateral -Directional  Roots  and  Characteristic 
Quantities  in  the  Complex  Plane 


The  imaginary  roots  of  the  characteristic  polynomial  can 
be  written: 

-5^Ru;n,,  ±  juJ^^  (5-19) 

or 

-So^COr,^  -  l^r,,^(^-Sl^y'^  (5-20) 

Other  quantities  of  importance  in  the  lateral-directional 
dynamics  are:   the  Dutch  Roll  period  denoted  by  T-dr  given  by: 

T„R=2n^,  (5-21) 

the  half  amplitude  for  the  Dutch  Roll  given  by: 


i 


V2 


J>B 


(5-22) 
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the   half  amplitude   for    the  Roll   Subsidence   mode    (if   stable) 

I  0.6°>51 

'^'«     "  ^«^                                                      (5-23) 

and  the  half  amplitude  for  the  Spiral  mode  (if  stable) 

,  0.G95J 


s  fs  (5-24) 

For  our  example,  the  characteristic  polynomial  of 
Equation  5-13    can  be   factored  in  the   form: 

Cs^+  0.0*55  5    +  S.5  505)  C  s- 0.0014)  Cs+ -I. T80!)  (5-25) 

From  this   by   inspection  we  can   determine: 

a)yijj^  =    C  3.5505)"^  r    iSl'^       rad/sec  (5-26) 

^^^    =    0.095/2.  I.8T9  =  0.0247  (5-27) 

■^/2dr   =     0.G95I/0.O247- -1.379    -   I4.<^51    sec  (5-28) 

-t'/2g5    =    0.6951 /jYsoi     "   0.3894  sec  (5-29) 

Many  of  the  basic  ideas  involving  longitudinal  stability, 
also  apply  to  lateral  and  directional  stability.  Lateral  and 
directional  stability  are  interrelated  due  to  cross  coupling 
effects  of  roll  and  yaw  motion,  that  give  rise  to  the  three 
modes   of   motion  discussed   so   far. 

In  the  Spiral  mode,  the  airframe  seems  directionally 
nearly  stable,  with  no  sideslip  but  it  is  excited  in  a  banked 
turn.  The  unstable  mode  is  found  usually  in  large  finned 
airplanes  with  no  dihedral,  so  that  the  side  forces  developed 
tend  to  turn  the  plane  into  the  relative  wind  but  since  the 
outer  wing  travels  faster,  generates  more  lift  and  the  plane 
will    roll    to    still    a    higher    bank    angle.      The   flight   of    the 
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airplane  is  a  heavy  banked  turn  of  ever-decreasing  radius, 
i.e.  a  tightening  spiral.  The  stable  mode  is  found  in 
airplanes  with  dihedral  where  the  flight  is  a  slight  banked 
turn  of  gradually  increasing  radius.  The  spiral  mode  is 
illustrated   in   Figure  5-2. 


\mhai    ^(jgM  po^h 


spiral   oonverjoici 


spinal  dinvr^nce 


Figure  5-2 
Spiral  Stable  and  Unstable  Modes 

In  the  Roll  Subsistence  mode  the  airframe  seems  laterally 
nearly  stable  with  no  bank  but  it  is  excited  in  a  sideslip 
turn.  The  unstable  mode  is  present  when  side  forces  generate 
yawing  moments  that  continue  to  increase  the  initial  sideslip 
disturbance.  This  condition  may  continue  until  the  airplane 
is  broadside  to  the  relative  wind.  Airplanes  with  successive 
weathercock  stability,  usually  have  this  mode  stable  where  an 
initial  sideslip  tends  to  die  out.  The  Roll  Subsistence  mode 
is  illustrated  in  Figure  5-3. 
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Figure  5-3 
Roll   Subsistence   Stable  and  Unstable  Modes 

Dutch  Roll  is  a  motion  which  has  the  characteristics  of 
the  Spiral  convergent  and  Roll  Subsistence  divergent  mode. 
It  is  found  in  airplanes  with  strong  lateral  stability  but 
weak  directional.  If  a  sideslip  disturbance  occurs  and  the 
plane  is  yawing  in  one  direction,  rolling  moments  are 
developed  in  a  counter  motion,  i.e.  yaw  to  the  left  is 
causing  roll  to  the  right  and  the  plane  wags  its  tail  from 
side   to   side.      Dutch  Roll  mode    is   illustrated   in   Figure  5-4. 

The  damped  frequency  of  the  Dutch  Roll  is  indicated  by 
the  imaginary  part  of  the  complex  root,  while  the  real  part 
indicates  the  damping.  To  reduce  the  effects  of  Dutch  Roll, 
an  increase  in  the  directional  and  a  decrease  in  the  lateral 
stability   is    required. 
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Figure  5-4 
Dutch   Roll  Mode 

5.3      DUTCH    ROLL,    SPIRAL,    AND   ROLL   SUBSISTENCE    MODE 
APPROXIMATIONS 

By  analogy  with  the  longitudinal  short  period  and  phugoid 
approximations,  we  can  reduce  the  degrees  of  freedom  to 
obtain  Dutch  Roll  Spiral  and  Roll  Subsistence  mode 
approximations. 

Dutch  Roll  is  a  motion  which  is  characterized  by 
negligible  variations  of  roll  and  sum  of  the  rolling  moments 
nearly  zero.  The  Dutch  Roll  quadratic  is  obtained  from  the 
matrix  Equation  5-1,  by:  (a)  eliminating  the  cp ,  and  p 
equations,  and  (b)  setting  all  rolling  moment  derivatives  to 
zero   obtaining   thus,    the   determinant: 


s-Yv 


i 

S-Nr 


(5-30) 


The  characteristic  Dutch  Roll  quadratic  comes  out  to  be: 

s^  -  (Yv+Mr)s  ^  (  YvMr +N,%)  =  O  (5-31) 
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For  the  example  given,  the  characteristic  Dutch  Roll 
quadratic   is 

s^4-O.I7S6s    +5.5579=0  (5-32) 

which  has   the   roots 

-0.0893   ±  jl.8841  (5-33) 

These   roots   approximate  the   exact   Dutch  Roll   roots. 

From  Equation  5-31  it  is  clear  that  the  frequency  of  the 
Dutch  Roll  is  heavily  dependent  on  the  weathercock  derivative 
^(^  and  the  damping  on  the  yaw  damping  derivative  Nr  and  the 
^^    derivative. 

All  the  derivatives  involved  in  the  Dutch  Roll  mode  are 
dependent  mainly  from  the  vertical  tail,  and  secondarily,  on 
the  wings  and  fuselage.  Thus,  proper  care  has  to  be  taken 
for  the  vertical  tail  as  far  as  the  Dutch  Roll 
specifications.  It  is  usually  required  that  Dutch  Roll 
oscillations  should  be  damped  to  half  amplitude  within  one 
cycle. 

A  first  approximation  to  the  Spiral  and  Roll  Subsistence 
mode  is  obtained  by  expanding  the  determinant.  Equation  5-30, 
by  the  roll  damping  term  (3-Lp)  or  s(s-Lp)4>.  The  determinant 
takes   the   form  of : 

s-Yv/  0  { 

0  sCs-Lp)         0  (5-34) 

-M|i  0  s-Hr 
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The  characteristic  equation  comes  out  to  be: 

S('5-Lp)[s'-  (Yv  +  Wrls  -t  CYvNr+N|i)]    =o  (5-35) 

Here,  the  free  5  corresponds  to  the  Spiral  mode,  the 
(s-Lp)  corresponds  to  the  Roll  Subsistence  mode,  and  the 
quadratic  to  the  Dutch  Roll  mode. 

By  this  approximation,  for  the  given  example,  the  Spiral 
root  is  zero,  the  Roll  Subsistence  root  is  -1.695  and  the 
Dutch  Roll  mode   as   previously  -0.0893   ±  jl.8841. 

A  more  precise  approximation  for  the  Spiral  and  Roll 
Subsistence  modes,    is   obtained   by  expanding   the   determinant 

0  -VVo  i 

-L|4  sCs-Lp)        -Lr 

-NjJ  s-Hp         S-Nr 

This   lead   to 


(5-36) 


9 


§ 


or 


M^s^  +  (HpL^-N^Lp  -f-L^)s  .  ^  CL^Nr-N^U]  =0 


(5-37) 


(^         "y       Vo    M(j  r      (Jo  ^    Nfi        "  >  (5-38) 

For    the    given    example    one    obtain    the    characteristic 
polynomial 

sH-l.7G59s    -0.0024    =0  (5-39) 

which  has    roots 

0.0014,    -1.7653  (5-40) 

These    roots    approximate    the    exact    Spiral    and    Roll 
Subsistence    roots. 
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5.4       SOLUTION   TO  THE    LATERAL-DIRECTIONAL    EQUATIONS 

Following  the  same  approach  as  in  the  preceeding  chapter, 
we  first  try  to  arrange  the  lateral -directional  equations  in 
the   form 

iii]   -  A  xCt)  t|  u[h)  (5-41) 

For  convenience,  the  unt ransf ormed  equations  are 
replicated  below: 

|g-(Yy)^    +r    -M0o)(^     =   W)^  (5-42) 

-if)^    +f  -CLp)p  -(Iyz/j,^)r   'Cir)r   ^LU)S  (5-43) 

-(W^)/i    -(Uz/U^)p   -(Np)p   +r  -CNr)r  rCM5)§  (5-44) 

Cp    =  f  (5-45) 

V    =    r  (5-46) 

Equation  5-42    is  written 

^    =  (Y.)(^  -r  f(9/a)4>   +CY5)5  (5-47) 

Equation  5-43    is  written 

9  =    C^)|3    +CLp)p   t(r«/i,,)r   ^Cir)ri-Ck)l  ^iis)^        (5-48) 
Equation  5-44    is  written 

f    =    (N^)^   fClxz/l^Jp    ^(Np)p  ^-CMr)r  +(:M5)^4<N5^5     (5-49) 
In    Equations  5-47,   5-48    and   5-49,    the   control   inputs   are 
expanded    in    terms  of   aileron    it)    and   rudder    (^)    deflections, 
as     indicated    by    Equation    3-143    neglecting    the    dotted 
derivatives   !-§  f  J^^    . 
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Substituting  Equation  5-49  into  Equation  5-48  we  obtain 


(5-50) 
Similarly  substituting  Equation  5-48  into  Equation  5-49 
we  obtain 


(5-51) 


We  now  define  for  convenience 


,  *     ^(i  ^  ^%x  M(i 

V  ""   ,  t2  / (5-52) 


etc.,    so   that  Equation  5-50    and   5-51   can  be  written 

?■   .(Mj)pCwpp^-fM?)r*W|)F  +<:iM?)S  (5-55) 
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Collecting  Equations  5-47,  5-54,  5-55,  5-45  and  5-46  in 
matrix  form  we  obtain  the  state  variable  representation  of 
the  lateral  directional  dynamics  model,  the  state  variables 
being   ^,    p,    r,   4>   and  t|>  . 
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(5-56) 


The  solution  that  follows  corresponds  to  the  given 
example  with  initial  conditions  and  control  inputs  indicated 
in   the   summary  of   specifications  tables. 

For  the  Dutch  Roll  approximation  the  following  model  may 
be   used 


(5-57) 


which  includes  all  the  parameters  of  the  determinant. 
Equation  5-30,  while  for  the  Roll  Subsistence  and  Spiral 
approximation 
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(5-58) 
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which    includes    all    the    parameters    of    the    determinant. 
Equation  5-36 . 


5.5       LATERAL-DIRECTIONAL   AERODYNAMIC   TRANSFER  FUNCTIONS 

Following    a    same    approach   as    in   the  preceeding   chapter, 

we    can    develope    lateral-directional    aerodynamix    transfer 

functions. 

For    the    lateral-directional    model,    the   input   variables 

were    the    aileron    control    deflection    (5)    and    the    rudder 

control   deflection  (?)  while   the   output   variables  were    f>,   p,    r, 

4?  and  ii> . 
Thus 


A^)  _ 


u 

0 


■? 


0 

s-Lp 

i 
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-5'%ri-v 


S-Mr 
0 
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-5/a 
o 
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-s 
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o 
-s 


SCs) 


5[   AS^   +35^  05^4-0)5  -hE] 


where 

1J, 


(5-59) 

(5-60) 
(5-61) 
(5-62) 
(5-63) 
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This  transfer  function  usually  takes  the  form 

Other    lateral -directional   transfer   functions  of    interest 
are 


Sfs)  As^'Bs'.-Cs^+T^s+e  (5-65) 

where 

A^  =    L5+^VjxvN5  (5-66) 

Q,   -  Y^'' CuNis-LpMr)  +U(;YvNr  +  Np)- M3  Ci^^YvLr)  (5-68) 

and  takes   the   form 

^(^^     '        As^  +  Bs'^Cs^Ds^-E  (5-69) 

Also 

Ks^     _       yir  5^4-  Br  S^-i-Cr  s  +^r 

where 

Ar  =    MS^^«/j^U  (5-71) 

%  =    Y?  CM^-I><yi,,l^)UsCNp-^%,Y.J-^J5a4Lp)  (5-72) 

Q   =   ^t  (L^Mp-Ni,jLp)  -LsY.Mp  ^  NsYvLp  (5-73) 

^r  =    %.  CuiM^-MsL(s)  •  (5-74) 

and   takes   the   following   forms 

rCs)     _        Ar  Cs^2r,)  C  S^^  ^^^UJr  S  +  CUr  ) 

^fs)  is-^fBs^Cs^^-DsA-E  (5-75) 
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or 

SC&5    "  ^s^^'Bs^Cs^.-^sfE  (5-76) 

It    is    further  desirable   to   include   the   transfer  function 
relating  the  acceleration  on  the  y-axis   to  the   control    input. 

The    acceleration   on   the    y-axis    is    given    from    the  original 
V-equation,    Equation  2-135,   which   for        =0  becomes 

ay  =    V   4- Dor  -00  (5-77) 
and   the   corresponding  transfer   function   is: 

a«  Ci)  Aa  s'*  +  'Ba  5^  +  Q  s%  T)«  s  -^  E^ 


where 

Aa   =    YsC^-^^/lxxIzz)  '              (5-79) 

-Be   =  -  Ys  CLp  f  Mr  +  ^^^j,^  Np  +  ^5izz  Lp )  ( 5-80 ) 

Ca   =  YsCN^+LpMr-MpLr  ^Vlz^lf)  -UaYv  C  Vi«U  ^  Ms)  (5-81) 

^a-M^,l^-kh   -W^^^P^'4'%x)]-f^-Y^[U(f^-Mp).N5(^^^,-^)]  (5-82) 

Ea  =    g  [Y*   (L^Nr  -N^Lr)  4-  Yv(M5Lr-LsWr)J  (5-83) 

and   takes   the   form 

5.6       LATERAL-DIRECTIONAL    EQUATIONS    IN    NON-DIMENSIONAL 
SYSTEMS 

In     the     non-dimensional     body     axes     system    the 
lateral -directional   equations  of   motion  are    [Ref.    2]: 
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-Crjfl  -sHg  -sGwp  SLc-Ci^r 

where   the    non-dimensional    massr    inertia  terras 

expressed  as 

|j  -  on/pSg 

?  being  now  the  half  of  the  wing  span  (6=  ^^2  ) 

L^  =  Iz2/pSe' 

In    the    non-dimensional    wind    axes    system 
directional   equations  are    [Ref.    3]: 


(3 

0        Cyj 

0) 

r 

X 

5 

(5-85) 

nertia  terras 

and   time   are 

the 


(5-86) 

(5-87) 
(5-88) 
(5-89) 
(5-90) 
lateral 


-2s+Cy^               2s                        C, 

^ 

0      Cj, 

PQ^               i^Qr            isQ,  -s^Ix 

V 

= 

Qs    <^, 

i 
5 

(5-91) 

pCtJ^             J^i^r"  5^2                   ^  ^^"^P 

cj) 

Cii|     Gjj. 

where    the    non-dimensional    mass,    inertia   terms 

3  and   time   are 

expressed  as 

fi   r    W  /pSb 

(5-92) 

(5-93) 

(5-94) 

x=   m^  /pSUo 

(5-95) 

^'  -  Vr 

(5-96) 
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TABLE  5-1 

SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE  LATERAL- 
DIRECTIONAL  EQUATIONS  DUE  TO  INITIAL  CONDITIONS 


A.  Variables  and  Initial  Conditions 
(S  =  0.05 

P  =  0.05 
r  =  0.075 
(t'  =  0.1 
-d?  =  0.2 
t  =  0.0 

B.  Stability  Derviatives   and  Constants 

Yy    =   -0.0829      L/s  =  -4.77  N^  =   3.55  Ixz/Ixx    =   0.0663 

¥5    =   7.656  Lp  =  -1.695  Np   =   -0.0025  Ixz/Izz   =    0.0370 

Lr  =  0.1776  Nr    =   -0.0957  g    =   32.714 

Lf  =  27.25  N|    =   -0.615  v„    =   660.0 

L5  =  0.666  Nj    =   -1.383 

C.  Special   Functions 

Y5*     ^   V5/U0 

L*^     -   L  L^  ^  C  I«/lx.)  N;j]  Id  ^    =  [  M^  ^  CT«/j^^  )  |_^]  /p 

Lp    =    [   Lp  ^-(IWl^x)  Np]  ID  ^p    =  [  Wp  ^(^Wlzz)  Lpl  1^ 

L"5    '    [    I5  tCl«/Txx)  Nf]/^  N|  ^    [Ki|  ^(^Wj2z)Lg]/i> 


D.      Control    Inputs 
5  =   0.0 
5=0.0 
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TABLE  5-1  (Continued) 


E,   Derivatives 


4^    =   ? 


? 


=     r 


F.      Outputs 

(ir    Pr    r,  4)  f  Tp    vs   time   at   interval    0.0125 
end  calculation  when   time    >   50.0 
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Figure  5-5 

/3 -Sideslip  Angle  vs  Time  for  Specifications  of  Table  5-1 
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Figure   5-6 
p-Roll  Rate   vs  Time   for   Specifications  of   Table   5-1 
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Figure  5-7 
r-Yaw-Rate   vs   Time   for   Specifications  of   Table   5-1 

220 


o 
o 


o 
o 


OJ 

o 
o 


o 
o 


o 
o 


uou 


aoi 


002 


003 


004 


FI    VS 


I  ME 


X-3ChLF=1 , GGE-Ol     UNITS     INCH= 
T-SCHLE-i  JJUE-Ul     UNITS     INCH, 

Figure  5-8 
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Figure  5-9 
T|^-Yaw  Angle  vs  Time   for   Specifications  of   Table   5-1 
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TABLE    5-2 

SPECIFICATIONS   FOR  COMPUTER   SOLUTION   OF   THE    LATERAL- 
DIRECTIONAL    EQUATIONS   DUE    TO    INITIAL    CONDITIONS   AND 

CONTROL    INPUTS 


A.  Variables   and   Initial   Conditions 
^  =    0.05 

p  =   0 .05 
r   =   0.075 
<P  =   0.1 
ID  =   0.2 
t  =   0.0 

B.  Stability  Derivatives   and  Constants 

Yv    =   -0.0829      Lfi   =   -4.77  N,a  =  3.55  Ixz/Ixx    =   0.0663 

Yj    =   7.656  Lp    =   -1.695  N?  =  -0.0025  Ixz/Izz   =    0.0370 

Lr    =   0.1776  Nr  =  -0.0957  y  =   32.174 

Ls    =   27.25  N?  =  -0.615  u„    =   660.0 

L5    =   0.666  N5  =  -1.383 

C.  Special   Functions 

Y5     =  Y5  IVo 

4      =    [l^  -CWlvx)    Np]/I)  N;   -    [Npt(Ix2/l^2)^]/3> 

I^     =    [  L^  +  CI«/i^O  ^f]  /!>  ^?^    f  ^i   fCTx2/lzz)L|]  /!> 


D.   Control  Inputs 

i  =  0.01  for  0<t<2.5        5  =  0.0  for  0<t<15.0 
f  =  -0.01  for  2.5<t<5.0      5  =  0.1  for  15.0<t<17.5 
f  =  0.0  for  t>5.0  5=  -0.1  for  17.5<t<20.0 

's  =  0.0  for  t>20.0 
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TABLE  5-2  (Continued) 

E.  Derivatives 

r   =   N(2  (2    f  N^p  +  Nr  r  f  Wj  ?  ^-  M$  ^ 

F,  Outputs 

(i  r    Pf    r,   4),-t>    vs   time   at   interval    0.0125 
end  calculation  when   time      >  so.o 
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p-Roll  Rate  vs  Time  for  Specifications  of  Table  5-2 
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Figure  5-10 
|S -Sideslip  Angle  vs  Time  for  Specifications  of  Table  5-2 
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CHAPTER  6 
STATE  VARIABLE  FEEDBACK  LAW 

6.1   INTRODDCTION 

The  automatic-flight  control  functions  of  an  aircraft  are 
provided  by  an  automatic  flight  control  system  which  is 
referred  to  as  autopilot.  Autopilot  statilizes  the  dynamic 
response  of  the  aircraft  providing  the  best  response  to 
disturbances  and  control  inputs.  Today,  almost  all  types  of 
aircrafts,  helicopters  and  missiles  are  equipped  with 
autopilots  with  complication  depending  on  the  particular 
application. 

The  equipment  comprising  the  autopilot  includes  sensors, 
controllers,  estimators,  computer  and  actuators.  To  damp  the 
vehicle  responses  to  disturbances  and  control  inputs,  the 
rates  of  the  angular  motion  in  roll,  pitch  and  yaw  are  sensed 
by  gyroscopic  instruments  and  fed  back  to  the  flight  control 
computer.  The  computer  combines  the  feedback  information 
with  the  desired  pilot  response  and  computes  commands  which 
are  applied  to  the  actuators  that  deflect  the  control 
surfaces. 

In  this  chapter  we  will  discuss  the  state  variable 
feedback  method   as   a  means   to  control   the   vehicle. 
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6.2      THE   FEEDBACK   CONCEPT 

The     negative     feedback     attempts    to    modify    the 

characteristic   equation  of   the   system   so   that  one   can   get  the 

specified  by   the    requirements   response. 

Consider    as    an    example    that    the    roll    motion    of    an 
aircraft   is   described  by  the   equation 

p   =  Lp^  4  L55  (6-1) 

The  open  loop  transfer  function  is 


u 


5-Lp  (6-2) 

Suppose    that  the   p-variable   is   continuously  measured  and 
fed   back   as    shown  in   Figure  6-1. 


Su 


=t( 


~\   eOs) 

Li 
s-  Lp 

^1?) 

J 

kf 

Figure  6-1 
Simple  Feedback  System  for  Roll  Stabilization 
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The  closed  loop  transfer  function  is 

J^    _        >^s/cs-Lp)         18 

^^'^  ~   >l^-l<pU/Cs-Lp)  '  s-nCtpU-Lp)  (6-3) 

It  is  clear  that  the  single  characteristic  root  changed 
from  L^  to  Kp-Lj-Lp,  Usually  L-p  is  negative  and  Ls  positive. 
So  the  system  will  have  a  faster  response  indicated  by  the 
greater  negative  root  of  the  characteristic  polynomial. 

In  automatic  control  systems  all  the  variables  of  the 
system  are  fed  back  with  corresponding  gains  K*s  chosen  to 
give  a  particular  desired  response. 

In  general,  the  system  is  described  in  matrix  form  as 

X  ^  A  X  i-  ?  u  (6-4) 

with  initial  condition  vector 

y  Co^  =  yo  (6-5) 

Applying  the  feedback  law 

where 

K   =    [Ki    K2    K3     ...    K„]  (6-7) 

is    the    feedback    gain   matrix    and    ua    the    desired    input,    one 
obtains  the   system 

X    =   CA-  BK  )  X  +  B  uj  (6-8) 

The  matrix  A-BK  now  determines  the  new  characteristic 
roots  of  the  system.  Note  that  the  characteristic  roots  of 
the   system.    Equation  6-4,    are   obtained  by  the   determinant 

Isl  -  Al    .  o  (6-9) 
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while  for  the  system.  Equation  6-8,  by  the  determinant 

I  si  -  CA  -  B^)  1  =  o  (6-10) 

Since  the  response  of  the  system  is  determined  by  the 
location  of  the  characteristic  roots,  any  desired  response 
will  specify  the  characteristic  polynomial  of  the  system 
which  if  equated  with  Equation  6-10  the  feedback  gains  are 
obtained. 

Feeding  back  all  the  variables  means  that  somehow  all  of 
them  must  be  measured  or  estimated.  Usually  not  all  of  the 
variables  are  measured  but  the  ones  that  are  not  measured  are 
estimated.  For  example,  the  pitch  angle  can  be  obtained  by 
integration  of    the  q-pitch   rate. 

e  --  J   ^  (6-11) 

Even  unstable  systems  can  be  stabilized  with  state 
variable   feedback  provided   that  the   system   is   controlable. 

Most  times  the  feedback  gains  are  selected  such  that  for 
the  system  to  meet  a  performance  criterion.  In  this  case 
optical   control   methods   are  used   to   select   the  gains. 

In  the  next  pages  a  computer  solution  is  shown  of  the 
given  example  for  the  longitudinal  and  lateral  directional 
case  using  the  state  variable  feedback  concept. 
Specifications  for  computer  solution  are  shown  in  Tables  6-1 
and   6-2. 
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TABLE    6-1 

SPECIFICATIONS   FOR  COMPUTER   SOLUTION   OF   THE    LONGITUDINAL 
EQUATIONS   DUE    TO    INITIAL   CONDITIONS,    CONTROL    INPUTS    AND 
STATE  VARIABLE    FEEDBACK   STABILIZATION 


A.  Variables   and   Initial   Conditions 
V  =   5.0    ft/sec 

w  =   2 .5    ft/sec 
q  =   0 .05    rad/sec 
o  =  0.07  5    rad 
h  =   10.0    ft 
t  =   0.0 

B.  Stability  Deriviatives   and  Constants 


Xy     =    -0.0097 

Zy    =   -0.0955 

My    =    0.0 

g 

=  32.174 

Xw    =   0.0016 

Zv.   =   -1.43 

Mw    =   -0.0235 

Vo 

=  660.0 

Xs    =    0.0 

Zi    =   -69.8 

Mw    =   -0.0013 
M.,,    =   -1.92 
Ms    =    -26.10 

C.      Special    Functions 


D.      Control    Inputs 

=   0.Oi+   f   for    0<t<10.0 

=  -0.01+   f   for   10.0<t<20.0 

=   0.0   +   f   for    t>20.0 

f    =    -0,005   u    -O.0OO5  vw    +  0.022  q 
-I-  aOC75  9    +  o.ooooooi  k 


E.      Derivatives 


=     UZo 


F.      Outputs 

V ,   w,   q,  9  f    h  vs   time   of    interval    0.0625 
end   calculation  when   time    >   250.0 
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q-Pitch  Rate  vs  Time  for  Specifications  of  Table  6-1 
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Figure   6-5 
6-Pitch  Angle  vs  Time   for   Specifications  of   Table   6-1 
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X-SCflLE=5= OOE+Ol    UNITS    INCH, 
T-SCnLE=5= OOE+01    UNITS     INCH. 

Figure  6-6 
h-Height  vs  Time  for  Specifications  of  Table  6-1 
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TABLE    6-2 

SPECIFICATIONS   FOR   COMPUTER   SOLUTION   OF   THE    LATERAL- 
DIRECTIONAL    EQUATIONS   DUE    TO    INITIAL    CONDITIONS,    CONTROL 
INPUTS    AND   STATE  VARIABLE    FEEDBACK    STABILIZATION 


A.  Variables   and   Initial   Conditions 
(S  =   0.05    rad 

p  =   0.05    rad/sec 
r   =   0.075    rad/sec 
ctJ  =   0.1    rad 
V  =   & .2    rad 
t  =   0.0   sec 

B,  Stability  Derivatives   and  Constants 
Yv    =   -0.0829      L^    =   -4.77        N,3   =   3.55 


^5 


=   7.656  Lp  =  -1.695 

Lr  =  0.1776 

L5  =  27.25 

L5  =  0.666 


C,      Special   Functions 

Y5    =   Y5/U0 
D     =    >!-  C   I'xz  /Ixxiza) 
L*^    =   [  L/J  +C^«/lKx)N,gI  /P 

L%    =    [  Lr  ^CI^2/t,,)N-J/^ 
L5    --   [  L?    tCIxz:/j^^)  ^^]/J) 


Np    =   -0.0025 
Nr    =   -0.0957 
Nj    =   -0.615 
N5    =   -1.383 


Ixz/Ixx    =   0.0663 
Ixz/Izz   =    0.0370 
q  =    32.174 
y^    =66  0.0 


Ni*  =  [  Mp  ^  Clxz/izO  Lp]  /!> 

M?    ^    [    Nr  +Clv2/j„)    LrJ/I> 

n||  ^  [  >^g  f  CI^/Jz?)  Lg]/J) 


D.   Control  Inputs 

?  =  0.01  +  If   for  0<t<2.5      5=  0.0  +  ?f  for  0<t<15.0 

i  =  -0.01  +  !f  for  2.5<t<5.0   t=  0.1  +  «f  for  15.0<t<17.5 

?  =  0.0  +  5f  for  t>5.0         ?=  -0.1  +  5f  for  17.5<t<20.0 

^=  0.0  +  '^f  for  t>20.0 

tf    =~0,0U8^     -0.222  p    -0,209  r    -0,3Ucp    -0.2(^7  H) 
'5  f    =       00515   $    +0.05<?6p    f0.4<^5K-   ^.  0.O9S2  4>  4.  o.  1 6-9©  ij? 
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TABLE  6-2  (Continued) 
E.   Derivatives 

4"    --  ? 


F.      Outputs 

l^r   Pf    r,  4?,  '4'   vs   time   at   interval    0,0125 
end  calculation  when   time    >   50,0 
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BETA    VS    TIME 


X-SCRLE=1 . OOE+Ol  UNITS  INCH, 
T-SCRLE=5= OOE-02  UNITS  INCH, 

Figure  6-7 
/3 -Sideslip  Angle  vs  Time  for  Specifications  of  Table  6-2 
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P       VS    TIME 


X-SCflLE=l .OOE+Ol  UNITS  INCH. 
T-SCRLE-K OOE-01  UNITS  INCH, 

Figure   6-8 
p-Roll  Rate  vs  Time   for   Specifications  of   Table   6-2 
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VS    TIME 


X-SCnLE=l .OOE+Ol  UNITS  INCH. 
T-SCflLE=5o OOE-02  UNITS  INCH, 

Figure  6-9 
r-Yaw  Rate  vs  Time  for  Specifications  of  Table  6-2 
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FI    VS    TIME 


X-SCflLE=1.00E+01  UNITS  INCH. 
T-SCflLE=l  , OOE-01  UNITS  INCH, 

Figure  6-10 
4^-Bank  Angle  vs  Time  for  Specifications  of  Table  6-2 
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PSI    VS    TIME 


X-SCRLE-KOOE  +  Ol    UNITS    INCH. 
T-SCflLE=l , OOE-01    UNITS    INCH. 

Figure  6-11 
^ -Yaw  Angle  vs  Time  for  Specifications  of  Table  6-2 
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CHAPTER  7 
SPECIALIZATION   TO   MISSILES 

7.1      INTRODUCTION 

The  preceding  chapters  have  been  oriented  primarily 
towards  the  dynamics  of  the  conventional  piloted  airplanes. 
The  equations  of  motion  for  the  missile  are  still  valid  since 
it   can  be   thought   as   an   automatically  controlled  aircraft. 

Missiles  are  found  in  different  geometry  considerations. 
Cruise  missiles  usually  have  no  rotational  symmetry,  i.e. 
rotation  through  360''  is  required  along  the  missile  axis  to 
leave  the  external  shape  unchanged.  Configurations  with  180° 
rotational  symmetry  are  seldom  encountered  while  120"  or 
higher-order    rotational   symmetry  are  more   typical. 

In  this  chapter  we  will  consider  the  cruciform  missile, 
i.e.,  with  90"  rotational  symmetry  to  take  the  advantage  of 
symmetries.  This  configuration  is  shown  in  Figure  7-1.  This 
missile  may  be  considered  as  typical  for  today's  large 
existing  missiles. 

In  the  next  sections  we  will  consider  the  simplification 
to  the  equations  of  motion  due  to  symmetry  and  the 
complication  of  concerning  the  coupling  effect  of 
longitudinal  and  lateral-directional  motions.  The  last 
happens  because  strong  rolling  moments  are  produced  when  the 
missile   is   simultaneously  pitched  and   yawed. 
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-x  yv 


Figure  7-1 
Cruciform  Missile  Configuration 

For  this  particular  model  we  will  assume  that  control 
inputs  are  coming  from  the  vertical  and  horizontal  fin 
deflectors  of  the  tail.  A  deflection  of  the  vertical  fins 
stands  for  a  rudder  deflection  while  a  deflection  of  the 
horizontal  fins  stands  for  an  elevator  deflection  .  These 
are  the  only  existing  control  inputs.  Roll  control  is 
obtained  by  differential  deflection  of  one  vertical  or 
horizontal    fin. 
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7.2      SYMMETRY  CONSIDERATIONS 

For  the  missile's  equations  of  motion,  the  reference  body 
axes  system  shown  in  Figure  7-1  will  be  used  as  in  the 
airplane  case,  and  straight  level  horizontal  flight  will  be 
assumed  as   the   nominal    flight   condition. 

Considering  the  existing  symmetries,  we  find  first  that 
in  addition  to  the  aircraft  symmetries  the  xy  plane  is  a 
plane  of  symmetry  and  that  makes  the  inertia  tensor  of  the 
missile  having  only  diagonal  terms  and  as  a  consequence  the 
body  axes  coinsides  with  the  principal  axes.  Further 
lyy  =    Izz. 

The  fact  that  the  missile  is  not  designed  as  an  airplane, 
gives  no  tendency  to  remain  in  the  same  roll  orientation. 
Vertical  and  horizontal  find  deflections  can  be  optimized  to 
give  zero  bank  angle  in  turning  maneuvers.  This  condition  is 
highly  desirable  in  missiles  as  well  as  in  airplanes  because 
the  lifting  forces  are  most  efficiently  generated  and 
lift-to-drag   ratio   is  maximum. 

However,  accidential  actuator  error  signals,  asymmetrical 
loading  of  the  control  and  lifting  surfaces  in  supersonic 
flight  and  atmospheric  disturbances  may  introduce  large  roll 
rates  and  large  roll  rates  cause  cross-coupling  between  the 
longitudinal  and  the  lateral -directional  motion.  This  makes 
the  equations  of  motion  not  separable  in  a  longitudinal  and  a 
lateral-directional  set  which  was  found  so  convenient  in  the 
aircraft   case. 
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We  therefore  find  the  equations  of  motion,  Equation  2-134 
and  2-135,  combined  in  one  set  of  equations  and  modified  as 
follows: 


ftr> 


^    -    ^    *9^  (7-1) 


%  (7-2) 


—  Z  =  w  -  u, 


i  ^    n  (7-3) 

l;:  *- "  ^  (7-5) 

—  N    =  r 

Iz2  (7-6) 

with  auxilliary  relations 

(7-7) 

(7-8) 

(7-9) 

Before    expanding    the    forces    and   moments    in    terms    of 

stability  derivatives  we  note   tht   due    to   symmetry  the  Y-force 

is   the   same   as   the   Z-force  which  makes  v=w  and 

Yv  =    Zw  (7-10) 

Y;,    =    Z;,  (7-11) 


? 

"^ 

4' 

% 

r 

9 

r 

- 

W 
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Also  the  N-yawing  moment  is  the  same  as  the  M-pitching  moment 

which  makes  r=-q  and 

Yr  =-Zq  (7-12) 

Nv  =  -Mw  (7-13) 

N;.  =  -MvO  (7-14) 

Nr  =  Mq  (7-15) 

because  lyy  =  Izz 

In  addition,  for  the  control  derivatives  one  has 

Yj  =  -Zn  (7-16) 

N5  =  Mn  (7-17) 

N5  =  M;,  (7-18) 

7.3      EQUATIONS   OF    MOTION   FOR   THE    NON-ROLLING    MISSILE 

If    we    expand    Equations    7-1    to   7-6    in   terms   of    stability 
derivatives  we  obtain 

V    =    Xu  u    +  Xv^,  w    -^e     +-  Vrt  0  (7-19) 

\t>i    -     Zy  L>    -t-  2i>»  vu    +  Voq  +   2y»  n  (7-20) 

ci     -      Mi»l;+^u)W»^jUw>^+M^c^^-jUr,y^  (7-21) 

or   substituting   Equation  7-20 

Cl    -    ('My  t    UviZu)  u    i-CMvw+/^'w2vw)lA<    4.fMa+WwDo)(^4-(Mn  +  )UiA)ZM))ir7-22) 

V  =    Yv  V    ~  Vor    +|<p   ^X^S  (7-25) 

9   =    Iv  V  -^  If  y    +  Ir  r    -^  If  ^  (7-24) 

f=     Mi/i/i-Mi/v-t.Np'p  +  Mrr  +  l^s^  (7-25) 

At    Equation    7-25    we*  would    like  to    keep    the    term   N;v 

because    it    is    the    same    in   magnitude  as    yi^i/.      Substituting 

Equation  7-23   for  v 
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r  r  (NvfMvYv)v  -i-Np?  ■t-rMr-Nv'U<,)r  ^^W;({;    -t-  (M2  + Nv Y^  )5  (7-26) 

If  the  missile  is  not  rolling  the  equation  is 
eliminated  as  well  as  the  p-stability  derivatives.  Also,  if 
the  v-velocity  perturbation  is  assumed  zero,  the  equation 
is  eliminated  as  well  as  the  v-stability  derivatives. 
Therefore,  one  obtains  the  following  model  for  the 
non-rolling  missile 

r   =r    CMv^Ni^Yv)v  +(;Mr'N^Uo)r  +  CN^  +  NvY^;)^ 
or  more  conveniently   in  matrix  form 


(7-2?) 
(7-28) 
(7-29) 
(7-50) 


w 

Zw 

V. 

0 

0 

w 

2rj 

0 

i 

Ht 

^l 

0 

0 

% 

4- 

< 

0 

n 

• 

V 

0 

0 

Yv 

-Vo 

V 

0 

Y5 

5 

r 

0 

0 

H* 

MT 

r 

0 

M;' 

(7-3i) 


This  model  may  be  augmented  with  the  auxiliary  equations 
and  transfer  functions  and  state  variable  feedback  can  be 
applied  for  proper  stabilization. 
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APPENDIX   A 
SOME    LAPLACE    TRANSFORMS 


tnlr^ 

K^:) 

-t>o 

i{m);  --  Ffe) 

1 

S(0 

1 

2 

uCO 

Cvnil 

slep) 

V^ 

S 

rCl)  or 

i    Cunil  rawp ) 

Vs^ 

4 

^:2 

Vs^ 

5 

c-^ 

VCsfa) 

6 

^e-^ 

ycsfa)^ 

7 

COS  (joi 

s/Cs^+co-') 

8 

sir)  wi 

tw/c  s^+u/j 

9 

0  e  ^^  cos  Ctci 

^•cp) 

^%.a..     *   ^■^Vs.c 

SOME    LAPLACE    TRANSFORM    THEOREMS 

Theorem   1:  i  [a  ^/O  ^b  f.Ci)  |     =     a  F^Cs)  f  b  F^Cs; 

Theorem   2:  i  {  ^I'a)]    =  ic  F(sj 

Theorem  3:  i  {  dW/dt  }    =     s  FCs) -  f[o) 

Theorem  4:  ^{d'm/di'}    =    %^  H^)  -    s  f  Co)   ~  <^^Mi\a 

Theorem  5:  .^  {  ciVCOM']     =    s' F(s)  -  s^-f(c)    -s^/dtl, 

Theorems:  i  |   f^*  ^(:T)dr  |      =     ^^^Vs 


cl'^/diN. 
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